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and,  in  some  cases,  FORTRAN  programs  are  described  for  both  deterministic 
and  statistical  analyses  of  the  inherent  errors  as  calculated  with  the 
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ABSTRACT 

This  report  describes  research  on  redundancy  management  in  digital  flight 
control  systems.  The  emphasis  is  on  the  properties,  techniques,  and  re- 
quirements associated  with  the  operations  of  monitoring  and  voting  and 
their  effects  on  the  closed  loop  system  operation  when  asynchronous 
sampling  is  used.  These  topics  are  among  those  being  studied  in  the 
real-time  simulation  facilities  of  the  Air  Force  Flight  Dynamics 
Laboratory's  Digital  Avionics  Information  System  (DAIS). 

Part  I is  concerned  primarily  with  the  monitoring  operation  for  quad- 
redundant  input  signals,  and  utilizes  assumptions  compatible  with  the 
DAIS  program.  The  redundant  signals  are  assumed  to  change  one-at-a-time 
and  the  monitor  uses  cross-channel  comparisons  and  a binary  test  for 
deciding  whether  two  signals  are  within  tolerance  of  each  other  as  the 
basis  for  its  fault-detection  algorithm.  All  possible  relationships 
among  quad-  or  tri -redundant  signals  subjected  to  such  comparisons  are 
tabulated  and  grouped  into  "patterns".  Patterns  are  in  turn  used  to 
deduce  "Keys",  which  are  useful  for  understanding  the  way  the  signal 
relationships  change  with  time.  A basic  algorithm  for  monitoring  based 
on  the  above  characterization  is  described,  tested,  and  compared  briefly 
with  other  currently  used  algorithms. 

Part  II  presents  three  extensions  to  a previously  reported  model 
for  closed  loop  flight  control  systems  that  have  dual-redundant, 
asynchronous  digital  controllers.  The  original  model  had  the  same  sample 
rate  for  each  controller  and  a fixed  time  skew,  or  offset  between  their 
respective  sample  times.  The  first  extension,  the  Multirate  Model,  allows 
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for  separate  sampling  rates  for  the  external  inputs  (pilot  input,  wind 
gusts,  etc.)  and  the  digital  controllers.  The  next  extension,  the  Delay 
Model,  allows  for  computational  delays  due  to  the  time  required  for 
data  conversions  and  control-output  computations.  The  third  extension, 
the  Output-Averaging  Model,  provides  for  averaging  the  control  outputs 

produced  by  each  of  the  redundant  controllers,  rather  than  always  selecting 

l 1 

the  output  of  the  same  controller,  which  is  the  scheme  followed  in  all 
the  other  models.  Equations  and,  in  some  cases,  FORTRAN  programs  are 
described  for  both  deterministic  and  statistical  analyses  of  the  inherent 
errors  as  calculated  with  the  models. 


I 

I * 

ft 

I " 


i. 


1-1 


• Part  I 

Characterization  and  Monitoring  of  Redundant  Signals 
1.0  INTRODUCTION  AND  SUMMARY 

The  inputs  to  redundant  digital  flight  control  systems  may  themselves 
be  redundant.1  In  such  cases  two  operations  may  be  performed  on  each 
set  of  these  redundant  input  signals.  First,  a monitoring  algorithm 
separates  the  faulted  from  the  unfaulted  signals.  Second,  a voting 
algorithm  produces  a single  best  estimate  of  the  true  value  of  the 
signal,  based  on  the  signals  that  the  monitoring  algorithm  declares 
as  unfaulted.  This  report  (Part  I)  is  concerned  primarily  with  the 
monitoring  operation. 

The  assumptions  made  herein  are  compatible  with  (but  not  neces- 
sarily identical  to  or  limited  to)  those  that  constrain  the  real-time 
simulation  facilities  of  the  Air  Force  Flight  Dynamics  Laboratory's 
Digital  Avionics  Information  System  (DAIS)  advanced  development 
program.  The  DAIS  program  was  established  jointly  by  two  Air  Force 
Laboratories  as  an  approach  for  reducing  the  spiraling  life-cycle 
costs  of  avionics.  The  overall  objectives  of  the  program  are  to 
develop  and  evaluate  a set  of  standardized  hardware  and  software 
modules  (core  elements)  that  can  be  configured  for  use  in  a wide 
variety  of  aircraft  types  and  missions.  The  modules  include  computers, 
multiplex  data  bus  hardware,  controls  and  displays  hardware,  and 
software. 

The  DAIS  hardware  modules  will  be  part  of  a Flight  Engineering  Facility, 
which  will  be  used  for  real-time  simulation  and  evaluation  of  digital  flight 
control  system  performance,  interactions  between  flight  control  and  avionics, 
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and  the  pilot  Interface  with  the  system.  The  basic  DAIS  flight 

control  system  is  quad  redundant.  Each  of  the  digital  flight  control 

processors  is  programmed  in  an  identical  manner  and  operates  on 

redundant,  sampled  Input  data,  but  the  sampling  is  asynchronous, 

which  means  that  the  common,  preprogrammed  sampling  rates  may  vary 

among  the  identical  processors  by  small  percentages.  As  a result, 

there  are  time  skews  in  the  input  data  sainpl  1ng--approaching  a full 

sample  period  in  the  worst  case--and  the  processors'  computations 

4 

are  close  to  each  other  but  not  bit-for-bit  identical. 

This  report  concentrates  on  the  monitoring  operation  for  quad- 
redundant  input  signals.  The  redundant  signals  are  assumed  to 
change  one-at-a-time  and  the  monitor  uses  cross-channel  comparisons 
as  the  basis  for  its  fault-detection  algorithm. 

Section  2 below  describes  a means  for  characterizing  all  possible 
relationships  among  quad-  or  tri -redundant  signals  subjected  to  a 
cross-channel  monitoring  scheme  that  uses  a binary  test  for  deciding 
whether  two  signals  are  within  tolerance  of  each  other.  States  are 
defined  and  tabulated.  Then  patterns  are  deduced  and  grouped  into 
Keys,  which  are  useful  for  understanding  the  way  the  state  changes 
with  time. 

Section  3 describes  a basic  algorithm  for  the  monitoring  operation. 
The  algorithm  is  based  on  the  signal -pattern  characterization  of  the 
previous  section.  Detailed  flowcharts  and  explanations  are  given 
for  the  algorithm  so  that  it  may  be  easily  implemented  in  software. 

A simpler  algorithm  is  also  mentioned  and  related  to  the  basic  algorithm. 
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Section  4 presents  an  evaluation  of  the  basic  algorithm,  based 
on  both  computer  simulation  studies  and  comparisons  with  other 
currently  used  algorithms. 

Section  5 contains  the  references. 


2.0  SIGNAL  CHARACTERIZATION 

2 . 1 Assumed  Sampling  and  Monitoring  System 

A typical  digital  flight  control  system  may  be  represented  as 
in  Figure  1.  When  redundancy  is  employed,  any  or  all  of  the  com- 
ponents depicted  may  be  replicated  or  multiplexed.  This  report 
assumes  quad-redundant  sensors  and  that  the  flight-control  processor 
implements,  in  software,  a monitoring  algorithm  to  isolate  the 
"good"  (unfaulted)  sensor  inputs  from  the  "bad"  (faulted)  ones. 

The  monitoring  algorithm  to  be  developed  employs  input-input 
comparison  monitoring  and  additional  logic  operations.  (See  Figure  2.) 
For  comparison  monitoring,  two  signals,  say  xi  and  x^,  are  compared 
according  to  | x,  " xi  | ^ Tolerance 

The  comparator  output  (one  of  the  f^  produced  by  the  C's  in  Figure  2) 
for  x^  and  Xj  has  a value  of  0 if  the  inequality  is  satisfied  and  1 
otherwise.  The  monitor  logic  implements  an  algorithm  that  uses  both 
the  signal  values  and  the  comparator  outputs  to  eliminate  those 
signals  to  be  faulted;  in  general,  the  algorithm  also  requires 
storage  of  past  history. 

2.2  States  and  Patterns 

Figure  3 shows  a possible  relationship  among  the  redundant  Input 
signals  x^,  x^,  x^,  and  x^.  The  bars  in  the  figure  Indicate  the  allow- 
able tolerance  between  signals.  (Generally,  this  tolerance  is  fixed 
at  a value  determined  by  analysis  of  the  system  dynamics  and  sampling 
rates.)  In  this  instance  x^,  x^,  and  x^  are  within  tolerance  of  each 
other,  as  are  x^,  x?,  and  x^;  x^  and  x^  are  out-of- tolerance  with 
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each  other.  A simplified  representation  for  the  signal  relationships 


in  Figure  3 would  be: 


X3X1X2X4 


We  call  any  such  set  of  relationships  a state. 

When  only  comparison  monitors  are  used,  their  outputs  cannot 
be  associated  with  a unique  state,  but  could  be  produced  by  several 
possible  states.  For  our  examole  the  6 comoarators  produce  the 

outputs  fi2=f13=f14=f23=f24=0,  f34=*’  but  these  outputs  could  also 
be  produced  by  any  one  of  the  following  three  states: 


X3X2X1X4 


X4X1X2X3 


x4x2x1x3 


A complete  discussion  of  the  correspondence  between  the  comparator 
outputs  and  the  states  is  given  below. 

The  6 comoarators  have  2b  or  64  distinct  output  combinations ; 
however  7 combinations  do  not  occur.  Table  1 shows  a complete  listing 
of  all  the  combinations  that  are  possible,  along  with  the  orderings 
of  that  can  occur  with  each  combination.  There  are  thus 

57  output  combinations  and  336  states. 

An  examination  of  the  table  shows  that  only  a small  number  of 
patterns  occur,  provided  that  slqnal  identification  numbers  are  ignored. 
This  fact  suggests  the  categorization  given  in  Table  2.  Here,  the 
circles  could  be  any  of  the  four  inputs  and  the  bars  again  show 
which  of  the  inputs  are  within  tolerance  of  each  other.  At  fault 
level  1 the  pattern  indicates  a so-called  "soft  fault"  in  that  the 
lowest-valued  and  the  highest-valued  signals  are  out-of-tolerance 
with  each  other  but  in-tolerance  with  their  other  neighbors.  At  level 
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Double  Hard  Fault 


Double  Hard  Fault 


Double  Hard  Fault 


Table  2 Fourteen  Fault  Patterns 
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3A  the  pattern  indicates  a "hard  fault"  in  that  the  lowest-valued 
signal  is  out-of-tolerance  with  all  three  remaining  signals. 

The  use  of  patterns  reduces  the  number  of  possible  relationships 
from  336  states  to  14  patterns.  However,  it  omits  some  of  the  infor- 
mation supplied  by  the  comparison  monitors.  For  example,  the  pattern 
OOOOoccurs  anytime  there  is  a single  soft  fault,  but  this  corresponds 
to  any  of  6 combinations  of  outputs  from  the  comparison  monitors  and 
to  24  states. 

There  are  24  possible  ways  to  order  the  4 signals;  hence  for 
14  patterns,  we  again  arrive  at  24  x 14  or  336  states.  Table  3 shows 
these  24  orderings  of  the  signals  and  the  comparator  outputs  for  each 
of  the  14  patterns.  This  table  shows  the  correspondence  between  a 
given  ordering  of  signals  and  the  fault  level,  the  pattern,  and  the 
comparator  outputs . 

2.3  Keys 

The  state  may  change  each  time  one  or  more  of  the  signals  change. 
Likewise,  the  pattern  may  change,  as  we  discuss  next. 

Assume  that  the  set  of  four  signals  is  examined  by  the  monitoring 
algorithm  each  time  any  one  of  the  signals  is  updated.  (In  those 
situations  where  more  than  one  signal  could  change  at  the  same  time, 
we  could  still  process  all  four  after  each  signal  change.)  With  this 
restriction  it  is  possible  to  develop  groups  of  patterns,  or  Keys, 
which  consist  of  all  patterns  to  which  the  present  pattern  may  move 
after  a change  in  one  signal.  Note  that  the  Keys  do  not  allow  us 
to  trace  the  sequence  of  states  that  occur--only  the  patterns. 


Level , Patterns  and  Comparator  Outputs 
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Table  3 Signal  Ordering  and  Fault  Level,  Patterns,  and  Comparator  Outputs 


rH  rH  rH  iH  rH  iH  » — 4 rH  r—4  rH  r — • r ( r~4  iH  r — 4 t— 4 O O i“4  r-4  r— 4 *— 4 O O 

* —4  r— 4 *— 4 r-H  f— 4 « — 4 rH  r — i t — 4 r— 4 0 0 r-4  r— 4 • — 4 t — 4 « — 4 r— 4 r-4  «— 4 O 0 r — 4 r— 4 

rHrHrHrHrHrHrHrHOOrHrHrHtHOOrHrHrHrHrHrHrHr-4 
r— 4 r— 4 f — 4 * — 4 0 0 r-4  r— 4 r — 4 • 4 r—4  r — 4 r-4  i — 4 » — 4 • — 4 i — 4 i — 4 0 0 i — 4 r— 4 *— 4 r— 4 

r— 4 r—4  0 C-*'  r—4  r—4  r—4  r—4  r — 4 r—4  r—4  r—4  0 0 r—4  r—4  r—4  r—4  r — 4 t — 4 r — 4 r—4  r—4  r—4 

0 O r—4  r—4  r—4  r—4  0 0 r—4  r—4  r — 4 t— 4 r—4  r—4  r—4  r—4  r—4  r—4  r—4  t — 4 r — 4 r—4  r—4  r-4 


r—4  r—4  r—4  0 r—4  0 r — 4 r—4  r—4  0 r — 4 0 r—4  r—4  r—4  r—4  r — 4 r—4  r—4  « — 4 r — 4 r—4  t — 4 r — 4 
r—4  0 r—4  r—4  0 r—4  r—4  r—4  r—4  r—4  r—4  r—4  r—4  r—4  r—4  0 r—4  0 r—4  r—4  r—4  r — 4 r—4  r — 4 
l i r—4  0 r—4  r—4  r—4  r—4  r—4  r—4  r —4  r—4  r “4  r—4  r—4  r—4  r—4  r—4  r—4  r—4  r—4  r—4  0 r—4  0 
r—4  r—4  r—4  r—4  r—4  r—4  r—4  0 r—4  r—4  0 r—4  r—4  O r—4  r—4  0 r—4  r — 4 r—4  r—4  r—4  r—4  r—4 
r—4  r—4  t— 4 r—4  r—4  r—4  0 r—4  O r—4  r—4  r — * r—4  r—4  r—4  r—4  r—4  r—4  r — 4 0 r 4 r—4  CD  ' — 4 
H H H i — • r—4  t 4 r—4  r—4  i — tr— 4r-Hr-HOrHCDrHr— 4rHCD»HCDr-Hr— 4r— 4 


0 CD  r—4  r—4  r—4  r—4  C~D  CD  r—4  r—4  r — 4 r—4  r-4  « — 4 r-4  r—4  0 (~P  r—4  r—4  r—4  r—4  CD  CD 

r—4  r—4  CD  CD  r—4  r—4  r—4  r — 4 r—4  r—4  CD  CD  CD  CD  r—4  i — I r — 4 r—4  r —4  r—4  CD  CD  r—4  r—4 

r—4  i — 4 r—4  r—4  CD  O ' — * r—4  0 CD  »— 4 r—4  r—4  r—4  CD  CD  r — 4 r—4  0 CD  r—4  r — * r—4  r—4 

r—4  r—4  r—4  r—4  O O rH  r— 4 O O rH  r- «r-4t— * CD  O r—4  r—4  O O r—4  r—4  r—4  r—4 

r—4  r—4  CD  CD  r—4  r—4  t— 4 r — 4 r—4  r—4  CD  0 CD  CD  « — * r—4  r — 4 r — 4 r—4  r-4  CD  CD  r-4  r—4 

CD  CD  r—4  r—4  r—4  r—4  CD  CD  r—4  r—4  r—4  i— 4 r -4  r— "4  r—4  r-4  CD  CD  r—4  r-4  r—4  r—4  CD  CD 


r—4  r—4  •— 4 CD  r—4  CD  r—4  r—4  r—4  0 r—4  0 r—4  r—4  r—4  r—4  CD  CD  r—4  r—4  r—4  t — 4 CD  Q 

r—4  O rHr-40rHiHrHr-4rH  CD  O rH  r—4  r-4  O r—4  CD  rH  rH  CD  CD  r—4  r-4 

O r—4  O « — * r—4  »— 4 r—4  CD  O r-4  r—4  r—4  r—4  CD  O r-  4 r-4  «— 4 r- 4 r— 4 O r-4  O 

r— 4r-4r- 4r— *00*~ 40r-4r— 40r—4r— 4CD*~4r--4  0*--400»— 4r— 4r— 4»— 4 
r—4  r—4  CD  CD  *-H  r—4  CD  t— 4 O r-4  r—4  r—4  CD  CD  »~H  rH  r—4  r—4  r-  4 CD  *H  t — 1 CD  r—4 

O O r-4  r— 4 r—4  i — 4 CD  CD  r—4  r-4  r-4  r—4  CD  rH  CD  r — I r—4  r—4  CD  rH  CD  r—4  r-H  r-4 


CD  CD  r—4  CD  r—4  0 CD  0 r — 4 CD  r—4  0 r—4  r-4  r—4  r H r—4  r—4  r—4  r—4  r—4  r—4  r—4  r—4 

r—4  O O O O r-H  r—4  r—4  r—4  r—4  r—4  r—4  O O r—4  O r-4  CD  rH  r—4  r—4  r-H  rH  r—4 

CD  r-H  CD  r—4  O CD  r—4  r—4  r—4  r—4  r—4  r—4  r-H  r-4  r—4  rH  rH  rH  CD  O r—4  CD  i — 4 CD 

r— 4r-Hr—4r—4r— 4r— 4r— 40000r-4r-H000CDr-Hr— 4r-Hr—4r—4r— 4r-4 
r— H rH  r—4  r—4  r—4  r—4  CD  rH  CD  r—4  CD  CD  * — 4 r 4 r H r-4  r—4  r—4  r-H  CD  CD  CD  CD  r—4 

r—4  rH  rH  rH  r—4  r—4  r— 4 rH  rH  rH  rH  rH  CD  rH  CD  r-4  CD  CD  O rH  CD  r-4  CD  O 


'O  Kf  (\J  0 (\J  0 r«  C*D  r—4  CDJ^J  r4  C\J  rH  V\l  Ol  H (\J  H 

^ (\)  ^ (\)  ^ ^ 0 H 4^  rH  0 (\J  rj  rH  H (\|  (\|  H 0 H i\| 

i\i  c\i  m ro  it  HHnnrj  h h m m r4  h (\n\i  0 0 

rH  H rH  H H r4  (M  (\|  (\J  (M  l\l  (\J  0 0 0 0 0 0 M M M 


to  o 


There  are  7 basic  Keys,  as  shown  in  Figure  4.  Key  A is  one  of 
the  simplest  and  most  important.  For  this  Key  the  symbol  f indicates 
a signal  which,  once  it  is  updated,  can  then  occupy  any  other  position 
marked  by  £ and  thus  produce  a new  state  having  a pattern  within  the 
same  Key.  The  symbol  O indicates  a signal  that,  when  it  changes, 
produces  a state  with  a pattern  outside  of  the  Key.  For  example, 
when  the  lowest-valued  signal  in  the  pattern  #(JO#  of  Key  A is 
updated,  the  new  pattern  could  be  unchanged  or  it  could  be  any  of  the 
other  6 patterns  in  Key  A;  the  new  pattern  is  guaranteed  to  be  in 
Key  A,  provided  that  the  given  signal  was  the  one  updated.  However, 
if  any  of  the  signals  denoted  by  Q is  the  next  one  to  change, 
the  new  state  will  either  remain  the  same  or  change  to  a state 
having  a pattern  in  one  of  the  other  Keys--but  not  Key  A. 

There  is  some  overlap  between  Keys  C and  E,  and  also  between 
Keys  D and  F.  For  Keys  C and  E consider  the  pattern  OOOO  • If 
the  second-lowest  signal  moves  to  be  left  there  are  two  possible 
results:  ei ther  QQQQ  or^VfpvQQ  , which  are  in  different  Keys 
(C  and  E,  respectively).  Similarly,  for  the  pattern  000~0  in 
Keys  D and  F,  if  the  second-highest  signal  moves  to  the  right,  this 
could  lead  either  to  or  t0  . In  these  two  cases, 

which  next  pattern  occurs  depends  on  the  magnitude  of  the  hard  fault 
that  had  been  present  and  not  just  the  fact  that  the  signal  differences 
exceeded  the  specified  tolerance.  Thus,  although  Keys  C and  E 
(D  and  F)  are  distinct,  they  cannot  be  distinguished  with  only  the 
information  supplied  by  the  comparators. 
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Another  way  to  present  the  information  of  Figure  4 is  shown  in 
Table  4.  This  representation  clearly  shows  the  ambiguity  present 
at  fault  levels  3,  4,  and  5.  As  the  number  of  hard  faults  increases, 
the  inability  of  the  comparators  to  indicate  fault  magnitudes  causes 
additional  overlapping  of  Keys.  For  example,  strictly  speaking  we 
cannot  move  from  CMxH'X'l  (Level  5,  Key  E)  to  QQQQ  (Level  3, 
Key  E).  However,  we  will  not  enumerate  additional  Keys  because  they 
are  not  used  in  the  monitoring  algorithms  presented. 
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3.0  BAFIC  ALGORITHM  FOR  MONITORING 

3.1  Relationship  to  the  Signal  Characterization 

In  this  section  we  present  two  algorithms  for  comparison  monitor- 
ing based  on  the  characterization  of  the  previous  section.  In  the 
first  algorithm  the  basic  idea  is  to  consider  all  four  signals  to  be 
good  as  long  as  the  fault  level  is  0 or  1,  but  to  eliminate  one  signal 
at  fault  levels  2,  3,  or  4.  With  these  restrictions  fault  levels  5 
and  6 are  not  needed  because  their  patterns  are  not  obtainable.  Further- 
more, only  Keys  A and  B are  needed  since,  from  fault  levels  0 and  1, 
they  contain  the  only  possible  routes  to  the  other  states  at  fault 
levels  2,  3,  and  4,  at  which  levels  a signal  is  eliminated  and  the 
redundancy  is  reduced  from  4 to  3. 

Figure  5 shows  the  combined  Keys  A and  B for  the  first  algorithm. 
Once  a signal  is  eliminated,  the  tri-redundant  case  becomes  applicable 
(Figure  6).  In  the  tri-redundant  case  we  again  consider  the  signals 
as  good  for  fault  levels  0 and  1 and  eliminate  the  hard-faulted  signal  at 
levels  2 and  3. 

Assume  that  at  system  start-up  the  signal  pattern  is  at  fault 
level  0 or  1.  If  the  pattern  stays  at  these  levels,  no  signals  are 
eliminated.  If  the  pattern  moves  to  fault  level  3 or  4 then  the  com- 
parison monitor  should  declare  the  most  recently  updated  signal  as 
being  faulted.  If  the  Dattern  moves  to  fault  level  2 the  comparison 
monitor  will  still  select  one  signal  as  being  faulted  but  the  logic 
for  identifying  the  faulted  signal  is  more  complicated  and  will  be 


I 


treated  below  during  the  discussion  of  the  flowchart  for  the  algorithm. 

The  second  algorithm  is  simpler  than  the  first  and  will  be 
treated  only  briefly.  The  simplicity  arises  from  the  constraint  that 
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a signal  is  declared  faulted  and  eliminated  when  the  pattern  is  at 
any  fault  level  other  than  0.  Since  Key  A contains  the  pattern  at 
fault  level  0,  this  Key  is  the  only  one  needed.  See  Figure  7. 

Once  a signal  is  eliminated  the  tri -redundant  case  becomes  applicable 
(Figure  8).  Since  this  second  algorithm  is  considerably  simpler 
than  the  first  one,  it  will  not  be  covered  in  detail. 

3.2  Flowchart  and  Detailed  Description  of  the  Algorithm 

The  first  algorithm  can  be  implemented  as  shown  in  the  flowchart 
of  Figure  9.  Table  5 provides  descriptions  of  the  subroutines  used 
by  the  main  routine.  These  are 

QSSD:  Quad  Signal  Selection  Device 
TSSD:  Tri  Signal  Selection  Device 
BSSD:  Bi  Signal  Selection  Device 
QUAD:  Quad  Redundancy  Routine 
QAUX:  Quad  Auxiliary  Routine 
TACPT:  Tri  Acceptance  Routine 
TRI:  Tri  Redundancy  Routine 
BACP1:  Bi  Acceptance  Routine 
BI:  Bi  Redundancy  Routine 
BACP2:  Bi  Acceptance  Routine 

The  first  3 of  these  subroutines  could  utilize  any  of  a number  of 
algorithms  for  choosing  the  best  estimate  of  the  signal  value  based 
on  redundant  measurements;  when  the  algorithm  requires  past  values  of 
the  signals  there  must  be  a way  to  assure  that  these  values  are 
available,  even  though  faulted  signals  may  have  been  eliminated. 
Figures  10,  11,  and  12  provide  flowcharts  for  the  last  7 subroutines 
used  by  the  first  algorithm. 


Table  5 Subroutines  Used  by  the  First  Algorithm 
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QSSD:  Quad  Signal 
Selection  Device 

Devises  a single  output  based  on  four 
redundant  inputs,  via  median  select,  mean 
select,  weighted  average,  random  select,  etc. 

TSSD:  Tri  Signal 
Selection  Device 

Same  as  QSSD  but  using  three  redundant 
inputs 

BSSD:  Bi  Signal 

Selection  Device 

Same  as  BSSD  but  using  two  redundant  inputs 

QUAD:  Quad 

Redundancy  Routine 

Updates  a signal  and  determines  whether  it 
is  definitely  faulted,  definitely  unfaulted, 
or  part  of  a pattern  to  be  analyzed  by  QAUX 

QAUX : Quad 

Auxiliary  Routine 

Used  in  conjunction  with  QUAD  and  TACPT  to 
determine  whether  a signal  other  than  the 
most  recently  updated  signal  has  faulted 

TACPT:  Tri  Acceptance 
Routine 

Determines  whether  a faulted  signal  now  has 
a value  that  would  qualify  it  as  an  unfaulted 
signal.  If  so  a counter  is  incremented  and 
once  it  reaches  a preset  limit  the  signal  is 
re-accepted  as  unfaulted  and  if  necessary 

QAUX  is  used  to  determine  the  status  of  the 
remaining  signals 

TRI:  Tri  Redundancy 
Routine 

Updates  one  signal  and  determines  whether  a 
fault  has  occurred.  If  the  updated  signal 
was  previously  declared  faulted,  then  control 
is  transferred  to  TACPT. 

BI:  Bi  Redundancy 

Routine 

Updates  one  signal  and  determines  whether  a 
fault  has  occurred.  If  the  updated  signal 
was  the  first  one  previously  declared  faulted, 
then  control  is  transferred  to  BACP1;  if  the 
second,  BACP2 . 

BACP1:  Bi  Acceptance 
Routine  1 

Determines  whether  a faulted  signal  now  has 
a value  that  would  qualify  it  as  an  un- 
faulted signal.  If  so  a counter  is  incremented 
and  once  it  reaches  a present  limit  the 
signal  is  re-accepted  as  unfaulted;  otherwise 
the  counter  is  reset. 

BACP2:  Bi  Acceptance 
Routine  2 

Similar  to  BACP1  except  this  subroutine  is 
for  the  second  failed  signal. 

UPDATE  SIGNAL  1 


NF1  > NACP1 


NF1  • NF? 
NACP1  » 
NACP? 
NF?  * 0 


NACP?  ♦ ) 


NF?  > NACP? 


RETURN 


RETURN 


4 ‘ L 
L?  ‘ L 
Li ’F 
NF?  • 

GOOD  • 

? 

1 

? 

0 

3 

PI  ■ 

NACP 1 ♦ 1 

F 

GOOD 

• ? 

Mf 
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The  first  step  in  the  flowchart  is  to  initialize  Lj,  L2.  L3,  and 
L^.  is  maintained  as  the  signal  number* of  the  most  recently  updated 
non-faulted  signal,  l_2  holds  the  signal  number  for  the  next  most 
recently  updated  one,  etc.  In  the  event  of  faults  the  first  faulted- 
signal  number  is  maintained  in  and  Fj,  and  the  second  In  Lg  and  F^. 
Since  the  system  may  be  asynchronous  the  sequence  of  the  signals  being 
updated  may  change  with  time.  This  fact  is  accounted  for,  as  the 
signal  numbers  are  shuffled  in  Lj  through  L^.  The  other  variables  that 
are  initialized  are  defined  in  Table  6,  which  contains  a complete 
listing  of  all  variables  used  in  the  flowcharts. 

A good  deal  of  the  processing  shown  in  the  flowchart  is  concerned 
with  startup.  The  main  loop  comprises  the  lower  half  of  the  flowchart 
and  will  be  discussed  first.  Afterwards,  startup  will  be  explained. 

For  the  main  loop  the  first  check  is  to  determine  whether  there 
are  4 unfaulted  channels,  i.e.,  whether  GOOD  = 4.  If  so,  QSSD 
selects  the  output  from  among  4 unfaulted  signals  and  control  passes 
to  QUAD. 

The  first  operations  in  QUAD  (Figure  10)  is  to  set  counters 
Nj  through  N4  to  0,  identify  the  signal  update,  and  reorder  the  contents 
of  Lj  through  to  reflect  the  latest  sequence  ordering.  Next, 

3 input-input  comparisons  are  performed,  using  the  most  recent  signal, 
x(Lj).  After  each  comparison  the  corresponding  counters  Nj  through 
are  incremented  whenever  the  corresponding  out-of-tolerance  condition 
is  found.  For  example,  if 

| xUj)  - x(L3)  > T 

then  Nj  and  N3  are  Incremented. 


*The  signals  may  be  thought  of  as  being  numbered  1 through  4. 


Variable 


Defini tion 


1 

1 

4 

Signal  number  of  the  most  recently  updated  non-faulted 
signal 

I 

( 

L2 

Signal  number  of  the  next  most  recently  updated  non- 
faulted  signal 

1 

i 

L3 

Signal  number  of  the  third  most  recently  updated  non- 
faulted  signal,  or  the  signal  number  of  the  second 
faulted  signal 

Id 

** 

• • 

L4 

Signal  number  of  the  fourth  most  recently  updated  non- 
faulted  signal,  or  the  signal  number  of  the  first 
faulted  signal 

- 

F1 

Same  as  for  G00D<  4 

F2 

Same  as  for  GOOD  < 3 

F 

1. 

NF 1 

Counter  for  the  number  of  times  faulted  signal  FI 
qualified  as  unfaulted 

f 

NF2 

Same  as  NF1  except  for  faulted  signal  F2 

1. 

r 

NACP1 

Maximum  count  for  NF1  before  faulted  signal  is 
re-accepted 

| 

NACP2 

Same  as  NF2  except  for  NF2 

f 

L 

K 

Counter  used  in  start-up 

GOOD 

Specifies  the  number  of  unfaulted  signals 

i. 

STRT 

Flag 

NSOFT 

Value  of  NACP1  or  NACP2  for  a soft  fault 

1 

NHAR1 

Value  of  NACP1  or  NACP2  for  a hard  fault 

1 

NHAR2 

Value  of  NACP1  or  NACP2  for  a hard  fault 

i 

N1 

Comparison-error  counter  for  x(L^) 

1. 

N2 

Comparison-error  counter  for  x(L2) 

f 

N3 

Comparison-error  counter  for  x(L^) 

t 

N4 

Comparison-error  counter  for  x(L^) 

* Table  6 Variables  Used  in  the  Flowcharts 

c 
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Next,  Nj  is  tested  for  a fault  in  the  most  recently  updated 
signal.  To  see  how  this  works  consider  the  examples  shown  in  Figure  13, 
in  which  the  signal  being  updated  is  and  the  values  of  N^,  N^, 


N3,  and  are  given  after  the  3 comoarisonsshown  in  Figure  10  are  made.. 
Figure  14  shows  the  counter  values,  after  3 comparisons,  for  the  most 
recently  updated  signals  (Key  A or  Key  B).  The  values  below  £ and 
correspond  to  and  those  below  the  unmarked  signals  are  for  whatever 
counter  (N£,  N^,  or  N^)  corresponds  to  that  signal. 

Examination  of  Figure  14  shows  that: 

(1)  When  Nj  is  0,  all  signals  can  be  declared  as  unfaulted 

(2)  When  Nj  is  2 or  3,  the  most  recently  updated  signal  can 
be  declared  as  faulted 

(3)  When  Nj  is  1,  either  all  signals  are  unfaulted  (Key  A, 
fault  level  1)  or  we  have  a pattern  where  a fault  may 
or  may  not  exist  (Key  B,  fault  level  2). 

If  the  third  case  is  present  (Nj  = 1)  then  all  6 comparisons  are  needed 
and  the  subroutine  QAUX  (Figure  10)  is  used.  Figure  15  shows  the 
counter  values  after  6 comparisons  are  made. 

Once  a fault  is  detected  in  QAUX,  F^  is  set  equal  to  the  signal 
number  of  the  faulted  signal  and  then  the  contents  of  through 
are  changed  so  that  now  corresponds  to  the  first  failed  channel; 
then  GOOD  is  set  to  3 so  that  the  next  time  signal  selection  is  per- 
formed TSSD,  rather  than  QSSD,  will  be  used. 

The  operation  of  TRI  is  similar  to  that  for  QUAD.  First,  counters 
Nj  through  N3  are  reset  and  a signal  is  updated.  If  that  signal  is 
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*QAUX  is  used  any  time  N,  ^ 1 but  only  in  the  case  N^l  are  more 
than  3 comparisons  needed. 
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Figure  14  Counter  Values  After  3 Comoarisons  with  the 
Most  Recently  Uodated  Siqnal 


the  one  that  was  previously  declared  as  faulted,  then  control  goes 
to  TACPT  to  test  whether  that  signal  can  now  be  re-accepted  as 
unfaulted.  If  the  updated  signal  was  not  previously  declared  as 
faulted,  then  the  values  of  Lj,  l_2,  and  l_3  are  re-arranged  and  the 
new  value  is  tested  for  a failure.  In  the  tri -redundant  case  a 
failure  occurs  if  any  counter  reaches  2.  (See  Figure  16.) 

Note  that  it  is  possible  for  all  signals  to  fail.  If  only  1 
failure  occurs,  and  become  the  signal  number  of  the  second 
failed  signal . 

If  TACPT  is  entered,  the  values  of  Lj,  L^,  and  L3  are  not  re- 
ordered. Instead  (to  save  time  later),  N^  through  N^  are  used  as  though 
the  re-ordering  were  done.  Thus,  in  TACPT  only,  Nj  corresponds  to 
L^,  N2  to  Lj,  to  l_2,  and  to  L^.  The  reason  for  this  is  that 
if  the  signal  corresponding  to  is  re-accepted,  then  the  signal 
numbers  will  be  re-ordered  and  the  counters  will  then  correspond 
exactly  (L^  to  Nj,  l_2  to  N2,  etc.);  then  if  QAUX  is  needed  to  test 
the  remaining  signals,  the  counters  will  not  require  re-ordering. 

Within  TACPT,  three  comparisons  are  made  with  x(L^)  and  the 
counter  Nj  is  tested  to  decide  whether  signal  is  still  faulted. 

If  it  is  not  faulted  (N^<  2),  the  counter  NF1  is  incremented  and 
compared  to  NACP1,  which  is  the  number  of  good  passes  necessary 
before  re-accepting  signal  L4  (Fj).  If  is  2 or  3,  signal  (F^) 
is  still  faulted  and  NF1  is  reset  to  0.  (NACP1  may  also  be  incremented 
to  increase  the  severity  of  the  requirement  for  re-acceptance.)  If  signal 
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is  re-accepted  (NF1  > NACP1),  then  the  signal  numbers  are  re-ordered 
to  show  signal  Fj  as  the  most  recent  update.  Also,  if  = 0 then 

we  are  assured  of  being  at  fault  level  0 or  1 and  we  can  return  with 
GOOD  = 4.  Otherwise,  Nj  = 1 and  we  must  enter  QAUX  to  further  test 
the  4 signals.  For  example,  if  eitherKO  ^or^Wi  is  the  pattern 

. O ' -M  r 

after  re-accepting  signal,^  , then  will  equal  1 and  we  may  or 
may  not  have  a new  failure;  QAUX  will  make  this  determination. 

The  BI  subroutine  receives  an  updated  signal  and  identifies  it. 

If  it  is  a faulted  signal,  then  either  BACP1  or  BACP2  is  used  to 
determine  whether  the  faulted  signal  should  be  re-accepted.  Otherwise 
the  updated  signal  is  compared  with  the  other  good  signal  to  determine 
whether  it  is  within  tolerance.  BACP1  and  BACP2  operate  much  like 
TACPT.  If  the  first  failure,  Fj,  is  re-accepted  before  F^,  then  F£ 
and  its  corresponding  parameters  are  transferred  to  F^. 

Note  that  in  all  the  re-accepting  routines,  the  oarameters 
NACP1  and  NACP2  can  vary  under  program  control.  For  example,  NACP1 
is  set  to  NHAR1  for  hard  faults  or  to  NSOFT  for  a soft  fault.  For 
the  second  failure  NACP2  is  always  set  to  NHAR2  since  only  hard 
faults  occur  in  the  tri-redundant  case  (Figure  16).  NACP1  and  NACP2 
may  also  be  incremented  each  time  a faulted  signal  takes  on  a new, 
out-of-tolerance  value.  (This  is  noted  by  the  optional  blocks  in 
Figure  11  and  Figure  12.)  NACP1  and  NACP2  are  also  used  during  start-up; 
here  they  are  initialized  to  0 to  facilitate  a fast  re-acceptance  of 
the  signals. 

This  completes  the  discussion  of  the  main  loop  and  associated 
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subroutines  for  the  first  algorithm.  The  startup  sequence  is  now 
treated. 

After  the  variables  listed  in  the  first  block  of  Figure  9 are 
initialized,  a signal  update  is  received  and  the  values  of  L1  through 
are  adjusted.  Then  Fj  is  set  equal  to  and  F^  to  lg.  Next, 
the  counter  K is  incremented  and  tested.  Since  this  was  the  first 
signal  updated  (K  = 1)  we  loop  back  to  update  another  signal  in  the 
same  manner.  Once  two  signals  have  been  updated  and  their  signal 
numbers  sorted,  we  enter  the  decision  block  at  the  bottom  of  the  flowchart. 
Then,  if  these  two  most  recent  signals  are  within  tolerance,  GOOD 
is  set  to  2 and  we  enter  the  main  portion  of  the  routine  with  NACP1 
and  NACP2  both  equal  to  0.  (This  facilitates  a fast  re-acceptance 
of  the  signals.)  If  the  two  most  recent  signals  are  not  within 
tolerance  we  return  to  the  top  of  the  routine  to  update  another  signal 
and  try  again  to  find  two  signals  within  tolerance  of  each  other; 
however,  once  K becomes  equal  to  7 the  flag  STRT  is  reset  and  NACP1 
and  NACP2  are  made  equal  to  their  nominal  values  before  continuing. 

Now  if  GOOD  is  still  1 (or  if  GOOD  becomes  1 during  normal  operation) 
the  test  associated  with  the  right-most  decision  block  will  indicate 
the  FAIL  condition  and  the  program  will  attempt  to  restart. 

The  routines  for  the  second  algorithm  can  be  obtained  by  simplifying 
those  for  the  first  algorithm.  In  the  second  algorithm  both  QAUX  and 
the  lower  portion  of  TRI  can  be  eliminated,  as  shown  in  Figure  17. 

Note  that  there  are  two  added  parameters,  NS0F1  and  NS0F2,  since  now 


there  are  two  levels  of  soft  faults.  Otherwise,  the  required 
operations  are  about  the  same  as  those  for  the  first  algorithm. 

The  second  algorithm  has  increased  sensitivity  to  soft  faults. 
One  way  to  overcome  the  sensitivity  is  to  increase  the  prespecified 
tolerance,  but  this  may  not  be  desirable.  Simulation  could  be 
used  to  resolve  this  question. 


4.0  EVALUATION,  COMPARISONS,  AND  CONCLUSIONS 

4 . 1 Computer  Si mul a tion  of  the  Algorithm 

The  basic  Algorithm  of  the  previous  section  was  programmed  and 
tested  using  Monte  Carlo  simulation.  All  4 redundant  signals  were 
given  a unit-step  input  plus  Gaussian  noise  of  various  levels.  The 
Gaussian  noise  was  generated  from  two  random  numbers  y^  and  y2  via 
the  transformation 

noise  = a noise  cos(2 'Tr  V 1n(l-y2) 
where  the  probability  density  function  of  yj  and  y2  is 

r 

i 1 0 < a < 1 

Py(a)  = { 

1 0 otherwise 

The  results  of  the  simulation  are  shown  in  Figures  18  and  19.  In 
Figure  18  the  values  of  NACP1  and  NACP2  are  1,  so  that  a faulted 
signal  is  re-accepted  the  first  time  it  comes  back  within  tolerance. 

In  Figure  19  NACP1  and  NACP2  are  increased  to  10. 

The  Figures  display  the  percentage  of  times  (out  of  1000  updates) 
that  QUAD,  TRI,  and  BI  were  used,  as  a function  of  the  noise  level 
on  the  input  signals  and  for  a tolerance  level  of  1.  Note  that  the 
rolloff  slope  of  QUAD  past  the  noise  level  of  0.2  is  greater  for 
the  larger  values  of  NACP1  and  NACP2.  The  program  was  set  to  restart 
when  a failure  occurs'.  (See  Figure  9.)  The  first  restart  occurs 
for  slightly  greater  than  0.3,  as  indicated  on  the  plots.  The 

restart  itself  begins  to  fail  for  ^n0-j se ^ ^ . 5 , but  the  program  was 
set  up  to  continue  to  try  to  restart. 


Tolerance  = 1 
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FAILURE  OCCURS 


Figure  18  Simulation  Results  for  Low  Re-Accentance  Level 
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F-inurp  io  Simulation  Results 

for  High  Re-AcceDtance  Level 

These  simulation  results  give  a general  idea  of  how  the  algorithm 


behaves  in  the  presence  of  noise  and  test  its  validity.  The  algorithm 
was  tested  further  by  allowing  the  inputs  to  be  updated  in  a random 
order.  Figure  20  shows  the  output  noise  variance  versus  the  input 
noise  variance.  Lower-median  select  was  used  to  select  the  output. 

The  values  of  NACP1  and  NACP2  do  not  appear  to  have  much  effect  on 
the  performance. 

4 . 2 Comparisons  With  Other  Monitoring  Algorithms 

The  algorithm  developed  here  is  applicable  to  the  DAIS  re- 
dundancy configuration,  which  currently  uses  a monitoring  and 

3 

signal  selection  algorithm  developed  by  Honeywell.  The  Honeywell 
algorithm  performs  the  same  6 input-input  comparison  monitoring 
steps  that  are  a part  of  the  algorithm  described  in  this  report. 
However,  the  Honeywell  algorithm  requires  that  a signal  exhibit  a 
hard  fault  for  a prescribed  number  of  algorithm-execution  intervals 
before  that  signal  is  declared  faulted.  Except  for  those  declared 
faulted,  all  signals  participate  in  the  signal  selection  process, 
which  selects  the  lower  median  (for  3 or  4 good  signals)  or  the 
average  (for  2 signals).  This  participation  of  faulted  signals 
(while  their  counters  are  being  incremented  up  to  the  prescribed 
limit)  in  the  signal  selection  process  may  result  in  degraded 
performance. 

Consider  Figure  21,  which  shows  three  instances  where  failures 


i 


occur  in  a quad-redundant  system.  Here  it  is  assumed  that  a signal 
must  exhibit  6 hard  faults  before  it  is  declared  as  faulted.  Case  (a) 
shows  a single  fault  and  good  behavior  of  the  algorithm.  Case  (b) 
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FIGURE  21  THREE  CASES  SHOWING  SIGNAL  SELECTION 
USING  THE  HONEYWELL  ALGORITHM 
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shows  2 faults— the  second  occuring  before  the  first  was  recognized 
as  a fault;  here,  the  algorithm  selects  a faulted  signal  as  the 
output  but  only  for  a brief  time.  Case  (c)  shows  a situation  where 
a faulted  signal  again  becomes  selected  as  the  output  but  this  time 
for  an  unlimited  time.  It  appears  that  when  the  fault-counter 
limit  is  set  high,  there  is  a greater  chance  that  the  occurrence  of 
situations  like  Case  (c)  will  cause  the  wrong  output  to  be  selected. 

Note  that  the  Honeywell  algorithm  detects  all  single  hard  faults 
but  does  not  eliminate  the  out-of-tolerance  signal  at  the  same 
time  that  it  is  detected.  The  algorithm  recognizes  all  soft  faults 
and  a split  fault  as  "glitches"  but  does  not  keep  track  of  which 
signals  are  involved.  The  procedure  is  to  count  the  total  number 
of  glitches  that  occur,  regardless  of  the  specific  signal  pattern. 

No  direct  use  is  made  of  the  glitch  count  or  glitch  occurrence  by 
the  signal  selection  algorithm. 

The  operation  of  Honeywell  algorithm  can  be  expressed  in  terms 
of  the  fault  patterns  developed  in  this  report.  Assume  that  all  4 
signals  are  participating  in  the  signal-selection  process.  (Perhaps 
only  because  the  hard-fault  counters  have  not  yet  reached  their 
prescribed  limits.)  Under  this  assumption  Table  7 shows  where  hard 
faults  and  glitches  are  detected,  which  signal  is  selected  as  the 
output,  and  which  fault  patterns  are  not  dealt  with  effectively, 
that  is,  poor  signal  selection  could  possibly  take  place. 

5-7 

Broen  gives  the  performance  of  several  voters  and  voter-estimators, 
both  of  which  are  designed  to  filter  random  variations  in  individual 
signals  while  simultaneously  discriminating  against  a faulted  signal. 
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The  voting  techniques  use  weighting  factors  chosen  so  that  the 
effect  of  one  faulted  signal  out  of  3 or  4 measurements  will  be 
minimal  in  comparison  with  the  remaining  unfaulted  signals.  The 
computations  require  multiplication  and  division.  Split  faults 
produce  outputs  that  may  straddle  the  split  and  thus  be  distant 
from  each  of  the  two  groups  of  signals.  This  difficulty  could 
possibly  be  reduced  if  the  quad  system  were  to  be  reduced  to  a 
triplex  system  by  some  simple  technique  such  as  removing  the  most 
positive  signal,  removing  the  last  failed  signal,  or  * r*ing  the 
last  signal  with  the  lowest  weight  value,  etc.  Broen's  voters  could 
be  used  as  the  signal  selection  device  in  the  algorithm  described 
in  this  report. 

Broen's  voter-estimators  require  a state  variable  model  to 
calculate  the  best  output.^ Here,  the  basic  idea  is  to  numerically 
isolate  a faulted  signal  and  to  apply  a least  squares  estimator  to 
the  remaining  signals.  The  total  computations  involved  probably 
are  not  feasible  for  the  DAIS  hardware  but  may  be  feasible  where 
additional  hardware  is  available. 
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Part  II 

Models  and  Software  for  Closed  Loop  Operation  of  Redundant  Systems 
1.0  INTRODUCTION  AND  SUMMARY 

Reference  1 describes  a model  for  closed-loop  flight  control  systems 
that  have  dual -redundant,  asynchronous  digital  controllers.  In  this  model, 
which  we  designate  the  Basic  Model,  the  digital  controllers  have  the  same 
sample  rate  but  there  is  a fixed  time  skew,  or  offset  between  their  res- 
pective sample  times.  This  same  skewed  sampling  scheme  is  used  throughout 
this  report  to  simplify  the  modeling  process. 

Three  extensions  to  the  model  of  Reference  1 are  developed  here, 
namely: 

(1)  Multirate  Model 

(2)  Delay  Models 

(3)  Output-Averaging  Models 

The  Multirate  Model  allows  for  separate  sampling  rates  for  the  external  in- 
puts (pilot  input,  wind  gusts,  etc.)  and  the  digital  controllers.  The  Delay 
Models  allow  for  computational  delays  due  to  the  time  required  for  data  con- 
versions and  control -output  computations.  The  Output-Averaging  Model  pro- 
vides for  averaging  the  control  outputs  produced  by  each  of  the  redundant 
controllers,  rather  than  always  selecting  the  output  of  the  same  controller, 
which  is  the  scheme  followed  in  all  the  other  models. 

The  Basic  Model  and  the  three  extensions  are  described  in  separate 
sections  below.  The  description  of  the  basic  model  is  in  summary  form  but 
that  of  the  others  is  given  in  detail.  Simplified  examples  are  presented 
for  which  the  calculations  can  be  done  by  hand. 


i 


Application  of  the  model  to  the  study  of  control  systems  for  realistic 
aircraft  requires  the  use  of  the  computer.  A software  package  for  the  Basic 
Model  is  described  in  Reference  1.  Corresponding  packages  for  each  of  the 
other  3 models  is  described  in  separate  sections  of  this  report.  FORTRAN 
Program  listings  and  example  output  are  given  in  the  Appendices. 


■ 
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2.0  BASIC  MODEL 

The  model  described  below  is  labelled  the  basic  model  because  the 
assumptions,  techniques,  and  style  of  analysis  are  the  basis  for  the  other 
models  described  in  the  remaining  sections  of  this  report.  The  basic 
model  was  developed  and  studied  in  detail  in  Reference  1. 

In  this  section  the  system  configuration  is  described  for  a closed- loop, 
dual -redundant  system  with  the  voting  rule  that  the  same  channel  (channel  1) 
is  always  selected  for  the  output.  The  system  state  equations  are  given 
without  showing  the  details  of  their  derivation,  as  these  details  are  avail- 
able in  Reference  1.  The  covariance  analysis  developed  in  Reference  1 is 
summarized,  as  are  two  examples,  which  are  also  considered  in  subsequent 
sections. 

2.1  System  Configuration  and  Dynamic  Equations 

The  closed-loop  system  configuration  for  the  basic  model  appears  in 
Figure  1.  The  system  input  is  assumed  to  be  the  continuous  vector  w , which 
is  applied  to  a continuous  time  model  of  the  aircraft,  the  sensor,  and  the 
control -actuator  dynamics.  The  system  output  yp  is  sampled  by  each  of  two 
digital  controllers;  they  use  the  same  sample  period  T,  but  controller  #2 
has  a skew  t.  The  inherent  error  e is  a piecewise-constant  function  of  the 
controller  outputs.  The  output  of  Controller  #1  serves  as  the  piecewise- 
constant  input  to  the  aircraft,  along  with  the  external  input  vector  wp. 

The  aircraft,  sensor,  and  actuator  dynamics  (plot  dynamics),  as  well 


as  any  dynamics  associated  with  the  external  input  take  the  form 


x = Ax  + B,  u + B0  w 
P P P 'P  P 2p  p 


yP  ’ cpxp 


(2-D 


(2-2) 


SAMPLE  PERIOD 


T : 

T:  SKEW 


FIGURE  1 BLOCK  OIAGRAM  FOR  THE  BASIC  MOOEL 
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where 

Xp  = plant  state  vector  (npxl) 

Up  = plant  input  vector  ( nupx 1 ) 

Wp  = disturbance  input  vector  (nwpxl) 

yp  = plant  output  vector  ( nopx 1 ) 

A„  = plant  state  matrix  (n  xn  ) 

P P P 

= plant  control  input  matrix  ( npxnUp ) 

B~  = plant  external  input  matrix  (n  xn  ) 
tp  p wp 

Cp  = plant  output  matrix  (nopxnp) 

Controller  #1  satisfies  the  following  vector  difference  equations 

xcl^k+l^  = Fcxcl(V  + Gcucl^k^  ( 

ycl(tk)  = Hcxcl(tk)  + Ecucl(V  ( 

for  k = 0,  1,2,...  . Controller  #2  satisfies 


(2-3) 


(2-4) 


■ FcV<Vt>  + 

yc2(VT>  = hcxc2(Vt)  + Ecuc2(VT> 

for  k = 0,  1,2,...  . Here, 

x , = controller  #1  state  vector  (n  xl) 

v*  I C 

u„,  = controller  #1  control  input  vector  (n  xl) 
ci  op 

y , = controller  #1  output  vector  (n  xl) 

Cl  (ip 

xC£  = controller  #2  state  vector  ( ncx 1 ) 

u^  = controller  #2  control  input  vector  (n  xl) 

yC2  = controller  #2  output  vector  (nupxl) 

F.  = controller  state  matrix  (n  xn  1 
c c c 


(2-5) 


(2-6) 
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G„  = controller  control  input  matrix  (n  xn  ) 
c r ' c op 

Hc  - controller  output  matrix  (static)  (nupxnc) 

E„  = controller  output  matrix  (inputs)  (n  xn  ) 
c r ' up  op 

The  plant  equations  and  the  controller  equations  are  related  by  the 
following  requirements 

(1)  The  control  input  to  the  plant  is  the  output  of 
controller  #1 . 

(2)  The  plant  output  is  the  input  to  both  controller  #1 
and  controller  i2. 

In  equation  form, 


up(tk>  ■ *c!(tk> 

(2-7) 

ud(tk>  ' W 

(2-8) 

uc2<V'>  = WT> 

(2-9) 

The  combined  equations  for  closed  loop  opet  ition  take  the  form 


x(tk+] ) = F(T.T)x(tk)  + G(tk,tk+1,f,w  (t ) ) 


where  x(t. ) is  a combined  state  vector  consisting  of  the  plant  variable  x 

k p 

and  the  digital  controller  variables  xcl  and  xc?,  as 


"xp(tk> 

’ xp(tkd> 

x(tk)  - 

xcl(tk> 

’ X^k+1^ 

xcl(lk*l) 

(2-10) 

I 

.xc2(tk*T,_ 

. xc2(Wt). 

and 


F(T.x)  = 


^k+l’V  + ^k+l’Wp  tk+l’tk^Hc  0 

GcCp  Fc  0 

L 6ccp[*(VT«tl)  + *(tk*T*tk,EcCp3  ^W^k^c  Fc. 


(2-11) 


J 


G(tk,tk+1,T,w  (t))  = 


{ k+1<t(tk+l’v)B2pwp(v)dv 


GcCp  { k *(tk+T’v)B2pwp(v)dV 


(2-12) 


where  4>(t,v)  is  the  state  transition  matrix  and  for  constant  Ap  is  given  by 


4>(t,v)  = exp[Ap(t-v)  ] 


The  controller  output  equations  are 


ycl(V  = Hlx(tk) 


yC2(tk+T)  = W + EcVt  *(VT’v)B2pwp(v)dv 


(2-13) 


(2-14) 


(2-15) 


where 


H1  ■ C£cCp  Hc  0] 


(2-16) 


H2  ■ [EcCp[«(tk*r,tk)  * +(tk+r,tk)EcCp]  EcCp„(tk*T,tk)Hc]  (2-17) 


The  piecewise-constant  inherent  error  e(t)  is  written  as  two  expressions. 


eft(t)  and  eB(t),  as 

e^(t)  - ycl(tk)  - yc2(tk^) 


(2-18) 


= (HrH2)x(tk)-EcC J k *(tk+r,s)B 2 w (s)ds 
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for  tk+i  < t < t^-p  0 < t < T,  k=0,l,...,  and 


egCt) 


= *cl(W  ' *c2(Vt) 


k+1 

= (H1F-H2)x(tk)  + H-jJ  *(tk+1,s)B2pwp(s)ds 

lk 


VT 


Wp  * *(tk+T’s)VP(s)ds 

lk 


VT 


* EcCp  ^ *(tk+T’s)B2pWP(s)dS 
1 1, 


(2-19) 


for  tk+i  £ t < + r»  0 < T < T,  k = 0,1,...  . 


2.2  Covariance  Analysi s 

The  inherent  error  can  be  characterized  in  a statistical  manner.  The 


external  input  wp  is  assumed  to  be  a sample  function  from  a Gaussian  white 


noise  random  process  with  0 mean  and  to  be  independent  of  x(0).  The  covariance 
matrix  of  the  states  is  defined  as 


Px(k)  = E[x(tk)xT(tk)] 


(2-20) 


where  E is  the  expected-value  operator.  P (k)  can  be  found  by  solving 


Px(k+1)  = F(T,T)Px(k)F' (T.t)  + V(T.t) 


(2-21) 


where 


V(T,r)  = 


* v0(t) 

0 

V (T)cV 
ov  ' p c 

0 

0 

0 

(2-22) 

G.C  V ( t ) 
c p 0 

0 

6 c r (t)cV 

c p o pc 

r 


l 


T 


r. 
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and 


v0(t)  = ; ♦(t,s)B2pW  bJp  »T(t,s)ds 


r 


The  matrix  W is  the  input-disturbance  covariance  matrix  defined  by 


E[wp(t)wp(T)]  = W ff(t-T) 


(2-23) 


for  all  t > tQ  and  x > tQ,  where  f(t-x)  is  the  Dirac  delta  function. 


The  steady-state  covariance,  designated  Pxs$  is  found  by  solving  the 
equation. 

Pxss  " F(T.t)PxssF(T,t)  * V(T.r) 


(2-24) 


Let  pe/\(t)  and  PeB^^  be  the  covariances  of  e/\  and  e^,  respectively; 
that  is,  for  example. 


PeA(t)  = EEeA(t) 


Then 


PeA(0  - (H,-H2)Px(k)(HrH2)T  + EcCpV0(T)cjEcT 


(2-25) 


(2-26) 


for  tk+x<t<tk+1  k=0,l,...,  0<t<T,  and 


PeB^)  = CHiF(T,t)-H2]  Px(k)[H1F(T,x)-H2]T  + EcCpCVQ(T)  - VQ(x)]cV 

(2-27) 

for  t <tk+]+T»  k=0,l,...,  0<x<T. 

All  of  the  above  equations  for  the  basic  model  are  developed  in 
Reference  1. 
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2.3  Examples 


As  the  first  example  consider  the  system  of  Figure  2.  For  this  system 


let  W = a . Then  the  analysis  reveals  that 

W 


PeAss  * KV<’  + 


(2-28) 


for  t.+T  < t < t.  ,,  0 < T < T,  k=0,l,...,  and 


PeBss  * * Sr1! 


(2-29) 


for  t,  i < t < t.  , + T,  0 < t < T,  K = 0,1,...  . 


^eAss  and  ^eBss  are  P^^ed  figure  3 as  a function  of  t.  The  diagrams 
to  the  right  of  each  plot  show  the  times  corresponding  to  the  values  of  the 


controller  outputs  used  to  calculate  e^  and  eg. 


As  the  second  example  consider  the  system  of  Figure  4.  Again  let 


W = o . Then  the  analysis  reveals  that 

W 


p = k2to  2 r2-K(J-T)(l+KT)i 
eAss  * w l2-KT(1+KTT  J 


(2-30) 


for  t.+T  < t < t.  .,  K=0,1,...,  0 < t < T and 


p = k2(t-t)  o 2 r2-^  (1+J^ 
KeBss  K u °w  L2-KT( 1+KT) 


(2-31) 


for  t.+1  < t < t.+.  + -,  K*0,l,...,  0 < T < T. 


PeAss  and  PeBss  are  P^otted  in  f:igure  5 as  a function  of 


iit 


I 


2-TK(UTK) 


FOR  EXAMPLE  2 


FIGURE  5 P( 
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i 


j 


ms 


r 


Note  that  the  results  obtained  for  both  examples  agree  with  intuition. 
For  small  T,  the  variance  of  eA  is  small  and  that  of  eB  is  large.  For  large 
x the  complementary  situation  holds. 


r 
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3.0  MULTIRATE  MODEL 

There  are  two  features  that  distinguish  the  Multirate  Model  from  the 
Basic  Model,  namely, 

(1)  The  external  input  is  sampled  using  a sample  period  T. 

(2)  The  controllers  operate  at  a sample  period  ^ where  n 
is  a positive  integer. 

In  all  other  respects  the  components  of  the  Multirate  Model  are  the  same 
as  the  Basic  Model.  Thus,  in  the  development  below,  the  same  state  variables, 
control  variables,  matrices,  etc.,  are  used  as  in  Section  2. 

3. 1 System  Configuration  and  Dynamic  Equations 

The  system  configuration  for  the  Multirate  Model  is  shown  in  Figure  6. 
Note  that  the  external  input  wp  is  now  sampled  and  at  a slower  rate  than  the 
rate  at  which  the  control  system  operates. 

For  the  aircraft,  actuator,  and  sensor  dynamics  we  may  write 


x = Ax  + B,  u + B9  w 
P P p IP  P 2p  p 

y = c„x 

yp  P p 


(3-1) 

(3-2) 


for  which  the  solution  is 


t t 

xp(t)  = *(t,t0)xp(to)  + / *U,s)81pup(t)ds  + / *(t,s)B2pwp(t0)ds 

fco  to  (3'3^ 


Let  tQ=tk  and  t = tk  + ^ . Then  (3-3)  becomes 

*p(tk  * K>  = *((k  * H'*t:k>Vtk)  * *l(tk  * n*tk)up(tk) 

* M'k  * 


where 


tfc«I 

k n 


*i<‘k  ■ ■ ! *(tk  ’ ff.‘k>BiPds 


(3-4) 


(3-5) 


r- 


DIGITAL 

CONTROLLER  * t 


DIGITAL 


CONTROLLER  2 


S/H  : SAMPLE- AND  - HOLD 


T:  PILOT  -INPUT  SAMPLE  PERIOD  AND  tk*l-tk»T 
T/  n : DIGITAL  - CONTROLLER  SAMPLE  PERIOO  ( n A POSITIVE  INTEGER) 


FIGURE  6 BLOCK  DIAGRAM  FOR  THE  MULTIRATE  MODEL 


AIRCRAFT 

AND 

ACTUATOR 

DYNAMICS 

11-18 


*2(tk  + H’tk) 


/ ®Uk  + n’VB2pds 


(3-6) 


The  discrete-time  equations  for  controller  #1  are 


'■cl^k  * n*  ’ Fcxclltk)  * Ccucl(tk* 


J'clf'k’  = Hcxcl(tk>  * Ecucl(tk> 


for  k=0,l,...,  and  for  controller  #2, 


xc2(tr  + n = + G.U^(VT) 


(3-7) 


(3-8) 


c c2l  k 


c c2l  k 


yc2(tk+T)  = Hcxc2(tk+T)  + Ecuc2(Vt) 


(3-9' 


(3-10) 


for  k =0 , 1 , . . . . 

The  plant  (aircraft,  actuator,  and  sensor  dynamics)  and  the  controllers 


are  related  by  the  equations 

up(tk>  * *cl(tk> 

“cl^k*  * V V “ cpVlk' 

uc2(tk*T)  = y v° ■ sy  V'> 

Substituting  (3-8)  for  yc^  and  (3-12)  for  u^  in  (3-11)  yields 
up(tk^  = Hcxcl^tk)  + EcCpVV 


(3-11) 


(3-12) 


(3-13) 


(3-14) 


Then  substituting  (3-14)  into  (3-4)  gives 


I 


f 

f 

i 

I 

i 


I* 

L 


4 

I 

i 

i 


i 


I 
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xp<‘k  * S>  * *<‘k  * K’WV  * ♦l<‘k  * S-VNcWV  * Eccpxp' 

* *A  * K’WV  = [»<S>  * *l<‘k  + S'‘k>Eccp]xP(tk> 

* *l(*k  *i-lk>Hcxcl(tk>  * *2<lk  ♦ii-WV 

(3-16) 

This  equation  for  xp(tk  + has  no  terms  involving  u (tk). 

In  a similar  manner,  substituting  (3-12)  into  (3-7)  and  (3-8)  gives 

xcl<tk  +S>  ■ Fcxcl(tk>  * GcCpxP(tk>  <3-'7> 

and 

yci(tk> 5 Hcxd(tk>  * EcVp<‘k>  <3-,8> 

From  (3-3)  with  t = tk  and  t = tj+x  , we  have 

xp(tk+x)  = [*(t)  + VVT-tk)EcCp]  xp(tk)  + ♦1(tk+T.tk)Hcxcl(tk) 

+ 4»2(tk+T,tk)wp(tk)  (3-19) 

Then  substituting  (3-13)  and  (3-19)  into  (3-9)  and  (3-10)  gives 

xc2^k  + n + = Fcxc2^tk+T^  + GcCp^T)  + 'h  ^tk+T,tk^EcCp-'xp^tk^ 

+ GcCp  ^*1  ^tk+T*tk ^ HcXc  1 ^ tk ^ + GcCp'l'2^tk+T’tk^wp^tk) 

(3-20) 

■ V ' Hc'c2,‘k*T)  * EcV»<T>  * ♦ l(VT-VEcCpVtk) 

• ' '.♦t<VT-,k>Hc«cl(tk>  * EcCp*2'VT'tk>"p(tk>  (3-?|) 


fa 
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The  state  and  controller  equations  can  be  put  in  compact  form  by  writing 
them  in  terms  of  a combined  state  vector.  Let 


*(tk) 


yv 

xcl*tk) 

xc2*V'* 


and  x(tj(+ii) 


V‘k x 

«ci'H 


xc2(tk  + n + 


The  state  equations  become 


x(tk  + ^)  1 AF  x(tt)  + AG  wjtj 


1 3-221 


where 


*(J)  + i|),4)E  C 
vn'  H 1 n c P 


v|».  (— )H 
M vn'  c 


(3-23) 


[ GCCPC"(T)  +*1(t)ECCp]  GcVl(T,Hc  Fc 


*,(i) 

AG  = 0 


(3-24) 


GCCpV  " 


T T ?T 

For  every  transition  interval  {(tk,tk  + -),  (tk  + tk  + 

(tk  + ^n'UT*  VT)}  *i^k  + W’V*  vl,i(tk  + f lk  + 1] 'MV'-H  + (iVr)J) 

are  equal.  Then  let  tk*0  in  equation  (3-5)  we  have  ^(tk  + tk)  = ^(t  n). 
Similarly  ^(t^x,^),  <p?(tk+T/n,tk)  and  4i?(tk+T,tk)  can  write  in  ^(t), 
vjip ( T/n ) and  \p0 ( t ) respectively,  and  we  use  ^(T/n),  <Pj(t),  and  vj>p( T/n ) instead 
of  <p1(tk+T/n,tk),  ( tk+T » tk ) and  \|>?( tk+T/n . tk > in  (3-23)  and  (3-24). 
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As  discussed  above,  and  E^.C  ,F  .G„  and  H are  constant.  Then  F and  G 

c p c c c 

T T 2T 

hold  for  every  transition  interval  {(tk,tk  + , (tk  + -,tk  + — ) , .... 

(tk  +^~^-  . tk+T)}.  Therefore 


x(tk  + £T)  = AF  x(tk  + J)  + AG  wp(tk) 


Substituting  (3-22)  into  (3-25)  gives 


Simi larly. 


and  so  on. 


x(tk  + jp)  = AF2x(tk)  + ( AF + 1 ) AG  w (tk) 


x(tk  + ~)  = AF3x(tk)  + (AF2  + AF+1)AG  w (tk) 


We  can  write  the  general  equation  for  these  equations  by 


x(tk  + Jp)  = AFm  <x ( t k ) + E ( aF ) 1 AG  wp(tk) 
where  m=l ,2, . . . ,n. 


Now  let 


xp(tk+l) 

x^k+l^  = xcl(tk  + l) 

. Xc2^k+1+T] 


Then  from  (3-28),  for  m=n 


where 


x(tk+])  = AFn  x ( t k ) + (AF)1^  wp(tk) 


x(tk+i)  = F(T,t)  x(tk)  + G(T,t)  wp(tk) 


F(T,t)  » (AF)' 


(3-25) 


(3-26) 


(3-27) 


(3-28) 


(3-29) 


(3-30) 


(3-31) 


I 
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n-1 

6(T,t)  = Z (AF)1AG 
i=0 


(3-32) 


From  (3-18)  and  (3-21)  we  can  write 


ycl-tk-  r H1  X^V 


(3-33) 


yc2(tk+r)  = H2x(tk}  + EcCp^2(T)wp(tk) 


where 


H,  - [E?Cp  Hc  0] 


H2  ■ [EcCp(»(T)  ♦ +-,(T)EcCp)  EcCp*,(T)Hc  Hc] 
As  same  as  (3-25).  The  controller  output  equations  are 
ycl^k+T/,n^  = Hix(tk+Tyn^ 


(3-34) 


(3-35) 


(3-36) 


(3-37) 


yc2(tk+T/n+T)  = H2x(tk+T/n)  + EcCp^2(T^  wp(V 


(3-38) 


The  general  equation  of  the  controller  output  equations  are 


yci(tk  +r}  = H^x(t^ 


T^k  1 n 


(3-39) 


yc2(tk  + f + T)  - H2x(tk  + ?)  + hWT)  W 


(3-40) 


where  m=0,l,...,n. 

Substituting  (3-28)  into  (3-39)  and  (3-40)  gives. 


mX  m m-  I • 

ycl(tk  + r]  = HlE(AF)m  x(tk)  + Z (aF)’aG  wp(tk)  (3-41) 


f 
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T m 1 i 

yc2(tk  + ™ + T)  = H2(AF)  x(tk)  + [H2  Z (aF)  AG  + E C vP3(t ) ]wp( tk ) (3-42) 


where  m=l,2,...,n. 

It  is  necessary  to  derive  (3-37),  (3-38),  (3-39),  (3-40),  (3-41)  and 
(3-42)  because  in  this  model  the  inherent  error  is  defined  as 


(3-42) 


k = 0,1,... 

+■  4-  + t ✓ t x t 4-  — T / T 

k n — k n *•  m=0, 1 , . . . ,n-l  and  0<r<- 


eB«,(t)  = *cl(tk  + ^ -*c2(tk  +f+T) 


(3-43) 


k=0,l , . . . 


4.  , in*’ if*  4 4.  _l  in*"  i j x f 

t,  + T < t < t.  + T + t { - . 

k n - k n vm=0, 1 , . . . ,n- 


1 and  0<t<- 


Figure  7 shows  the  skewed  sampling  and  inherent  errors  of  the  multi- 
rate model.  Channel  1 produces  the  sampled  outputs  at  times  tk,  tk  + ..., 

t^+1,  for  k=0,l,...  and  channel  2 produces  the  sampled  outputs  at  times 

tk+T»  \ + ^ + tk+1+T* 

Substituting  (3-39)  and  (3-40)  into  (3-42)  and  (3-43)  gives 

<W‘>  = x(tk  «•  f)  - EcCp*?(T)wp(tk)  (3-44) 


eBm(t)  = (VF_H2)x(tk  + 21)  + [HlAG-EcC  ,2(t)]w  (tk)  (3-45) 


where  m=0, 1 , . . . ,n-l . 


•j 


CHANNEL  I 
t ycl  ) 

CHANNEL  2 
<»c2  ) 

‘Am  (tl  : Vcll^lT  ) -*c2  K*-5?  +T) 

FOR  thv5^  +T  < t < tk  ♦ ^ T k =0 , i , 

m=0,l,...,n-l,and  0<  r 

*Bm(tl:  »CJ  (V5?") 

FOR  tk+I!^iT  <t  <t  +S2ii-T+  k = 0,l,... 

m =0,1,  . . . ,n-l  and  0<  r< 

FIGURE  7 SKEWED  SAMPLING  AND  INHERENT  ERRORS 
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Substituting  (3-41)  and  (3-42)  into  (3-42)  and  (3-43)  gives  e^t)  and  eBm(t) 
in  terms  of  x(tk) 

m-1  . 

eto(t)  « (HrH?)AFx(tJ  + C(H1-H2)  E (AF)1AG-ErCn4i9(T)]wn(tl>)  (3-46) 


cV2'l,JV  kJ 


for  t,  + H!l+  T <t<t  + T{k  = ®»^»*** 

k n - k n 1 m = l,2,...,n-l  and  0<T<T/'n 


eBm(t)  = (H!AF-H2)(AF)  x(tk)  + C ( H] AF-H^ ) 


for  ^ + ^ T <t<t.  + T+t  {k  “ 

k n - k n = 1 7 . 


(3-47) 


m = l,2,...,n-l  and  0<x<T/n 


We  can  write  e^m( t ) and  e0m(t)  when  m=0  in  terms  of  x ( tk ) by  using  (3-44)  and 
(3-45) 

eM(‘)  - (H,-H2)  x(tk)  - EcCp  2(t)wp(tk)  (3-48) 

for  tk+T<t<tk  + ^ , k = 0,1,...  and  0<T<T/n 


eB0^  = (H]aF"H2)  x(tk)  - CH]AG-EcC  i)»2(t)]w  (tk) 


(3-49) 


f0r  \ + H -t<t;k  + n + T’  k=0’1’---  and  0<T<T/n. 

From  (3-42),  we  define  the  inherent  error  e^t)  where  m=0, 1 , . . . ,n-l  for 
{(tk+T,  tk  + J),  (tk  + J + t,  tk  + (tk  + ^ + t,  tk  + T)}  so 


let  us  define  EA?  the  average  error  of  e A in  period  (tk,tk+T), 


Ea  = n ,n=0 


(3-50) 


Substituting  (3-46)  and  (3-48)  into  (3-50)  gives 


. n-i  , n-i  m-i  . 

Efl  n E (HrHJ(AF)1 "x(tj  + ^ E [(HrH  ) E (AF/AG]w  (t.  ) 
A n m=0  2 * n m=l  1 i i=0  P 


' K m=0 
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H, -H_  n-1  _ H.-H,  n-1  m-1  . 

Ea  = E (AF)mx(tJ  + [-L-2(  E E (AF),)AB  - E„CJ,,(t) ]wn(t t)  (3-51) 


m=i  i=0 


'c  pv  1 J p'  k' 


Similarly  we  define  Eg;  the  average  error  of  eg  in  period  (t^.t^+T+T)  so 


Eb  = n m-n  eBm  ' 
m=0 


Substituting  (3-47)  and  (3-49)  into  (3-52)  gives 

1 n_1  m 1 n~1  m-1  i 

E.-  I (H,AF-H  )(AF)m  x(t.)  E [(H.AF-HJ  E (AF)  ABVftJ 

0 n m=0  ' £ n m=l  1 1 i=0  p 


H.4F-H,  n-1 

- 4_  v 


H.AF-H-  n-1  m-1 

1 1 v vi 


E (AF)  x(t.  ) + [—■■■--(  Z E ( AF ) 1 ) AG  + H,/ 
n m=0  K n m=i  i=0  1 


(3-52) 


- EcCp*2(T)]  wp(tk) 
n-1  m-1  ^ 

Note:  The  factor  E E (AF)  can  be  expressed  in  closed  form: 
m=l  i=0 


n — 1 m—  1 • n — 1 •.  /*c\^  n— 1 i / *r\^ 

E t (4F)’  . t • i . ML, 

m=l  i=0  m=l  1 v m=l  Ah  1 Ah 


n-1  n;]  (AF)m  1 _ n 1 - ( AF ) n 

]-AF  ' m=0  ^ ‘ ]-AF  " (T^F)7 


(3-53) 


f 
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3,2  Covariance  Analysis 
3.2.1  Covariance  of  the  States 

Let  the  input  wp(tk)  be  a Gaussian  white  noise  random  process  with  zero 
mean  which  is  independent  of  x ( 0 ) (Reference  2).  Let  {•}  indicate  the  ex- 
pected value;  then 

E[wp(tk)]  = 0 (3-54) 

E[x(t0)wpT(tk)]  = 0 (3-55) 

and  let  Wk  be  the  covariance  matrix  of  wp(tk).  Then 

E[wp(tk)«pT(tk)]  = Uk  (3-56) 

Let  us  define  the  covariance  matrix  of  the  states  as 

Px(k,m)  = E[x(tk  + Jl)  xT(tk  + f-n  (3-57) 

and 

Px(k,m+1)  = E[x(tk  + lH!±il-I)xT(tk  + (3-58) 

= E{[AFx(tk  + 21)  + AGwp(tk)][AFx(tk  + J)I)  + AGwp(tk)]T} 

= AF  Px(k,m)(AF)T  + AGwk(AG)T  (3-59) 

where  m=0,l,...,n-l. 

Substituting  (3-28)  into  (3-57)  gives  P (k,m)  in  term  of  P (k)  where 

A A 

Px(k)  is  the  covariance  matrix  of  x(tk)  and  then 

Px(k)  = E[x(tk)xT(tk)]  (3-60) 

Px(k)  = Px(k,0) 


and 
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Then  we  have 


m t m m—  1 • » ID“  1 y : 

P (k,m)  = (AF)mP  (k)(AFr)  + E (aF)1AGW,  (AG)1  Z (aF1)1 
x i =0  i=0 


(3-61) 


Let  us  define  the  covariance  matrix  of  the  states  P (k+1)  as 


Px(k+1)  = E[x(tk+1)  xT(tk+l )]  . 


(3-62) 


Substituting  (3-30)  into  (3-62)  gives 


Px(k+1)  = F(T,-)Px(k)FT(T,T)  + G(T,-)WkGT(T,T) 


(3-63) 


The  steady-state  covariance,  designated  pxss  is  found  by  solving  the  equation 


Pxss  ■ F<T,t)  Pxss(T-t)  + G<T-t>  “kGT(T-’> 


(3-64) 


3.2.2  Covariance  of  the  Errors 

The  covariances  of  e^t)  and  egm(t)  are  calculated  using  the  same  pro- 
cedure as  in  the  previous  developments.  Let  PpA  be  the  covariance  of  eA  then 


W>  ■ 


(3-65) 


For  m=0,  substituting  (3-48)  and  (3-60)  into  (3-65)  gives 


W1’  " <HrH2>Px<k»Hl-H2)T  * EcCpVT>“k  2T(t)CpTecT  ,3-66’ 


for  tk+T<t<tk  + - , k = 0,1,... 


For  m>l,  substituting  (3-46)  and  (3-60)  into  (3-65)  gives 
PeAm(t)  = (HrH2)(AF)mPx(k)[(HrH2)(AF)m]T 


m-1  : m-1  i t 

+ [ ( H , -H  ) E (aFTaG-E  C *,(T)]W.[(HrH.)  E (aF)’aG-E  C _<M-)]' 
l 2 i=0  c p 2 K £ 1*0  C P 


(3-67) 


, * A mT  , , y y , . m+1  T rk  * 0,1,... 

for  lk  n~  T - 1 *k  n T’  {m  = 1,2 n-1  and  0<T<T/n. 


I 
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i! 


-1 


r 


r 

i. 


i 

i 

i 

i 

i 


Similarly,  let  PgB  be  the  covariance  of  eB  then 

Peta(t>  = E»i<J‘>  «n,T(‘)] 


Bnr  ' Bm 


(3-68) 


For  m=0,  substituting  (3-49)  and  (3-60)  into  (3-68)  gives 

PeB0(t)  = (Hi^F-H2)px(,<)(H1AF-H2)T 

+ [H1AG-EcCpiP2(T)]  Wk[H1AG-EcCpH»2(T)]T 
for  tk  + ^ <t<tk  + ^ + t,  k "0,1,...,  0<x<T/n. 

For  m 1,  substituting  (3-47)  and  (3-60)  into  (3-68)  gives 

PeBm(t)  = C(H1AF-H2)(AF)m]Px(k)[(H1AF-H2) (AF)m]T 

m-1 


(3-69) 


m-1 


+ [ (H-jAF-Hg)  Z (AF) 1 AG  + H1AG-E(;C  4»2(t)]W|cC(H1AF-H2)  Z (AF)1  G 


i=0 


i=0 


H1AG-EcCpiP2(T)]1 


(3-70) 


for>  t + T <t<t  + T+  0,1,... 

k n - k n m=l ,2, . . . ,n-l  and  0<x<T/n. 


From  Figure  6,  if  n=l  then  this  model  is  almost  the  same  as  the  basic 
model  except  that  wp  in  this  model  is  sampled  and  zero-order-hold  whereas 
wp  in  the  basic  model  is  continuous.  With  n=l,  there  are  only  one  value  of  m 
that  is  zero.  Therefore,  if  we  let  n=l  then  the  covariance  errors  of  this 
model  should  be  close  to  or  behave  the  same  as  the  covariance  errors  of  the 
basic  model.  In  the  next  section  (3.3)  we  will  apply  the  data  from  the  first 
example  in  section  2 (the  basic  model)  with  the  equations  of  this  model  and 
show  that  if  n=l  then  the  covariance  errors  of  this  model  behave  the  same  as 
the  covariance  errors  of  the  first  example  of  the  basic  model  as  expected. 

3.3  Example 

By  using  the  data  from  the  first  example  in  Section  2,  Ap=0,  B -j p=  1 , 
B2p=l,  Cp=1,  Fc=0,  V0’  Nc=0  and  Fc=-k.  For  the  zero-order  hold  of 
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N 


o2„ 


are  calculated  to  be 


F = 


1 


and 


G = 


(3-72) 


Substituting  (3-71)  and  (3-72)  into  (3-31)  and  (3-32),  the  matrix  F and  G 
are  calculated  to  be 


F = 


(1  - {p)"  0 0 


0 

0 


0 0 
0 0 


(3-73) 


and 


T 
n 

0 + 

0 

S' + n 


1(1  -Jsl, 

nv  n ' 
0 

0 


— ) + (1  - — )2  + 
0 
0 


1(1  . iSl)n 

nv  n ' 

0 

0 


+ (1  -£!)" 


I r1^  i 
0 
0 


(3-74) 
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where 


P = 1 - 


(3-75) 


The  steady-state  covariance  of  the  states  is  found  by  solving  (3-64) 


2-.,  n,2  r>  0 0 ' 

°w  Hl-P  ) 

pxss  = — — Tri ? 0 0 0 

XSS  n( 1-p  ; ( 1-p/ 

0 0 0 


(3-76) 


P . and  P 0 can  be  calculated  from  (3-66),  (3-65),  (3-67)  and  (3-70). 
eAmss  eBmss 

From  (3-35)  and  (3-36)  we  get 


H1  = [-k  0 0] 


(3-77) 


H2  = C-k(l-kr)  0 0] 


Therefore,  we  get 


kVa  2T  [1-Pn]2  k2x2a  2n 
w c J w 

PeAOss  n[l-P2n][l-P]2  T 


(3-78) 


(3-79) 


k4o  2T(1-Pn)2  T 2 k2o  2n  T 2 
p _ w /T  _ y /T  _ y.d 

eBOss  n(1.p2n)(1_p^Z  (n  T n 


(3-80) 


p . k4x2ow2TP2m  (1-Pn)2  k2ow2T2n  k2J2(1-lf)2  _ 2kT(  1-Pm)  -■ 
eAmss  n( 1-P2n)  (1-P)2  T n2(l-P)2  ^ 

(3-81) 

where  m=l ,2, . . . ,n-l . And 


eBmss 


4 T ,2  p2n’t’^T(,-Pn)2  k2°„2n.T  ,2  rl*k2T2(l-Pm) 


2kT( 1-P 

nTNP'T 


(3-82) 
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where  m=l ,2, . . . ,n-l . 


As  discussed  in  the  previous  section,  let  n=l.  Then  (3-79)  and  (3-80) 


becomes 


eauss 


,2  2 r kx  , T -j 

k to  [^-tv  + j ] 
W C-K  I 


(3-83) 


P . k2(T.T)  0 2 ^ + Hzii3 

eBOss  v ' w L2-kT  T 


(3-84) 


PeA0ss  and  PeAss  (e9uation  (2-28),  first  example  in  Section  2)  are  plotted 
in  the  same  graph  in  Figure  8 as  a function  of  x.  pegQss  and  PeBss  (e9uation 
(2-29),  first  example  in  Section  2)  are  plotted  in  the  same  graph  in  Figure  9 
as  a function  of  x.  From  the  graphs  of  Figure  8 and  Figure  9,  pe^Qss  and  PeB0ss 
behave  the  same  as  PeAs$  and  PeB$s  as  expected.  PeAOss  and  PeBQss  are  also 
equal  to  PeAss  and  PeBss  at  x=0  and  x=T,  but  they  are  less  than  PeA$$  and  PeBss 
for  0<x<T. 

In  this  model,  we  assume  to  use  the  same  value  of  the  covariance  matrix  of 
wp(tk)  for  every  interval  of  time  {(tk+l,tk),  (tk+l,tk+x),  (tk+x,tk)}.  It  isn't 
true  because  for  n=l  E[wp(tk)  wp^(tk)  = y (w  is  the  input  disturbance  covariance 
matrix)  is  derived  from  the  noise  in  period  T(tk+l,tk).  That  means,  we  can  use 
this  value  y in 

^2  ^2 
E[  / 4>(t 2,o)  Blpwp(tk)do  / wpT(tk)Bjp4>T(t2,s)ds] 
t1  t1 

if  the  limit  of  this  integration  t-|  and  t2  are  tk  and  tk+l.  But  in  the 
equations  of  this  model  we  use  y for  every  interval  of  time  {(tk+l,tk), 
^tk+1,tk+T^’  (tk+T»tk^  whl’ch  makes  the  difference  between  PeAss  and  PeAOss« 


eBOss 


and  P p r r * Tlle  reason  for  these  are  as  follows; 
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FIGURE  8 P,AoMAND  PeA„ 


FIGURE  9 P 
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With  data  from  this  example. 


k k T 

Ef  *(tk+T,a)B1pwp(tk)da  f wpT(tk)Blp  *(tk+T’s)ds] 


k k 

= / 4>(tk+T,o)B1  da  E[w  (tk)w  T(tk)]  / <J>(tk+x  ,s  )ds 


(3-85) 


* WT(y) 

where  E[wp( tk )wpT( = y 


(3-86) 


wt(j)  isn't  the  true  value  of  (3-85)  and  from  reference  1 the  true  value 
of  (3-85)  is  wt.  Because  t<T  then  wx(y)  < wx  and  tfi-’+T  then  wx(i~)*  wT;  that 
means  at  x=0  and  x=T  of  (3-86)  is  the  true  value  of  (3-85).  Therefore  we  can 
say  that  the  covariance  errors  of  this  model  are  the  approximate  values  and 
they  are  less  than  the  true  value  which  they  should  be. 

For  n>l,  E[wp(tk)wp^(tk)]  isn't  equal  to  wT  because  in  this  case, 
^twp(tk^wp^(tk^  is  derived  from  ttie  n°ise  in  period  ^[tk+^,tk]  then 
E[Wp(tk)Wp^(tk)].  The  following  expressions  are  the  covariance  errors  of 
this  model  when  n=2  (it  is  difficult  to  show  the  covariance  errors  when  n>  2]. 


For  n=2  (3-74),  (3-80),  (3-81)  and  (3-82)  become 


eAOss 


k4  w*t2T(1-P2)2  + 2k2x2aw2 
2[l7r][l-P]2  T 


(3-84) 


3 .fcVTM-P*)2  T t)2, 

eBOss  2[1-P4][1-P]T  2 


0.2  2 

2k  o T r\ 

T-®-  tj  - 


(3-85) 
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eAlss 


k4o2P2x2T(l-P2)2  2k2o  2t2  2,2 

— w~~  ~~a o~ + n + V1’ kT^ 

2(1-P  )(1-Pr 


(3-86) 


eBlss 


.4  2n2 


k o CP‘  T(  1-P  ) T ? 2k  o , o .2,2 
= w (T  _ t)2  + w_  (T  . t)2[1  + kT  . kT]  (3_87) 

2( 1-Pj (1-P)2  2 T 2 4 


PwA0ss  and  PeB0ss  are  plotted  in  Figure  10  and  PeAlss  and  PeBlss  are 
plotted  in  Figure  11.  All  of  them  are  plotted  as  a function  of  t.  The  dia- 
gram to  the  right  of  each  plot  shows  the  times  corresponding  to  the  values  of 
the  controller  outputs  used  to  calculate  eA  and  eg. 

The  result  obtained  for  this  example  agree  with  intuition.  For  small  t, 
the  variance  of  eA  is  small  and  that  of  eg  is  large.  For  large  t,  the 
complementary  situation  holds. 


KT(4-KT) 


FIGURE  10  P, 


KT(4-KT) 


FIGURE  II 


4.0  THE  DELAY  MODEL 


The  Delay  Model  is  close  to  the  basic  model,  except  there  are  time  delays 
in  the  control lers,  A/D  conversion  at  the  input  of  the  controllers  and  D/A 
conversion  at  the  input  of  the  plant.  Figure  12  illustrates  the  relationship 
between  the  plant;  the  controllers  and  the  time  delays  of  this  model. 

There  are  three  variations  (A,B,  and  C)  for  this  model  and  these  variations 
have  different  amounts  of  time  delay.  For  the  first  variation,  there  are 
three  time  delays;  at  the  controller,  the  A/D  converter  and  the  D/A  converter. 
For  the  second  and  third  types,  there  is  only  one  time  delay  in  the  system 
and  it  is  the  time  delay  in  the  controller.  The  time  delays  in  the  A/D 
converter  and  the  D/A  converter  are  neglected.  The  value  used  for  the  time 
delay  in  the  system  is  the  difference  between  the  second  variation  and  the 
third  variation. 

4.1  Variation  A:  System  Configuration  and  Dynamic  Equations 

There  are  three  time  delays  in  this  variation.  At  the  controllers,  the 
A/D  convertors,  and  the  D/A  converter.  Define 
<$c:  controller  computation  time 
£,:  A/D  conversion  and  transfer  time 

a 

6^:  D/A  conversion  and  transfer  time. 

The  following  assumptions  hold: 

1.  The  sum  (A)  of  three  time  delays  must  be  less  than  the  value  of  the 

skew  time.  That  is  A = £ +6  +£  . and  0< A<-r. 

a c d 

2.  The  computation  times  for  the  output  of  the  controllers  must  be 
complete  in  the  period  T. 

3.  From  2,  the  value  of  the  skew  time  must  satisfy  A<t<T-A. 


A OC  : ANALOG  -TO-  DIGITAL  CONVERTER 
OAC  : DIGITAL -TO-  ANALOG  CONVERTER 
S/H  : SAMPLE- AND -HOLD 

T : PILOT  -INPUT  SAMPLE  PERIOD  AND 
T : SKEW 


tk-l-tk=T 


FIGURE  12  BLOCK  DIAGRAM  FOR  THE  DELAY  MODEL 
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Figure  13  illustrates  the  time  diagram  of  the  state  variables  and  the 
output  of  the  controllers  of  the  system.  Table  1 explains  the  kev  events 


involved  in 

the  computations. 

Table  1.  The  Steps  of  the  Computed  Time  of  the  System 

y sampling  initiated  for  controller  1 and  w 
p sampling  initiated.  p 

V«a 

V5a 

yp  reaches  controller  1. 

VV5c 

Wsc 

ycl  computation  is  completed 

*kVc,ld 

ycl  reaches  plant 

VT 

VT 

yp  sampling  initiated  for  controller  2 

VT*5c 

VT*fa 

y reaches  controller  2 

P 

VT*V'c 

V"+V5c 

yc?  computation  is  completed 

For  the  aircraft,  actuator,  and  sensor  dynamics  we  may  write 


x„  = Ax  + B,u  + B0  w 
P P P IP  P 2p  p 


y = C x 
yP  P P 


(4-1) 

(4-2) 


for  which  the  solution  is 


xp(t)  = ,*,(t,t0)xp(t0)  + / 4>(t,s)Blpup(t)ds  + / 4>(t,s)B?pwp(t0)ds 

to  to 


(4-3) 


As  in  the  basic  model,  the  control  input  for  the  plant  is  the  output 
of  the  controller  1 and  the  plant  output  is  the  input  for  each  controller  at 
the  different  times.  Thus 


*cl(tk)  ■ V'k1 


(4-4) 


. 


1 


I 
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Vl^a*  " VV 

= Cpxp<tk) 

uc2(W6)  "W** 

= CpXp(tk+T) 


(4-5) 


(4-6) 


The  discrete-time  equations  for  controller  #1  are: 

xcl(VuV6c)  ' Fc*cl(  W6c>  + Gcucl<Vsa'  <4‘7) 

ycl<  Wsc>  ’ Hcxd<  Wsc>  + Ecucl(tk+<a>  (4-S) 

for  k = 0,1,...,  and  for  controller  #2. 

xc2(V'+’t4a+5c>  = Fcxc2(Vx+VSc>  + Gcuc2(tk+x+5a>  (4'9> 

yc2(VTF5atscJ  = Hcxc2(tk+T*sa*sc)  * Ecuc2(V,+V  (4'10) 

for  k = 0,1 ,. . . , . 


Substituting  (4-5)  and  (4-8)  into  (4-4)  gives  the  plant  input  in  terms 
of  the  plant  and  controller  1.  State  variables  as 

up(tk>  * Hcxcl(VHVac>  + Eccpxp(tk-’)  <4- 

Define  new  variables  xhpl  and  xhcl , in  order  to  have  the  variables  that 
change  in  a piecewise  constant  manner  (otherwise  we  would  have  had  to  include 
x (t.-l)  and  x 1(t.-l+6,+6  )).  There  are  two  additional  equations  that  are 

p K C I K a C 

requi red 
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and 


xhpl(tk+1)  ’ VV 


= xcl(V6a+6c) 


(4-12) 


(4-13) 


Then  (4-11)  becomes 


up(tk}  = Hcxhcl(tk+6a+6c)  + EcCpxhpl(tk) 


(4-14) 


At  time  tk+A,  a new  value  of  u^  takes  effect.  So  let  t0=tk  and  t 
then  (4-1)  becomes 


VA- 


X (tk+A)  = *(a)x  (tk)  + *1(A)up(tk)  + *2(A)wp(tk)  (4-15) 


where 


and 


(a)  = eAA  Je'At1B1  do 


/ \ _ A A " Ac1  p i 

^p(  A ) Q G d< 


Substituting  (4-14)  with  (4-15)  gives. 


xp(tk*M  = ♦U)xp(t|t)  * ♦,(»)HcXhcl<tk+sa*sc)  * *l(s)tcVhpl(tk> 

+ 4'2(A)wp(tk)  (4-16) 

At  t = tk+l , xp( 1 ) depends  on  the  value  of  *p(tk+A),  up(tk+A)  and 
wp(tk+A).  So,  let  tQ  = tk+A  and  t = tk+l.  The  equation  (4-1)  becomes: 

Xp(tk+1)  = i(T-A)xp(  tk+A)  + li-j  (T-A)up(tk+A)  + .2(T-A)wp(tk)  (4-17) 


By  (4-4)  and  (4-8); 
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up<Va>  = Hc*cl(  W5c>  * cccp,p(tk) 


Substituting  (4-16)  and  (4-18)  into  (4-17),  gives 

*p<V'>  ■ WT)  * '^-‘IVpVV  * ■»<T-a)*1<fl)tcCpxhp1(tk) 

♦ ♦,(T-i)Hcxcl(V‘,t«c)  * ♦(T-i)*,(4)Hcx|lc,(VJa«c) 

* [»( T-i  )^2  ( A ) ♦ *2  (T-a  ) Vp  ( tk) 

Substituting  (4-5)  into  (4-7)  and  (4-8)  gives 


and 


xcl(V'*sa*5c)  ■ Fcxcl(V5aV  +GcCpxP(tk) 


*cl(tk*5aV  * Hcxd<  W«c>  * EcCp«p<V 

Similarly  for  y 0 and  x 0,  substitutino  (4-6)  into  (4-9)  and  (4 

Ct.  Cc. 


gives 


and 


xc2(VU'4Vsc>  = Fcxc2(V'+V5c>  * GcCpxp(tk*’! 


l,c2<VTt‘atJc)  * Hcxc2<Vx*Vsc>  * EcCpxP<V’> 


The  quantity  x (t^+t)  can  be  written  using  the  solution  to  equation 
P K 

Wt)  = ^(t)  + *l(T_a)[:cWtk)  + *(T'A),i'l(A)EcCpxhpl(tk) 

+ *i(^)Hcxcl(tk\+jc>  + ♦(^A)*i(4)HcxhCi(V6a+«c) 

+ [<t>( T- a) 4^2 ( A ) + g;2( T- A) 3w  ( tk) 


1 

(4-18) 

(4-19) 

(4-20) 

(4-21) 

-10) 

(4-22) 

(4-23) 

(4-19)  as 

(4-24) 


By  substituting  (4-24)  into  (4-??)  and  (4-23).  The  controller  2 
equations  are  obtained  as 

xc2<V'+Ttsa+sc>  ' GcCp[,(t)  * *tcCp  (,-4)*1(4)CcCpxhpl(tk) 

♦ OcCp*,(,-s)Uc*c,(tk+Je«c)  + GcCp»(T-4)»,(s)Hcx„cl(tkt«atJc) 

+ Fcxc2^tk+T+6a+5c^  + GcCp^^T“A^2^A^  + MT-A ) ^wp(  V 

(4-25) 

and 

yc2(tk+T+6a+6c)  = EcCP[*(t)  + MT“A,EcCp]  + EcCp<I>(T_A)^l(A)ErcCpxhpl(tk) 

+ EcCp^l(T-A)"cxcl(tk+6a+6c)  + EcCp^T-A)t1(4)Hcxhci(tk+4a+6c) 

+ Fcxc2(tk+T+6a+6c)  + EcCpi:<t>(T-A)^2(A)  + ^(T-A)3wp(tk)  (4-26) 

The  inherent  error  is  defined  in  two  parts  as  follows: 

e„(t)  » yc,( VV'c>  - *c2(tk+’+Vsc>  <4'27> 

for  t|<+T+l5a+lSc  .1  t < ^k+^+4a+l5c’  ^ = 0,1,...,  0<t<T-a  and 

e„(t)  - yc,(VH«a*sc)  - yc2(y-T+6a+«c)  (4-28) 

for  tk+l+6a+6c  <.  t < tk+1  + T+6a+6c,  k = 0,1,...,  A<tj<T-A. 

Figure  14  shows  the  skewed  sampling  and  inherent  errors  of  the  Delay  Model. 
Channel  1 produces  the  sampled  outputs  at  times  t,+6  +6„,  t. +1  +6 , . . . , for 

K a C K a C 

k = 0,1,...,  and  channel  2 produces  the  sampled  outputs  at  times  tk+T+<$a+£c, 


V1+T+V6c'*" 
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th*T*®o*®c 


H CHANNEL  2 

lk  ‘2  »Ct 


i4tt) 

* »c 

. Ith»8«*8c  ) 

»c* 

Uh»T 

*V*c> 

FOR 

lk* 

T *9  *8  ^ 1 

Q C 

< t 

. ♦ i ♦ 8 

k a 

♦8c 

k = 0,l  ...  , 

A< 

r < T - A 

i,lt  ) *• 

»c. 

(tk*u8a  *V  - 

»ci 

(<„♦  r ♦ 

W 

FOR 

'■I.  * 

1 ♦»„ *8C  < t < 

*k» 

l * t ♦ SQ 

♦8c 

k = 0,  i . - . 

, 

r < T - A 

WHERE 

A 

8a  * 8c  * 8<f 

FIGURE  14  SKEWED  SAMPLING  AND  INHERENT  ERRORS 


r 


i 

\ 


11-47 


These  equations  can  be  put  in  compact  form  by  writing  them  in  terms  of 


a combined  stated  vector. 


x(tt) 


j xhpl (tk ^ 

wvw 

vi'vw 

Vj'V’W 


The  state  equations  become 


and  x(tkfl) 


«pV> 

xhptOk*" 

WV'W 

vdVW 


»(tk*l)  ■ F(T.t)  x(tk)  * G(T.t)  w (tk) 


where  F(T.-t)  is 


fn  - *(T)  + ^(T-A)  EcCp 
*12  " *(T-A)  ^(A)  ErCp 

fl3  3 Hc 

f]4  - 1>(T-A)4>](A)Hc 


(4-29) 


f?2  3 *23  * «24  3 *25  * 0 


31  c p 


f34  3 *35  3 0 
f41  " *42  * 0 
f43  3 1 


44  ‘ 1 45 


f5l  ■BcCPC*(T)  * V^Vpl 
fS!  " 6ccp»(T-»)*,(A)  CcCp 
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f53  " GcCp 

f56  * GcCp*  (T‘  >Va)Hc 


fC7  = F 

57  c 


and  G(T,t)  is 

9!  * (T-  )<MA)  *i|^(T-A) 

g2  = g3  = g4  = ° 

g5  = GcCp[  (t-a)^(a)  + «*(T-A)] 

The  controller  output  equations  are 


and 

>cl<VW  ■ Hl*(‘k> 

(4-30) 

*c2(VT+W  = H?x(tk)  + 0wp(V 

where 

(4-31) 

and 

H,  « [EC  0 H,  0 0] 

l c p c 

(4-3?) 

H2  » [EcCp(*(T)  ♦ *1(T-A)EcCp)  EcCp<K(T.A)l|i1(A)EcCp  E^  (t-AJH. 

Ecy(T-A)4»1(A)Hc  Hc] 

(4-33) 

P * EcCp[4>(T-A)(l.2(Al  + ^(t-A)] 

(4-34) 

The 

expressions  of  eA(t)  ep(t)  become 

eA(t)  - (HrH?)  x(tk)  - P wp(tk) 

(4-35) 

for  ti+i+S  +A  < t < ti+1+6  +A  , k * 0,1,...,  A <t<T-A  and 

K d C — K d C 


LI 
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eB(t)  = [H^d.i)  - H2]x(tk)  + [H1G(T.T)-n3wp(t|c)  (4-36) 

for  t|(+l+^a+,5c  < t < tk+l+T+$a+6c,  k = 0,1,...,  A< t<T-A 
Let  be  the  average  error,  so 

EA  = 2^eA^  " 

= ^H^H^d.T)  - 2H2]  x(tk)  + ]{H1G(T,t)  - 2n]wp(tk)  (4-37) 

4.1.2  Covariance  Analysis 

As  in  the  multirate  model,  the  input  wp(tk)  is  a Gaussian  white  noise 
random  process  with  zero  mean  which  is  independent  of  x(0).  Then  (3-54), 

(3-55)  and  (3-56)  are  repeated  here  as 

E[wp(tk)  = 0 (4-38) 

E[x(tk)wpT(tk)]  = 0 (4-39) 

E[wp(tk)wpT(tk)]  = Wk  (4‘40) 

The  covariance  matrix  of  the  states  is  defined  as 

Px(k)  = E[x(tk)xT(tk)]  (4-41) 

and 

Px(k+1)  = E[x(tk+l)xT(tk+l)]  (4-42) 

= F(T,T)Px(k)FT(T,T)  + G(T,T)WkGT(T,r)  (4-43) 

The  steady-state  covariance,  designated  Px$s  is  found  by  solving  the  equation 

Pxss  = F(T’t)  Pxss  fT(T*t)  + g(T’t>  wkGV.T) 


(4-44) 


JT  " W*.  ‘-"I"  ■■■tl  liU  M.U^  M ^ _ I.L  l_. 
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[ I 

h 

r Covariance  of  the  Errors 

The  covariances  of  eA( t ) and  eg ( t ) are  calculated  using  the  same  pro- 
cedure as  in  the  previous  development.  Let  PpA  be  the  covariance  of  e^, 

' then 


PeA^)  - E[eA(t)eAT(t)] 

(4-45) 

for  tk+x+Sa+6e  < t < tk+1+fa+fc.  k = 0,1,...,  A<t<T-a 

Substituting  (4-35)  with  (4-45)  gives 

* 

Pefl(t)  ■ HAPx(k)HflT  ♦ "»/ 

for  tk+T+lSa+<sc  < t < tk+l+lSa+<f;c>  k = 0,1,...,  a<t<T-A 

where 

(4-46) 

Ha  = HrH2 

Let  Pfig  be  the  covariance  of  e^,  then 

(4-47) 

PeB{t)  = E[eB(t)eBT(t)] 

for  t|c+l+,5a+fc  < t < tk+1+T+,5a+,5c»  k = 0,1,...,  A<^<T-A 

Substituting  (4-36)  into  (4-48)  gives 

(4-48) 

PeB(t)  • HBPx(k)HBT  * eBwkeBT 

for  tk+1+fa+5c  < t < tk+l+T+6a+<$c,  k = 0,1,...,  A<t<T-A 

where 

(4-49) 

Hb  = H^d.x)  - H2 

and 

(4-50) 

"B  * H-|G(T,t ) - o 

(4-51) 
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From  Figure  12,  if  the  value  of  the  time  delays  (6^6,.,$^)  are  equal  to 
zero.  The  model  is  close  to  the  basic  model  except  that  the  input  wp  in 
this  model  is  sampled  and  zero-order  hold.  Then  if  the  time  delays  of  this 
variation  are  equal  to  zero,  the  equations  of  this  variation  should  be  close 
to  the  equations  of  the  basic  model  or  the  equations  of  the  multirate  at  n=l. 

Let  A=0,  then  the  term  x^  and  x^  become  xp  and  x^  or  x(t^)  and 
x ( t^-*- 1 ) of  this  variation  become 

VV  + xhpl(V  VV1*  + xhpl^k+1^ 

x(tk}  = xcl(tk)  + xhcl(tk}  • 4nd  x(V')  a xcl(tk+1)  + xhcl(tk+1) 

Xc2(tk+T)  xc2(V1+t) 


Then  F(T,i)  and  G(T,t)  becomes 


♦ (T)  + i|»1(T)EcCp 


*i(T)Hf 


HT,t)  = 


GcCU(t)  tflh)EC  ] GC^(t)Hc  0 


^(T) 

G(T.-r)  = 0 


GcCp  2(t) 


where  ^,(«c)  = t2(^c)  = 0. 


The  matrix  F(T,t)  is  identically  to  F(T,t)  in  the  Basic  Model  and  the 
matrix  G(T,t)  is  close  to  the  second  matrix  of  x ( t ^ + 1 ) in  the  Basic  Model.  This 
is  because  the  noise  of  this  model  is  sampled  and  zero-order  hold  but  the  noise 
of  the  Basic  Model  is  continuous.  Now,  examine  H1  and  H^,  with  A=0.  The 
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equations  of  Hi  and  form  the  equations  (4-32)  and  (4-33)  become 

H,  = [EcCp  Hc  0] 
and 

H2  - tEcCp[*(x)  * *,(T)EcCp]  EpCp*,(T)Hc  Fc] 

These  two  matrices  are  the  same  as  Hi  and  H^  in  the  basic  model.  By  this 
technique,  the  equations  of  this  variation  (time  delay  = 0)  are  close  to 
the  equations  of  the  basic  model  and  they  are  the  same  as  the  equations  of 
the  multirate  model  as  expected.  By  using  the  data  from  the  first  example 
of  section  2 (Basic  Model),  we  can  study  the  characteristic  of  the 
covariance  errors  of  this  variation.  As  discussed  above,  with  data  from  the 
first  example  of  section  2,  if  we  let  the  time  delays  in  the  expressions  of 
the  covariance  errors,  the  results  must  be  equal  to  the  covariance  errors  of 
the  multirate  model  (n=l). 


4.1.3  Example 


The  data  from  the  first  example  of  section  II  are 


Ap  = 0 
Blp  = 1 


C = 1 
P 

:et  wk  = ow2/T. 
♦(t,t0)  = 1 


Fc  = 0 
Gc  =0 
Hc  = 0 
Ec  = -k 


ihCO  ■ t2(T)  = T 

^ 1 (T-a)  = ||»2(T-A)  = T-A 

l|>1  (T)  = \p2(T)  = T 

(t-a)  = 4>2(t-A)  = T-A 


* 


r 

i. 


i. 

[. 


[ 

t 

i 
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By  using  the  equations  of  matrices  F(T,t)  and  G(T,x)  with  these  data, 
we  have 


F = 


l-k(T-A) 

1 

0 

0 

0 


-ki  0 0 0 

0 0 0 0 


0 0 
0 1 


0 

0 


0 0 0 


0 

0 

0 


T 
0 

G = | 0 
0 
0 


The  steady  state  covariance  of  the  states  if  found  by  solving  (4-44) 


xss 


o 2T 
w 


l-(l-kT)* 

0 

0 

0 

0 


0 0 

0 0 
0 0 
0 0 
0 0 


ow2T 

T-(l-kT)O-kT)  0 

0 0 

0 0 

0 0 

0 0 


By  using  equations  (4-32),  (4-33)  and  (4-34)  we  have 


H1  = C-k  0 0 0 0] 

H2  = [-k+k2(T-A)  k2A  000] 


n = -kx 


Therefore  we  get 


(4-52) 


(4-53) 


(4-54) 


(4-55) 

(4-56) 

(4-57) 


■XI 


k3(T-4)202  kVo? 


eAss  (2-kT) 


and 


W * ‘'(T-xy,.*  [|SfH  * ^ 
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(4-58) 


(4-59) 


1 

f 


Let  A=0.  Then  we  have  the  expression  of  PeAss  and  PeBss  of  the  Basic 
model  with  discrete  noise. 


and 


Peflssl  ■ 'W  tFFT  * T] 

PeBssl  * SV-'K2  ^ 


(4-60) 


(4-61) 


These  covariance  errors  are  equal  to  the  covariance  errors  of  the  multi- 
rate model  when  n=l  (equations  (3-83)  and  (3-84))  as  expected.  As  the  same 
reason  in  the  multirate  model,  the  covariance  errors  of  this  variation  are 
the  approximate  values  and  they  are  less  than  the  true  value. 

PgAss  (4~58)  and  PeBss  (4-59)  are  plotted  in  Figure  15  as  a function  of  t, 
The  diagrams  to  the  right  of  each  plot  show  the  times  corresponding  to  the 
value  of  the  controller  outputs  used  to  calculate  e and  eg.  From  (4-58)  and 
(4-59),  PeAss  depends  on  the  value  of  the  time  delays  but  Pegss  doesn't. 

Under  the  specific  case  of  the  zero  value  of  Ap  makes  Ppgss  does  not  depend 
on  the  time  delay,  but  for  other  values  of  Ap  both  PpAss  and  Ppgss  depend 
on  the  time  delays.  However,  in  this  example,  the  results  obtained  for  this 
example  agree  with  intuition.  For  small  t,  the  variance  of  e^  is  small 
and  that  of  eg  is  large.  For  large  r,  the  compl imentary  situation  holds. 


k * «r  * ( T - 2 A >*♦  h*(  T - A )acr 


tfc  tk**Q*8 


t^*l  ^)j*2 


FIGURE  15 
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4.2  Variation  B 

Physically,  the  value  of  time  delay  at  the  controller  (fQ)  is  a high  value 
and  it  is  greater  than  the  values  of  time  delay  at  the  Vt>  converter  (6  1 and 

a 

at  the  d/a  converter  (5^).  Then  in  this  variation  and  the  next  variation  of 
this  model,  and  6^  are  neglected  and  the  value  of  $ is  assumed  to  be 
close  to  the  value  T.  For  this  variation,  is  still  less  than  t (value  time). 
Then  t is  also  a high  value  and  the  limit  of  t is  <$c<t<T.  Because  of  the  high 
value  of  5c  then  the  computation  of  y^  is  not  complete  in  the  same  period  of 
the  computation  of  yc^  but  it  is  complete  in  the  consecutive  period  (t+<$c>T). 
Figure  16  shows  the  computations  of  yc^  and  ycp. 

From  Figure  16,  if  6C  is  changed  to  be  A,  this  variation  is  close  to 
variation  A,  except  that  the  computation  of  yc?  is  not  complete  in  the 
same  period  of  the  computation  of  ycj.  However,  the  input  of  the  plant  of 
this  model  is  the  output  of  the  first  controller  which  is  the  same  as 
variation  A.  Therefore  the  equations  of  variation  A can  be  used  in  this  varia- 
tion by  using  6 instead  of  A and  using  5 and  in  the  equations  equal  to 
zero.  The  equations  shown  below  are  just  for  calculating  the  covariance  errors 
of  this  variation. 


'xp(tk> 

"yv 

x„p(tk> 

xi.Pi(V'> 

xd<W 

,x(tk+l)  = 

xcl(V,+*c> 

xhcI<W 

Vi'VV 

,xc2(tk*T*5c). 

_xc?<VlMV 

and 
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x(tk+l)  = F(T,T)x(tk)  ♦ G(T,t)  wp(tk) 


where  F(T,t)  is 

f11  a»(T)  + *,(T  «c)EcCp 
f12  -♦(t-«c)*1(«c)EcCp 

f13'*l(T-6c>Hc 
f14  -•(T-«C)*1(«C)HC 

f15  = o 
f21  - 1 

f22  = f23  = f24  = f25  = 0 

f31  = GcCp 

f32  = ° 

f33  3 Fc 

f34  3 f35  = 0 

f41  = f42  = 0 

f43  = 1 

f44  * f45  ’ 0 

fsi  ■ W#(T)  + 

f52  3 W^VVWp 
f53  = Wl(T‘6c)Hc 
f54  3 GcCp*(T"6c),pl(T)Hc 
f55  3 Fc 
and  G(T,t)  is 

g1  = *(T-6c)4>2(6c)  + 

g2  = g3  3 94  = ° 

95  3 GcCp  t>(T-6c)  W + *2(t"8C)] 


The  controller  output  equations  are 


ycl^W  ” W 

and 

1'c2(tktI+{c)  ""  H2x(tk>  + 0 "p^k1 

where 

H,  > [EcCp  0 Hc  0 0] 

and 

«2  - [EcCp(»(T)  * *,<T-*c)EcCp)  KcCp*(t-«c)*1(«c)  EcCp 
EcCp(*,(T-«c)Hc  EcCp»(t-«c)*,(Sc)Hc  Hc] 

o = EcCpC<I>(T-6c)vp2(  + vMT"fiC^ 
eA(t)  = (HrH2)  x (tk)  - n wp(tk) 

for  tk+r+6c  < t ct^+l+fi^,  k = 0,1,  A<t<T  and 

eB(t)  = (H^T.t)  - H2)  x(tR)  + [H-|G(T,t ) - n]wp(tk) 

for  t. +1+6  < t < t. +1+6  , k =0,1,  A<t<T.  The  average  error  of  e.  and  e 
k c — k c n 

EA  = 2 ^A^  + 

The  covariance  of  states 

Px(k+1)  = F(T,t)  Px(k)  FT(T,t ) + G(T,r)wkGT(T,T) 

The  covariance  of  errors 

PeA<t>  ■ HAPx(k’HAT  * pVT 

for  t.+l+6  < t < ti+1+6  , k =0,1,..., 
k c — * c 


A<t<T. 
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where 

HA  " H,  - H? 
and 

PeA(t)  " MBPx(k)  HRF  + nBwk°BT 
for  t^  + l+6c  < t < t^  + 1 + t + (Sc  , k * 0,1,...,  A<i<T 

where  HR  * H^F  - Hp 
and 


Recall  the  expressions  of  P 
we  have 

P , a 
eAss 

for  t.+i+A  <t<t.+l+A  , k *0,1, 

K C K C 


eAss  and  PeBss 


3 (^c>V 
' 2-kT 


A <t<T. 
c— 


from  the  example  in  variation  A, 


,2t2„  2 
k T o 

w 


And 


for  t.+l+A  <t  t.+l+T+A  , k * 0,1,..  ,,A  <t<T. 

^ c~  k c c 

PeAss  and  PeBss  are  P^ot*ed  ’n  Figure  17  as  a function  of  t.  The  dia- 
grams to  the  right  of  each  plot  show  the  times  corresponding  to  the  value  of 
the  controller  outputs  used  to  calculate  e^  and  e^.  Variation  A,  PeBss 
depends  on  the  time  delay  Ac  but  P(,pss  floes  not.  The  results  obtained  for 
this  example  agree  with  intuition.  For  small  t,  the  variance  of  e^  is  small 
and  that  of  e^  is  large.  For  large  t,  the  complimentary  situation  holds. 
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4.3  Variation  C 


The  same  as  variation  B,  we  have  only  one  time  delay  $c  but  in  this 
variation,  the  value  of  $c  is  greater  than  the  value  of  other  time  (t). 

There  are  two  cases  of  this  variation: 

1.  The  computation  of  yc?  is  completed  in  the  same  period  of  the 
computation  of  yc]  (f+^T) . 

2.  The  computation  of  y ^ isn't  completed  in  the  same  period  of  the 
computation  of  yc;>(T+<Sc>^  • 

4.3.1  Variation  C (Case  I) 

Figure  18  shows  the  computations  of  ycl  and  The  limit  of  the  skew 

time  is  0<t<6c.  In  this  case  the  output  of  the  first  controller  is  the  in- 
put of  the  plant.  Therefore  the  equations  of  this  case  are  almost  the  same 
as  the  equations  in  the  variation  B except  the  equation  of  XpU^+T).  Because 
the  output  of  the  first  controller  is  the  input  of  the  plant  then  the  equa- 
tion of  xp(tk+l)  in  this  case  is  the  same  as  the  equation  of  xp(t)(+l)  in 
variation  B and  from  variation  B. 


y V>  ■ C*<T>  * 4 *(T-«c)*,(«c)EcCpxhp,(tk) 

* *,(T-«c)Hcxcl(tk«c)  * ♦(T-«cH,(«c)Hcx|lc)(V«c) 

♦ C*(T-«c)*2(«c)  * *j(T-«cn  wp(tk) 

The  controller  equations  are 

xd<v'45c>  ■ Fc*ci<vy 4 GccPv(*k> 

*ci<vy  ■ Vci'w 4 yyyy 

xc?(tk4,«4fc)  ■ Fcxc2(tk*T*5c)  * GcCpxp(tk*T) 


and 


COMPUTATION  OF 


COMPUTATI 


FIGURE  18  THE  COMPUTATIONS  OF 


t 
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1 


h 


i 


>c2(tk+T+6c)  = Hcxc2(tk+T+6c)  + EcCpxP(tk+T) 

In  this  case,  the  equation  xp(tk+x)  is  derived  by  letting  t Q=t ^ and  t=tk+x. 
Then  from  the  plant  equation: 

xp(tk+T)  = *(x)xp(tR)  + ^(t)  up(tk)  + t2(x)  wp(tk) 

Substituting  up(tk)  into  the  equation  of  xp(tk+T)  gives 

Xp(tk+T)  3 *('r)xp(tk)  + V^VhclW  + 'P1(T)EcCpxhpl(tk)  + ^2(T)wp(tk) 


By  substituting  the  value  of  xp(tk+x)  into  the  equation  of  xc2(tk+l+6c+x)  and 

yc2(tk+T+,5c^  we  hdve 


xc2<V1+6c)  = Hcxc2<V^c>  + W(T)xp(tk)  + GcVl(T)Hcxhcl(V5c) 
+ GcVl(T)EcCpxhpl(tk)  + 6cCpVT)wp(tk> 

= Hcxc2(VT+V  + + Wl,T)Hcxhcl(VfC) 

+ EcVl(T)EcCpxhpl(tk)  + EcVMT 


’W 

‘ xp'V> 

xhpT(tk> 

xhpi(V'> 

WV5c* 

and  x(tk+l)  = 

xcl<tkt,t5c> 

xhcl<V*c> 

xhcl*tk*',Sc) 

xc2(VttSc> 

xc2(VUt+Sc' 

The  states  equation  is 

x(tk+l)  = F(T,x)  x(tk)  + G(T,x)  wp(tk) 


Then  F(T,t)  is 


f11  =*(T)  + ^l(T-6c)EcCp 

f12  *«T‘*cWWp 

f 13  ■♦l<T‘6c>Hc 
f,4  - •(T-«C)»1(«C)HC 

f15  * 0 
f21  = 1 

f22  = f23  = f24  = f25  = 0 

f31  * GcCp 

f32  = 0 

f33  = Fc 

f34  = f35  = 0 

f41  = f42  = 0 

f43  = 1 

f44  = f45  = 0 

f51  = GcV(T) 

f 52  = GcCp^l^x)EcCp 
f53  = 0 

f54  = GcCp*l(T,Hc 
f55  = Fc 
and  6(T,t)  is 

* <t>( T-6C )\(/2 ( 5C ) + 

g2  = 93  = 94  = 0 

g5  = GcV2(T) 

The  controller  output  equations  are 


»cl<V5c>  ■ 


and 

yc2(VT+fc}  = + °W 

where 

H1  ■ [EcCp  0 «c  0 0] 

and 

H2  ■ [£cCp«(T)  EcCp*,(T)EcCp  0 EcCp*,(T)Hc  Hc] 

The  inherent  errors  in  this  case  is  defined  the  same  as  variation  B. 

Then  from  variation  B 

eA(t)  = (H1-H2)x( tk)  - owp(tk) 

for  tk+T+«c  < t < tk+l+«c,  k =0,1 0^  <fc 

and 

eA(t)  = (H1F-H2)x(t|<)  - (H1G-o)wp(tk) 

for  tk+l+6c  < t < tk+l+T+6c,  k = 0,1,...,  0<t<Sc. 

Figure  19  shows  the  skew  sampling  and  inherent  errors  of  this  variation. 
The  average  error  is 

EA  = 7 ^A(t)  + eB^^ 

The  same  as  variation  B,  the  covariance  of  states  is 
Px(k+1)  = F(T,x)Px(k)FT(T,T)  + G(T,t)wpGT(T,t) 

and  the  steady  state  covariance  is 

Pxss  = F(T’t)PxssfT(T*t)  + G(T,T)wkGT(T,x) 


| 


k 


r 


! 

I 


f ' 

! 


I. 
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The  equations  of  the  covariance  errors  are 


pe»<‘>  * >Wk>HAT  * -V 


for  tk+t+6c  < t < tk+l+5c,  k = 0,1,...,  (kK{c 


where  = Hj-H^  and 


Pe8(t)  = HBPx^k^HBT  + WbT 


for  tk+l+fc  < t < tk+l+T+6c,  k = 0,1,...,  CKt<6c 


where  Hg  = H^F  - H? 


ng  = H-jG  - o 


The  same  as  the  variations  A and  B,  if  6 is  equal  to  zero,  the  equations 


of  this  case  must  be  equal  to  the  equations  of  variation  A or  B which  the  time 


delays  of  these  variations  are  equal  to  zero.  If  Sc  is  equal  to  zero  then  xhp1 


and  xhcl  become  to  xp  and  xcl  and  x(tk)  and  x(tk+l)  becomes 


VV  + '"'hpl^V  VV1*  + xhpl(tk+1) 

xUk)  = *cl(tk)  + xhcl(tk)  and  x(Vl)  = xcl(tk+l)  ♦ xhcl(tk+l) 

Xc2(tk+T)  xc2<V1+T> 


and  the  matrices  F(T,t)  and  G(T,t)  are 


*(T)  + (T) EcC 


»h  (T)H, 


ScCp[»(T)  * *,(T)EcCp]  GcCp*,(t)Hc  Fc 


1 
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<MT) 


G = 0 


GcV2(t) 


where  4»1(6(;)  = ^(^c)  = 0 and 


H,  - [EcCp  Hc  0] 


K2  ■ CEcC  [♦{t>  * *,<t)EcC  ] EcCp*,(T)Hc  hc] 


I 


The  equations  of  this  case  when  6 is  equal  to  zero  are  the  same  as  the 
equations  of  the  variation  A or  B when  the  time  delay  of  these  variations 
are  equal  to  zero  as  expected. 


4.3.2  Example 


From  the  data  of  the  first  example  of  Section  II  ,, 

V 0 Fc  • 0 " ■ °J  then  wk  ■ -f- 


Blp  ° ' 
B?P  * 1 

V 1 


Gc  -0 
Hc  - 0 
Gc  = -k 


l-k(T-«c)  -k«c  0 


0 0 


0 0 
0 0 


0 0 
0 0 
1 0 
0 0 


xss 


a 2T 

w 


l-(l-kT)‘ 

0 

0 

0 

0 


Therefore, 


for  0<T<fc  . 
and 


eAss  ’ T 


o 2t 

w 


l-( 1-kT) ( 1-kT) 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 


H1  =«  [-k  0 0 0 0] 

H_  * [-k  k2t  000] 


k4(T-«J2ow2T  k2(T-T)2ow2 

~7 


p * c;  w + 

eBss  ' kT(2-kT) 


for  0<t<{ 
“ c 
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i 


( 


I 

r 
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PeAss  and  PeBss  are  P1otte<1  in  figure  ?0  as  a function  of  T.  The  dia- 
grams to  the  right  of  each  plot  show  the  times  corresponding  to  the  values 
of  the  controller  outputs  used  to  calculate  eft  and  eg.  The  same  as  varia- 
tion A f*egss  depends  on  but  Pegss  does  not.  The  results  obtained  for 
this  example  agrees  with  intuition.  For  small  T,  the  variance  of  e^  is 
small  and  that  of  eg  is  large.  For  large  r the  complementary  situation 
holds. 

4.3.3  Variation  C (Case  II) 

Because  the  input  of  the  plant  is  the  output  of  the  first  controller, 
therefore  there  isn't  any  effect  in  the  system  when  the  computation  of 
isn't  complete  in  the  same  period  of  the  computation  of  yc^.  Therefore, 
the  equations  of  this  case  are  the  same  as  the  equations  of  Case  I.  The 
difference  is  the  value  of  the  sum  of  fc  and  t,  in  this  case  $c+t<T. 


r 


I: 

I! 

I! 
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5.0  THE  OUTPUT-AVERAGING  MODEL 

This  model  is  the  same  as  the  Delay  Model  except  that  the  input  of  the 
plant  is  the  average  outputs  of  the  controllers.  There  are  three  variations 
in  this  model,  variation  D,  E and  F.  Variation  D is  nearly  identical  to  the 
basic  model  (the  value  of  the  t'ime  delay  is  equal  to  zero).  Variations  E and 
F are  nearly  identical  to  variations  B and  C,  respectively.  Figure  21 
illustrates  the  relationship  between  the  plant  and  the  controllers  of  this 
model . 

5.1  Variation  D:  System  Configuration  and  Dynamic  Equations 

Figure  22  shows  the  time  responses  of  y^  and  y^.  They  are  the  same  as 
the  basic  model,  the  aircraft,  actuator  and  sensor  dynamics.  Therefore 


and 

Xp(t)  . ♦<t,t0)xp(t0)  + *,(t,t0)„p(t0)  * *?(t.t0)«p(t0) 

(5-1) 

V1'  ■ Wl) 

(5-2) 

In  this  model  the  input  of  the  plant  is  the  average  output  of 

the  con- 

trol lers. 

V'k)  t*c2<tk>J 

(5-3) 

uci<V  ■ V*k>  ■ cPxP(tk> 

(5-4) 

uc2<V>  = WT)  = CpWt) 

(5-5) 

The  discrete-time  equations  for  controller  number  1 are: 

xcl(V>  ’ Fc*cl<‘k>  * Gcucl(tk> 

(5-6) 

»Cl<*k>  ■ Hc*Cl'lk>  * EcUcl<V 

(5-7) 

for 

k = 0,1,...,  and  for  controller  number  2, 

FIGURE  21  BLOCK  DIAGRAM  FOR  THE  OUTPUT  - AVER  AGE  MODEL 


r— 


11-76 


xc2<V,tT>  ■ Fcxc2<Vt)  + 

yc2(tk*r)  ■ Hcxc2(tk+T)  * EcUc2(tk«) 
for  k *0,1,...  . 

From  Figure  22  = yc2(tk'1+T) • Then  by  using 

*c2(tk>  * *c2(VltT)  ’ Hcxc2(tk-,tT>  + EcUc2(V1tT> 
Substituting  (5-5)  into  (5-10)  gives 

yc2(V  = yc2(  V1+t)  = Hcxc2(tk'1+T)  + EcCpXP(V1+T) 


(5-8) 

(5-9) 


(5-10) 


(5-11) 


Define  new  variables  x^  and  x^  j in  order  to  have  variables  that  change 
in  a piecewise  constant  manner.  There  are  additional  equations  that  are 
required,  namely, 

Xhc2(tk+1+T)  = xc2(Vt)  (5~12) 

and 

Xhp2(tk+1+T)  = Xp(tk+T)  (5_13) 

Then  (5-10)  becomes 


yc2^V  = yc2^tk"1+T^  " Hcxhc2^k+T^  + EcCpxhp2^k+T ■ 
Substituting  (5-5)  into  (5-7)  gives 


ycl(V  = Hcxc1*V  + EcCpxP(V 


(5-14) 


f 5-15) 


Substituting  (5-10)  and  (5-15)  into  (5-3)  gives 


W = Hcxcl^V  + EcCpVV  + Hcxhc2^k‘T)  + EcCpxhp2^k+x)  ^5"16^ 


f 

. 
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At  time  t^+r,  the  value  of  upchanges  and  lets  t=tk+T  and  tQ=tk, 

(5-1)  becomes 

WT)  =$(T)xp(t|<)  +'*'i(T)up(tk)  +'^2^T^wp^tk^  (5-1?) 


Substitution  of  equation  (5-16)  into  (5-17)  gives 


Vt) 


*,(t) 


xp(tk+T)  = [*(t)  + ~2  EcCp]xp(tk)  + -y—  Hcxc](tk)  + 
*i(t) 

*~i Hcxhc2(tkM>  + *2<T>Vtk) 


*,(t) 


EcCpxl,p2(tk4T> 


(5-18) 


Now,  let  t=tk+l  and  t0=tk+r.  Then  (5-1)  becomes 

xp(tk+l)  ■ *(T-t)x  (tk+r)  +'I'1(T-t)u  (tk+r)  + ^2(T-T>wpCtk) 


From  Figure  22: 


u p(VT)  "I^c^V  +yc2(tk+T)] 


Substituting  (5-9)  an u (5-15)  with  (5-20),  gives 


Substitution  of  (5-11)  into  (5-19),  gives 

♦,(T-t)E 

xp(tk+l)  = [•(T-x)  ♦ -V-  EcCp]  WT>  + 


♦ ^T-t) 


Hcxcl (tk5 


i ( T-t  ) 


VT-t) 


+ 2 EcCpW  + -V-  Hcxc2<VT>  + VT-T>wp(tk> 


Substitution  of  equations  (5-18)  with  (5-22)  gives 


(5-19) 


(5-20) 


VVO  ■ 7 CHcxc,(tk)  ♦ Ecc  x (tk)  * Hcxc2(tk*T)  + Ecc  X (tk+T)J  (5-21) 


(5-22) 
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W'>  ' 4 "’I1!11  Eccp][t(,>  * *1  ¥ Eccp]  * *1  ^T1  Eccp>xp(tk> 

4 {[♦(T-t)  * ♦,  ^ EcCp]  *,  M Hc  * *,  Hc)  xc]{tk) 

4 (*(T-r)  + *,  Ufl  *,  ^ EcCp  xwp2(tk+T) 

+ 4 ♦,  ^ W +1  ¥ Hc  *hc2<V'>  4 *1  ^ Hc  *c2<Vt> 

(5-23) 


+ {[4>(T-t ) + ^ EcCp]  4>2(t)  + *2(T-t)}  wp(tk) 


Substituting  (5-4)  into  (5-6)  and  substituting  (5-5)  into  (5-8)  gives 


xcl(tk+1)  = F<-  ViUiJ  + GrCnxn(tJ 


c Aclltk'  T ucVp'  k' 


and 


Xc2(tk+1+T)  = Fcxc2(Vt)  + GcCpxp(tk+1) 


(5-24) 

(5-25) 


Substituting  (5-78)  into  (5-25)  gives 

*e2<V,tT>  - Fcxc2(tk)  ♦ Scc  C*(T)  * *,  ^ EcCp]x„(tk)  * 6cC  *,  £Vcl<tk> 


PJ  P' 


* GcCp  ^ (t)  EcCpx„p2(VT)  ♦ GpCp  ^ MHcxhc2(tkT) 


4 GcCp  2<T)"p(tk> 

The  inherent  error  is  defined  in  two  parts  as  follows: 

eA(t)  = yci(tk)  - yC2^h+j) 


(5-26) 


(5-27) 


for  tk+x<t<tk+l,  k = 0,1,...,  0<t<T,  and 


eB(t)  = ycl(tk+D  - yc2(tk+T) 


(5-28) 


for  tk+l<t<tk+l+x,  k = 0,1,...,  0<t<T. 


* 

r 


14 
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Figure  23  shows  the  skewed  sampling  and  inherent  errors  of  this  case. 

These  equations  can  be  put  in  compact  from  writing  them  in  terms  of  a 
combined  stated  vector. 


W 

x(tk>  ■ \>2(Vt) 

\:2(lk‘T) 

xl,c2(tk*T) 


, X(tk+1) 


Vv]) 

xcl(tk+1) 

xhp2(V1+T) 

xc2(V1+t) 

Xhc2(tk+1+T' 


The  state  equations  become 

x(tk+l)  - F(T,' ) x(tk)  ♦ G(T,t)  wp(tk) 
where  F(T,t)  is 

f„  ■ [»(T-T)  * (T)  * *,  -^Ccp]  * «.,^£CC( 

fl 2 * C»(T-t)  * *,^cCp]  Hc  ♦ Hc 

f13  ■ [♦(T-t)  * *l^fcCp 

fl4  = *1,T-t)  Hc 

f,5  ■ C«(T-t)  * Hc 

f21  * GcCp 
f22  = Fc 

f23  = f24  = f25  = 0 
f31  = + ^l^cCp 

f32  = *l‘^€cCp 


(5-29) 
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, CHANNEL  I 

V v2  yci 

CHANNEL  2 

I 1 1 1 1 y 

tk  tk*'2*T  tk*2  tk>2  + T C2 

= yci  UK)  - yc2  <VT) 

FOR  tk*  T < t < tk+ I k = 0,l,...  ,0<  T < T 

*B  (1)  = *CI  <V'>  ‘VC2  <tk  + T) 

FOR  tk*l  < 1<  iktUT  ksO.I, . . . , 0 < T < T 

FIGURE  23  SKEWED  SAMPLED  AND  INHERENT  ERRORS 
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K,  - [EcCp  Hc  0 0 0] 


H?[EcCp[«(T)  * *,^EcCp)  EcCp*,ili«c  EcCp*,iliEcCp  H( 


(5-32) 


(5-33) 


n = EcCp  2(t) 


(5-34) 


The  expressions  of  eA(t)  and  eg(t)  become 


eA(t)  = (H1-H2)x(tk)  - nw  (tk) 


(5-35) 


for  tk+r<t<tk+l,  k = 0,1,...,  0<t<T,  and 


eB(t)  = (H]F-H?)x(tk)  + (H,G-n)Wp(tk) 


(5-36) 


for  tk  + l<t<tk+l+T,  k = 0,1 0<t<T. 


5.1.1  Covariance  Analysis 

The  analysis  is  the  same  as  the  delay  model,  then 


Pk(k+D  = F(T,T)Px(k)FT(T,T)  + G(T,x)wkGT(T,T) 


(5-37) 


The  steady-state  covariance  is 


Pxss  = F(T-T)PxsspT(T’T)  + G(T,T)wkGT(T,T) 


where 


Px(k)  - E[x(tk)xT(tk)] 


wk  - E[«p(tk)wpT(tk)] 


(5-38) 


PeA<*>  ■ HAPx<k)HAT  * "VT 


(5-39) 
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for  tk+T<t<tk+l,  k = 0,1,...,  0<t<T 


where 


HA  " H1  • H2 


(5-40) 


PeB^  = HBPx^k^HBT  + nBwk°B 


(5-41) 


for  tk+l£t<tk+l+T , k = 0,1,...,  0<t<T.  Where 


Hg  = h]F-h2 


(5-42) 


Og  = H^G-r> 


(5-43) 


5.1.2  Example 

The  data  from  the  first  example  of  Section  II  are: 


Ap-0 

“ip” 

B2P=  1 

V ’ 


Fc  =0 
Gc  -0 
Hc  = 0 
Ec  = -k 


w = w 

wk  "T 


By  using  the  equations  of  matrices  F(T,t)  and  G(T,t)  with  these  data,  we  have 

[l-k^ft[l  - ft  - k-^-  0 -[1-k-^jft  £ 0 0 


0 0 


0 0 
0 0 


1 0 


(5-41) 


FP  F = 
xss 


[l-kT+kM^^]2  0 [l-kT+k2T-^~^-](l-kT)  0 0 


0 


[1-kT+k^T-^— -^](  1-kT)  0 


0 

0 


0 

0 


O-kT)6 

0 

0 


and 


Gwk6T  = 


[T-kT(T-T)]2  \ 


T0 


0 [T-kt(T-T)]  -rfL 


t a 


[T-kT(T-r  )y^L 
0 
0 


0 

0 

0 

0 


0 9 
a 2t^ 


w 


T 

0 


Therefore  by  equation  (5-38)  we  have 


0 

0 

0 

0 

0 


0 0 
0 0 
0 0 
0 0 


0 

0 

0 

0 

0 


(5-46) 


(5-46) 
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P 


xss 


[T-kt(T-T)]2ow2 

0 

[T-kx(T-T)]xow2 

0 

0 

T[l-{l-kT+k2x-iIlll}Z] 

T[l-{l-kT+k2x(T"T)}{1'kT}-1 

0 

0 

0 

0 

0 

[T-kx(T-x ) ] xo  2 

o 2t2 

w 

0 

w 

0 

0 

T[l-{l-kT+k2xiI^-}{l-kT}] 

T|>(l-kT)Z] 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

From  (5-32),  (5-33),  (5-34),  (5-40),  (5-42)  and  (5-43),  we  have 


H1  = C-k  0 0 0 

0] 

(5-48) 

H2  = [-k(  1-  I1)  0 

0 

0] 

(5-49) 

Ha-c-^  0 

0 

0] 

(5-50) 

H^F  = [-k[l-  -k  (T-x)  + k2xll^-]  0 

k2I 

2 

[l-killll]  0 0 

(5-51) 

Hb  = [ (k2(T-x ) -k2x  0 - (T- 

■t) 

0 0] 

(5-52) 

= -k(T-kx(T-x)] 

(5-53) 

°B  = "k(T~T) 

(5-54) 

From  (5-39)  and  (5-41)  we  have 


p* 

eAss 


[ 


[T-kx(T-x)]2  aj 
T[l-{l-kT+k2x (^}2] 


+ 


+ 


2[T-kT(T-T)]TO|)g  ^ 

T[M)-kT*  i|l(T-T)}l-kT~ ] T[kT(2-lT)] 

k2T2o  2 

— , — — (5-55) 


and 


«"■  upw 
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i2.  2 


eBss 


r,  2/t  , ,3,  (T-t)A  tT-kT(T-T)3‘ V , 

- [k  (T-x)-k  T A-  J.]  [ ] 

4 ni-n-kT+k^T-^y^r] 

2 3 iT-xl  3 [T-kT(T-T)]  to  ‘ 

- 2[k^(T-T)-kJT-LLli-]  kJT(T-x)[-  w 


T[l-{l-kT+k  x(T-x)}(l-kT) 

2n.f  1 -,2_  2 


„ _ .2  o Px ^ 7 (T-t )P[kT-l ]^o 

+ \<2t2  (t't)  r - w 1 + v2  1 W 

k ' 16  [T[kT(2-kT)]]  k T 


(5-56) 


PeAss  and  PeBss  are  Plotted  in  fi9ure  24  is  a function  of  x.  The  dia- 
grams to  the  right  of  each  plot  show  the  times  corresponding  to  the  value  of 
the  controller  output  used  to  calculate  and  eg.  The  results  obtained  for 
this  example  agree  with  intuition.  For  small  x,  the  variance  of  eA  is  small 
and  that  of  eg  is  large.  For  large  x the  complimentary  situation  holds. 

5.2  Variation  E:  System  Configuration  and  Dynamic  Equations 

Figure  25  shows  the  time  responses  of  yc-|  and  y^.  They  are  the  same 
as  the  basic  model,  the  aircraft,  actuator  and  sensor  dynamics.  Therefore, 


xp(t)  = *(t,t0)xp(t0)  + ♦1{t,t0)up(t0)  + 4>2(t,t0)wp(t0) 


and 


y (t)  = C x (t) 
•'p  P p 


(5-57) 

(5-58) 


In  this  variation,  the  input  of  the  plant  is  the  average  output  of  the 
controllers  and  the  output  of  the  plant  is  the  input  to  each  controller  at 
the  different  times.  Then 

VV  = + yc2(tk,] 


(5-59) 
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“ci'H1  = W ' SW  (5'60) 

Uc2<V,)  “ »p<  VT>  ‘ Cpxp<VT>  (5-6,) 

The  discrete-time  equations  for  controller  number  one  are: 

Vl<V*sc>  * Fcxcl<VSc>  + scucA>  (5-6?) 

ycl<Vfc>  5 Hcxcl<VSc>  * Ecuc)(tk’  (5-63) 

for  k = 0,1,...,  and  for  controller  number  two, 

xc2<VUltSc>  » Fc«cj(tk^«c)  * CcUc2(tk*'>  (5-64> 

yc?(tk+Tt5e)  = Hcxc2(V,*5c>  * Ecuc2(tk+,)  (5_65) 

for  k = 0,1,...  . 


From  figure  25  yc2(tk)  = yc2^tk'2+T+,Sc^  then  by  using  (5_65^ 

yc2(tk}  = yc2(tk’2+T+fc)  = Hcxc2(tk-2+T+fic)  + EcUc2(V2+t)  (5'66) 

Substituting  (5-61),  with  (5-66)  gives 

yc2<V  * »c?ttk-z*T+«c)  ■ Hcxc2(V2txt{c>  + VpVV^1  (5-67) 

From  figure  25  ycl(tk)  = ycl ( tk-l+$c)  hy  lJsin9  (5-63) 

ycl<*k>  a ycl<*k-,t4c>  ■ WVV  * Ecep,P,tk-,)  t5'68) 

Define  the  new  variables  *hcl.  *hpl.  \hc?.  «hc2.  «hhp?  and  *hp2  10 

have  the  variables  that  change  in  a piecewise  constant  manner.  There  are 
additions  that  are  required,  namely. 
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hcl(tk+1+V  " xcl(W 

(5-69) 

hpl(tk+1)  = xp(V 

(5-70) 

hhc2(tk+1+T+V  = xhc2(tk+T+6c) 

(5-71) 

hc2(tk+1+T+6c)  = xc2<VT+V 

(5-72) 

hhp2^tk+1+T^  = Xhp2^k+T) 

(5-73) 

hp2(tk+1+x)  = Xp(tk+T) 

(5-74) 

By  using  (5-69),  (5-70),  (5-71),  (5-72),  (5-73)  and  (5-74),  (5-67)  and  (5-68) 
become 

*c2<V  - *c2<V2tTtSc>  ' Hcxl,hc2<VTtSc>  * EcCpxhhp2<VT>  (5-75) 
»d(tk>  - *d<V'*V  ■ Hcxhd(V!c>  * EcCpxhpl(tk>  (5-76> 

Substituting  (5-75)  and  (5-76)  into  (5-59)  gives, 

W = + Vp'hpl^kJ  + Hcxhhc2(tk+T+6c)  + EcCpxhhp2(tk+T)] 

(5-77) 

At  time  t^- 1 +t+6c  the  value  of  Up  changes,  so  let  t=t|(-l+T+6c  and  tQ=t|<. 

Then  (5-57)  becomes 

xp(tk-l+T+6c)  = $(T+VT)xp(tk)  + ^1  (T+6c“T)Up(tk)  + 4'^T+VT)wp(tk)  (5_78) 

Substituting  (5-77)  into  (5-78)  gives 
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(t+«.-T) 

xp(tk-1+T+fic}  = ♦'  +VT)xp(tk)  + *1  2 Hcxhcl(W 

(t+6  -T)  (t+6  -T) 

+ 4*!  2 EcCpxhpl(  Vwp(  V + ^1  2 Hcxhhc2(tk+T+V 

(t+6  -T) 

+ *1—1 EcCpxhhp2(  Vt)  + *2(t+VT) 

At  time  tk+6c  the  value  of  Up  changes  again  so  let  t=tk+6c  and  tQ=tk-l+T+6 
then 

XP<W  = <t(T-T)xp(tk-1+T+6c)  + 4/1(T-T)Up(tk-l+T+6c)  + ^2(T-T)wp(tk 
Consider  Up(tk-l+T+6c) . From  figure  25 

Up(tk-'+™c!  = 4tycl(tk-1+T+sc)  * ’ 3*>cl(V,Mc> 

+ yc2(VHT+sc)] 

By  using  (5-61)  and  (5-65) 

yc2(tk-'+TMc)  = Hcxc2(tk-UT*5c>  * EcCpxp(tk-UT) 

By  using  (5-72)  and  (5-74),  becomes 

yc2<tk-HTMc)  - Hcxhc2(tk+,*Sc)  + EcCpxhp<;(tk*T) 

Substituting  (5-76)  and  (5-83)  into  (5-81)  gives 

V V1+t+6c)  = ^HcVl*W  + EcCpxhpl(V  + Hcxhc2^k+T+6c^ 

+ EcCpxhp2(tk+T):i 


(5-79) 


(5-80) 


(5-81) 


(5-82) 


(5-83) 


(5-84) 


Substituting  (5-79)  and  (5-84)  into  (5-80)  gives 
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t+6  -T 


xp(V«c)  = *(«c)xp(tk)  + )HC  + ^1(I?)HC](V«C) 


+ [*(T-t)v1>1  ( 


x+6  -T  / T 

~r~)EccP  + WVi  + 


♦l (t+6  -T) 

+ ♦(T-t)— — y— H x 


4»-«  (t+6  -T) 

\ I Li 


cxht,c2(V™c>  + t(T-T> 1 


+ *,(T-T)Hcxhc2(tktT«c)  t ♦v^T^cVV'V*1 
+ [♦  ( T-t  )'P2(t+<5c-T)  + ^(T-1)]  wpUk) 

Now  let  t=tk+l  and  t0=^|<+5c*  Then  (5-57)  becomes 
V V’>  * '>(T-6c)xp(tkrtc)  k *,(T-«c)Up(tk+5c)  * *2(T-«c)wp(tk) 
Consider  U (t. +6  ) 

P *>  C 

From  Figure  25 

W«c>  = + yc2(tk+6c)]  = ?CyCl^V6c) 


(5-85) 


(5-86) 


+ yc2<V1+T+6c>] 


By  using  (5-58)  and  (5-63) 


(5-87) 


(5-88) 


»cl<W  ■ Hcxd<Vsc>  * Eccpxp(*k> 


Substituting  (5-88)  and  (5-83)  into  (5-87)  gives 


Wsc>  - 2^c\\^W  * Eccpxp<tk>  + Hcxhc2<V”V  * Eccp’Wtk*T>  <s-89> 


Substituting  (5-85)  and  (5-89)  into  (5-86)  gives 
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(T-6  ) (t+6  -T) 

xp(tk+l)  = O(T)  + EcCp]xp(tk)  + *(T-«c)[*(T-t)^1 Hc 

f T T \ 'l'i(x+<sc-T) 

+ Hc]xhcl(tk+6c)  + ♦(T-gc)[»(T-T)  - /■  ■ EcCp  + 

d t i (t+6  -T) 

* t,11?1  EcCp]xhp1(tk)  * *{T-«c)*(T-t)*, gS — Hcxhhc2{tfc-«+ac) 

(t+5  -T)  fT-r) 

t *(T-«c)*(T-t)+, j^EcCp  xhhp2(tk«)  * t*(T-«c)+1Air-L  „c 


(T-«c) 


* [ ♦<T-«c>*lIEri  EcC[ 


+ ^^2^cCp]xhhp2(tk+T)  + ^1^2^CxC1(V6C)^(T-t)'I'2(t+VT) 


+ 4.2(T-t)]  + 4>2(T-6c)xh  2(tk+x)w  (tk) 


(5-90) 


Substituting  (5-58)  and  (5-60)  into  (5-62)  and  substituting  (5-58)  and  (5-61) 


into  (5-64)  gives. 


xcl(tk+1+6c)  = Fcxcl(tk+6c)  + GcCpXP(tk) 


xc2(tk+1+T+6c)  = FcXc(tk+T+6c)  + GcCpXp(tk+T) 


(5-91) 


(5-93) 


In  this  case  the  quantity  x (t.  +t)  can  be  written  using  the  solution  to  equation 

D K 


(5-91): 


(t~6  ) 


(t+6  -T) 


x (tk+T)  = C*(t)  + 4»i  ~ g - EcWV  + *(T-6c)|’*(T"T)'h f 


(t+6  -T) 


*, ^r1  Hc]xhc1(tk)  ♦ ♦(t-<c)[»(T-t)»1  f EcCp  * *,^cCp]xhp,(tk) 


* ♦(T-Sc)*(T-t)*)1^-  H£xhhc2(tk*-*fc)  * ♦(T.f,.)  ^cC( 
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i <T-V 


(--«_) 


EcCp]xkp2(VT>  * *1-1^-  Vcl'VV  * <*<T"0*?(t«c-T) 


+ ^(T-t)]  + \P2(T-fc)}  w (tk) 


(5-92) 


The  inherent  error  is  defined  in  two  parts  as  follows: 


eA(t»  - *ci<vsc>  ■ ycj(v,tSc> 


for  t|<+'''+'5;c  £ t < t|<+T+<sc»  k = 0,1,...,  <5c<t<T 


(5-93) 


eB(t)  = ycl(tk+Hfc)  - yc2(VT+6c) 


for  t. +1+6  < t < t.+l+t+6  k = 0,1,...,  6 <t<T. 


(5-94) 


Figure  26  shows  the  skewed  sampling  and  inherent  errors.  These  equations 
can  be  put  in  compact  form  by  writing  them  in  terms  of  a combined  stated  vector 


x(tj 


’xP(tk> 

xhpllV 

xcl<V5c> 

xhcl*W 

xhp2(VT> 

Xhhp2^k+T^ 

xc2<VtMc> 

xhc2<VT+Sc) 

>c2(VTt5c> 


xp<V> 
xhpl*V1* 
xcl'V'Wc> 
xhd<V,*T> 
X<tkt,)  = xhp2<VU,) 

Xhhp2(tk+1+T^  | 

xc2(tk+,+T+fic) 

xhc2(tk+1+T+(Sc) 

Xhhc2(tk+1+T+V 


l 


The  states  equations  become 
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l 

f r 

I 

r 

i: 
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”9nre 


x(tk+l)  = F(t,x)x(tk)  + G(T,x)wp(tk)  (5-95) 

where  F(T,x)  is 

(T-5J 

f„  - ♦(T)  + 2-S-  EcCp] 

-T)  (T  , 

f12  = ♦CT-«c)C*(T-t)*1 EcCp  * EcCp] 

(T-«c) 

f 13 = ♦it1-  Hc 


(T+VT)  IT-*) 

f14  = *(T-*c)[*(T-x)»1— -r—  HC  + HC3 


* 


r 


i 


, 


51 


(T-6  ) 

*(T)  + ^ EcCp 


52 


(T+fc-T)  (T-t) 

•(t-«c)C«(T-x)*1-t£—  eccp  ♦ EcCp 


(t-6c) 


53 


1 2 


H 


54 


(t+VT)  (T-t) 

^ Hc  + Hc 


55 


fT.T\  (T“0 

•I'-Wrr1  EccP  + *iV-  EccP 


56 


(t+6  -T) 

*(t-«c)*(T-t)i|.1 f EcC 


57 


= 0 


58 


H + 


*(T-«c>^r 2 "c  T T 


(t-O 


59 


(r+6  ~T) 

*(t-6c)*(T-t)^ £ Hc 


61 


f62  = f63  = f64  = 0 


65 


= 1 


'66 


f67  = f28  " f69 


71 


G C f,, 
c p 51 


72 


GcGp  fc2 


73 


G C f, 
c p c3 


74 


GcCp  f54 


75 


GcCp  f55 


76 


GcCp  f56 


77 


= F . 


78 


G C_  fcB 

C p 58 


79  = 


GcCp  f59 
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r 


N 


f87  = 1 


f88  = f89  = 0 


f98  = 1 

f99  = 0 
and  G(T,t)  is 

g1  = *(t-«c)[*(t-t)4»2(t+«c-t)  + vp2(t-t)]  + <i»2(t-«c) 

g2  = 0 = g3  = g4 

g5  = 4>(T-6c)[4>(T-t)4»2(t+«c-T)  + *2(T-t)]  + *2(r-«c) 

g6  = 0 

g7  = GcCp  g5 
g8  = 0 = g9 

The  controller  output  equations  are 

ycl(Vsc)  ’ V<V 


and 


then 


and 


yc2(VT+V  = H2x(V  + n 

H,  = [EC  0 H 0 0 0 0 0 0] 

I c p c 

H2=^cCpf51  EcCpf52  EcCpf53  EcCpf54  EcCpf55  EcCpf56  Hc 

EcCpf58  EcCpf59^ 

0 = EcCpg5 


(5-96) 

(5-97) 

(5-98) 

(5-99) 

(5-100) 


1 


„ - r 

' W 

1 

! i 

1 : 

11-99 

r 

b 

i t* 

The  equations  of  ( t ) and  eg ( t ) become 

«- 

eA(t)  = (HrH2)x(tk) 

= nwp(tk) 

(5-101) 

! 

for  tk+T+6c  < t < tk+l+6c,  k = 0,1,..., 

6c<t<T 

and 

| I. 

eB(t)  = (H1F-H2)x(tk 

+ ( H i G-o ) Wp ( tk ) 

(5-102) 

for  tk+l+6c  < t < tk  + l+T+<?c,  6c<t<T. 

| 

; 1 #■- 

5.2.1  Covariance  Analysis 

The  analysis  is  the  same  as  Variation  D.  In  this  variation  of  the  value  of 
time  delay  (6  ) is  equal  to  zero.  Then  this  variation  is  the  same  as  variation  D. 

i - ^ 

Therefore,  if  the  time  delay  (<5  ) is  equal  to  zero,  the  errors  of  this  variation 

c 

must  be  equal  to  the  covariance  errors  of  variation  D. 

I 

5.2.2  Example 


With  data  from  the  first  example  in  Section  II.  Section  F(T,t)  and 
G(T,t)  are: 


I 


I 
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.2  {t-6c}  k2 


hA  = t‘k"  2 ^ ' F 6 c 0 0 'k  ~~7 


? (T-O  ? <t+VT) 

-k  0 0 0]  (5-107) 


n = -kT 


(T-«  ) k\  .2  .2 

HtF  - C-k[l-k  — g-S-3  -fS-  o 0 Y (2T-t-«c)  f-  (t+6c-T)  0 0 0] 


k2  'k?6c 

hb  ■ tf  <T-'>  -T-- 


0 0 


-k‘ 


T 


(T-t)  0000] 


H|G  = -kT 


nB  = 'k^T"T^ 


By  using  equations  of  PgA  and  PpB  we  have 


4 (t-6  )2o  2T  4 (T-6  )(r-«  )o  2T  R4  (T-6  )2o  2T  k2T2o  2 

eAss  r KT(2-Ct) — F riF^TO-kTy]  F TOFFT]-  “T — (5-108) 


9 

i,**  9 o T T o I 9 /y  o 

PeBss  = T (T~t)  fkT(2-kT)  + l-(  l-kT)( 1-kTT  + kT(2-kT)]  + k T °w  (5_109) 


2o 


o 2T 


The  same  as  the  delay  model,  PeAss  depends  on  time  delay  6c  but  PepSs 
doesn't.  PgAss  and  PeBs$  are  plotted  in  Figure  27  as  a function  of  t.  The 
diagrams  to  the  right  of  each  plot  show  the  times  corresponding  to  the  values 
of  the  controller  outputs  used  to  calculate  eft  and  e^.  The  results  obtained 
for  this  example  agree  with  intuition.  For  small  t,  the  variance  of  eA  is 
small  and  that  of  e^  is  large.  For  large  t the  compl imentary  situation  holds. 


r 


i 

l 


P_ — —7- 

I 
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5.3  Variation  F 

The  same  as  variation  C,  there  are  two  cases  in  this  variation. 

1.  The  computation  of  y ^ completed  in  the  same  period  of 

the  computation  of  yc1 (t+5c<t) • 

2.  The  computationof  ycp  is  not  completed  in  the  same  period 
of  the  computation  of  ycp(T+fc<T)* 

5.3.1  Variation  F (Case  I):  System  Configuration  and  Dynamic  Equations 

Figure  28  shows  the  time-responses  of  ycl  and  yc?.  The  same  as  the  basic 
model,  for  the  aircraft,  actuator,  and  sensor  dynamics.  Therefore: 


xpU)  = *(t,t0)xp(t0)  + ^1(t,t0)Up(t0)  + H*2(t,t0)wp(to) 

and 

(5-112) 

yp(t)  ■ CpXp(t) 

(5-113) 

In  this  model,  the  input  of  the  plant  is  the  average  outputs  of 

the  controllers 

so 

Up(tk>  " ? tycl ttk)  tyc2(tk) 

(5-114) 

“cl<lk>  ■ W ■ cpxp(tk> 

(5-115) 

uc2(tk*T>  ' VVT>  ’ CpxP(VT) 

(5-116) 

The  discrete-time  equations  for  controller  number  one  are: 

xcl<Vmc>  = Fcxcl<V5c>  ‘W'k1 

and 

(5-117) 

ycl ( tk+<Sc J * Hcxcl(W  + GcCpW 

(5-118) 

for  k = 0,1,...,  and  for  controller  number  2, 

■■■ 
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xc2(tk+1+T+6c)  = Fcxc2(Vt+6c)  + GcCpxp(tk+T) 


(5-119) 


yc2(tk+T+5c)  = Hcxc2(tk+T+fc)  + EcCpxp(tk+T) 


(5-120) 


for  k = 0,1,..., 
From  Figure  28 


W "I^cl^kJ  + *c?W 


where 


= ycl(tk+6c)  = "cWV+V  + EcCpV  V]> 


(5-121) 


(5-122) 


yc2(tk5  = yc2(V1+6c+x)  = Hcxc2<V,+T+4c>  + EcCpxp(tk-1+T>  (5'123) 

Define  new  variables  xhcl,  xhpl,  xhc2  and  xhp2  as  in  type  E,  namely, 


xhcl^k+UV  = xcl(W 
xhpl(tk+1)  = V V 
xhc2(V1+x+5c>  = xc2(tk+T+fc) 

xhp2(V1+x)  = yvT> 

Then  the  equations  of  yc^  and  yc2  become 

ycl(V1+V  3 HcWW  + EcCpxhpl(V 


(5-124) 


(5-125) 


(5-126) 


(5-127) 


(5-128) 


yc2<V1+T+6c>  = Vh^V^V  + EcCpxhp2^tk+T^ 


Substituting  (5-125)  and  (5-129)  into  (5-121)  gives 
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V V = 2 ^cVl^VV  + EcCpxhpl(V  + Hcxhc?(tk+T+V  + EcCpxhp2(tk+T^ 

(5-130) 


at  t = t ^+6c , the  value  of  up  changes  or  let  t=tQ=tk.  Then 

xp<  VSc>  ■ ♦(VV‘k>  * *l(*c,up(tk>  * *2(JcVtk) 

Substitution  of  (5-130)  into  (5-131)  gives 

xp(tk«c)  ■ ♦(«c)xp(tk)  * EcCpxhp,(tk)  * «cxhc,(tk«c) 


(5-131) 


M5,-) 


<M«  ) 


Cpxhp2<V5c>  * -V-  Hcxhc2<VTtSc>  + WV*<] 


(5-132) 


at  t»tk«  the  value  of  up  changes,  or  let  t=t|c+T+l5c  and  t0=t|<+^c*  Then 

x (tk+T+5c)  = *(T)xp(tk+«c)  + Vt)uP(W  + Vt)VV  (5-133) 


From  Figure  28 


where 


up(tk+6c}  = 2 [ycl(W  + yc2(tk+6c)] 


ycl(tk+6c)  = "cWW  + EcCPxP(tk> 


(5-134) 


(5-135) 


yc2(tk+5c}  = yc2(V1+T+V  = Hcxhcl(tk+T+V  + EcCpxhp2(tk+T)  (5‘136) 
Substituting  (5-135)  and  (5-136)  into  (5-134)  gives 

V V!C>  ■ ? tHcxcl<V{c>  * EcVp^k'  + HcV2(VT*fc'  * Eccpxhp2<VT” 

(5-137) 


— r rr 
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Substituting  (5-132)  and  (5-137)  into  (5-133)  gives. 


'Mt) 


^(t) 


XpVxt5c>  - r*(T«c)  * -y-  EcC  ],  (tk)  * -V  Hcxcl<W 


+ *<T>  -V1-  EcCpXhpl(tk)  + ^)  “V-  Hcxhcl< W 


4>i($  ) 'Pi  (t)  v~  / 

+ [*(t)  1~  EcCp  + 2 EcCp]xhp2{tk+T)  + + *1  ~r~  Hc 


+ HcKc2(tk+T+6c)  + C*(t)  2(  V + '*'2{T)]wp(tk)  (5_138) 


The  equation  of  xp ( t^+1 ) can  be  derived  by  letting  t=tk+l  and  to=tk+T+<!;c.  Then 


xp(tk+1)  = f(T-T-6c)xp(tk+T+6c)  + vp1(T-T-6c)up(tk+T+(Sc)  + ^^"WV 

(5-139) 

Consider  up(tk+T+«5c) . 


From  Figure  28 

WtMc>  -•  \ 

Substituting  (5-118)  and  (5-120)  into  (5-140)  gives 

VVT+5C>  ‘ I tHcxcl<VSc>  * EccpxP(tk)  + Hcxc2(VT,sc> 

* EcCpxp(  Vt)] 


(5-140) 


(5-141) 


Substituting  (5-138)  and  (5-141)  into  (5-139)  gives 
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x (tk+l)  = {^(T-t-6c)[4.(t+6c)  + EC]  + 


(T-t-6  ) 

— H1-  eccp}  W 


(T-t-6  ) 

+ r* — £- 


(T-x-«c) 


(6C) 


(5-145 


yc2<VT+V  = Hcxc2(tk+T+4c)  + EcCp  (T>xp(tk}  + EcCp*r  2 EcCpxhpl 
* EcCp*i¥  Hcxhcl(tk+8c>  * VM"1 

+ EcCp*i¥  Hcxhc2<VTt8c>  * EcEp^2(T^wp(Ek^ 


<‘k» 


These  equations  can  be  put  in  compact  form  by  writing  them  in  terms  of  a combined 
stated  vector. 


I 

I 

I 

r 

1 

I 

I 

I 

I 

I 

I 


x(t. ) = 


xp(tk> 

Xp(tk+]) 

xhpl(V 

xhpl(tk+1) 

xcl{W 

xcl(tk+1+6c) 

xhcl(W 

,x(tk+l)  = 

Xhcl(tk+1+V 

Xhp2^k+T^ 

Xhp2(tk+1+T) 

xc2(tk+T+{c) 

xc2(tk+1+T+fc} 

xhc2^k+T+6c^ 

_ xhc2(V1+T+V 

The  state  equations  become 

x(tk+l)  = F(T,x)x(tk)  + G(T,r)wp(tk) 

where  F(T,t)  is 


(5-146) 


,T)  (T-t-S  ) (T-t-S  ) 

fll  * *(t-t-«c)C»(t«c)  * *,Yeccp]  * ♦, — j-^ccp  * *l-T-^ccp  <T> 

(S  ) (T-t-S  ) , , 

f,2  = t(T-Sc)*,^  EcCp  + j-*-  EpCp*,-if  EcCp 

f,3  - ♦(T-T-«c)t,Ill  Hc 

(6  ) (T-t-S  ) , , 

fI4  - «(T-SC)»,-^  Hc  + EcCp*,^i  Hc 


ii-m 


f61  = GcCp  f 51 
f62  = GcCp  f 52 
f63  = 0 
f64  = GcCp  f54 
f65  = GcCp  f55 
f66  ’ Fc 


f67  = GcCp  f57 

f71  = f72  = f73  = f74  = f75  = ° 


f77  = o 

and  G(T,t)  is 

(T-t-«) 

g]  = 4>(T-t-6c)[<1>(t)\P2(^c)  + 4»2(t)]  + ^2(T-t-6c)  + ^ ^ EcCp  2(x) 

92  = 0 

g3  = o 
g4  = o 

g5  = *2(t) 

96  = GcV2(t) 
g7  = 0 


The  controller  output  equations  are 


ycl(V«c>  - H,x(tk) 


(5-147) 


then 


yc2(tk+T+fic)  = H2x(tk)  + 0 


(5-148) 


and 


" ■ EcCp  2(t)  (5-151) 

The  same  as  variations  D and  E 

eA(t)  = (HrH2)x(tk)  - owp(tk)  (5-152) 

for  t|<+T+<sc  1 t < t|(+l+|5c,  k = 0,1,...,  0<t<T-6c,  and 

( 

eB(t)  = (H1F-H2)x(tk)  + (H1G-o)w  (tk)  (5-153) 

1 

for  tk+l+6c  < t < tk+l+T+6c,  k = 0,1,...,  0<t<T-6c> 

Figure  29  shows  the  skewed  sampling  and  inherent  errors. 

5.3.1  Covariance  Analysis  and  Example 

The  analysis  is  the  same  as  that  for  variation  D. 

5.3.2  Example 

j 

With  data  from  the  first  example  shown  in  section  II,  the  matrix  F(T,t) 


and  G(T,t) 
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I i 


vrt8c  V 


CHANNEL  1 


-+-  CHANNEL  2 

V2  y<:2 


ftAm=  *ci(V8c,-»c2(Vt*8c) 
FOR  tk*T+Sc  < t <th  + l*Sc 

A (t)  = yCi  (tk+,*8c)  - *C2  (tk+T<'8c) 
FOR  Lk  + I+Sc  < t<t  <tk  + l ♦T48c 


k = 0,1,  ....  0 < r <T-S 


k = 0 ,1,  ...  0 < r < T - 8 


FIGURE  29  SKEWEO  SAMPLING  AND  INHERENT  ERRORS 
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,2T  .2t  k«.  -k«r  .2t 

(h1f-h2)  = [yl  + V + T^-  + V(t-t-sc)]  0 0 

-k«c  k2T 

■ kC~r£  + V{T"t"*c)]  ° 0] 


n = -kT 


H G = -k ( T - (T-x-6  ) 


H-iG-o  = -kT  + V-  fT-T-f  ) + kT  = -k(T-x)  + V-fT-x-*  ) 


2,4  2 


■ V rr  tT-  r <T-^»  . 

A Lt  O ’ T 


kVo  2 


eAss  4 lT  . 2t  , 

l-[  l-kT+  — T-t-5c)  j 


j-  ( l-( 1-kT)  j ] + 


2 7 7?  k^T  le^T  k ^ r le^T 

4 [T+6  -rHT-kT  2(T-x-6r)r  ° k%  - ~(T-x+f  ][V  + 

= * — £ £_  ji_  + ?r c 2 , c ? ^ <L 


eBss  T-  I _r -i  bT+  k^T/j  T k^T 

L ' 2 U VJ  1 -[  1 -kT  + ^p(T-x-fc)][l-kT] 

• r TO4  2 

[T  . |I(T-t-«c)] 

k26  . 3t  o a 2T2  . 2t  , o 2T 

[-  __c  , jS_I(T-x-6c)f  JL-.  + [-k(T-x)  + iy(T-x-6c)]2  -f- 

l-[ 1-kT]2 


kl  .3, 
c k 1 


PeAss  and  PeBss  are  Plotted  in  Figure  30  as  a function  of  x.  The  diagrams 
to  the  right  of  each  plot  show  the  times  corresponding  to  the  values  of  the 
"f roller  outputs  used  to  calculate  e^  and  e^.  The  results  obtained  for  this 
• agree  with  intuition.  For  small  x,  the  variance  of  e^  is  small  and 

,s  large.  For  large  x the  complementary  situation  holds. 


WHERE 


FIGURE  30  P 
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5.3.3  Variation  F (Case  II):  System  Configuration  and  Dynamic  Equations 


Figure  31  shows  the  time  responses  of  ycl  and  y^.  From  this  figure, 
xp(tk+x)  must  be  greater  than  xp(tk-l+T+tfc)  because  <$C<T.  For  example  if  t=0 
then  xp( 7 +x+<5c ) becomes  xp(tk-l+fc)  which  is  less  than  xp(tk).  If  we 
compare  this  figure  and  Figure  25  of  variation  E.  These  time  responses  are 
identical  except  that  the  position  of  xp  is  tk+T.  So  we  can  use  all  the 
equations  from  type  E except  the  equation  of  Xp(tk+T)«  Therefore,  from 
variation  E,  the  equations  of  Xp(tk-1+T+5C)»  xp(tk+‘f'c^  and  xp(tj<+^)  are: 

(T+«  -T)  (t+«  -T) 

*p(Vl+r+«c)  - •(*+VT)xp(tk)  + *1 2“  Hcxhcl(W  + *1—7 ■EcCpxhpl(tk) 

(t+6  -T)  (t+6  -T) 

+ 4*!  5 Hcxhhc2(tk+T+V  + *1  2 EcCpxhhp2^k+T^ 


+ VT+VT^VV 


(5-154) 


(t+6  -T) 


^ l 1 / it  T1 

(V«J  = ♦(Sjxjtj  + [♦(T-t)^ f-  Hc  + V2V  xhc1(tk+«c) 


pv  k c;  k 


(t+VT)  (T-t) 

+ [<i>(t-T)^1  | eccp  + *1nrj€ccp]xhpi(tk) 

(t+6  -T)  (t+6  -T) 

+ *(T"t)4,1 f-«cxhhc2(VT+fc)  + *<T-t>*1 \ EcCpxhhp2(tk+T) 


+ ♦l<T-T)HcxhC2<tk+T+6c)  + '^l1IrlEcCpxhp2(tk+T) 


+ [4>(T-t)4»2(t+(Sc-T)  + vP2(T-t ) ]w  (tk) 


(5-155) 
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(T-«  ) (t+«c-T) 

V V1J  = [*(T)  + ^-^-EcCp>p(tk)  + •(T-«c)C*(T.t)*1— g£ Hc 

* Hc>hcl(tkt{c>  + «(T-«c)t»(T-T)*,— ^2ecCp 

(T-t)  (t+«  -T> 

* Wxhpl(tk>  * <'(2T-5c-T>*l-7£—Hcxhhc2(tktTtfc) 

+ «(2T-«c-t)+1— EcCpxhhp2(tk*T)  * [♦  (T-«c)*,^c 

* + WT-V*!11?1  EcCp 

(T-6J  (T-«_) 

+ 'h~2 EcCp]xhp2(tk+T)  + ^1 — 2 HcXcl(tk+6c) 

+ {*(T-6c)[*(T-tH2(t+«c-T)  + xP2(T-t)]  + 4»2(T-«c)}wp(tk)  (5-156) 


The  same  as  variation  F (case  I)  the  discrete-time  equations  for  controller 
number  one  are: 

xcl<V>*sc>  * Fcxcl<W  * GccpxP(tk>  (5-157) 

*cl<  W ’ Hcxcl(tk*sc>  + Eccpxp(tk)  (5-,58) 

for  k = 0,1,...,  and  for  controller  number  2, 

xc2(tk+1+T+6c)  = Fcxc2(tk+T+V  + GcCpxp(tk+T)  (5‘159) 

yc2(VT+V  = Hcxc2(tk+T+fc)  + EcCpxPUk+T) 


for  k * 0,1 


(5-160) 


ppF"»wyi»i»i  « L IIMPWPHH 
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Let  t=tk+r  and  t0=tk-l+T+6c  then  by  using  (5-112) 

Xp(tk+T)  =<D(T-6c)xp(tk-l+i+6c)  +'h(T-Sc)xp(tk-l+T+6c)  + (5-161) 

Substituting  (5-79)  and  (5-84)  into  (5-162)  gives 

(t+6  -T)  (T-*  ) 

xp(tk+T)  =*(x)xp(tk)  + [4>(T-6C)4»1— Hc  -Hc]  xhci(tk+6c) 

(t+«  -T)  (T-«.) 

+ [*(T.«c)^_^_  Eccp  + ^_c_  EcCp]  xhpl(tk) 

(t+«  -T)  (t+6  -T) 

+ *(T~6c)'pl r Hcxhhc2(tk+T+V  + ♦CT-«C>*1 r~EcCpxhhp2(tk+T) 

(T-6J  (T-«  ) 

+ *1-7“  Hcxhc2(tk+T+6c)  + ^1  2 ^ cCpxhp2^tk+T^ 

+ [*(T-«c)<I»2(t+6c-T)  + *2(T-«c)]wp(tk)  (5-162) 


The  inherent  error  is  defined  into  two  parts  as  follows: 

W1’  - W - T-Sc^t<{c  (5-,63> 

for  tk+T+Sc  ^ t < tk+l+fc»  k = 0,1,...  . And 

eBl(t)  = yCl(V1+6c)  - yc2(tk+T+6c)  T'V^c*  (5-164) 

Figure  32  shows  the  skewed  sampling  and  inherent  errors. 

These  equations  can  be  put  in  compact  form  by  writing  them  in  terms  of  a 


combined  stated  vector. 


CHANNEL  t 

V'*8C  V2  <,c,) 


L|j*l  Lk»T  + 8c 


s*2 


CHANNEL  2 

nct) 


»*(t)  = yCi « V8C>  * Vet  (tk*r<.8c) 

FOR  tk  »t»8c  < t < tk«  I *8C  k : 0 , 1 , . . . , T - 8C  < r < 8C 

o e it)  = yc,  <tk*u8c)  - yctitk*T*8c ) 

FOR  tk*u8c<t<tk4l*r*ic  k = 0,1,..,  T-8C<T<8C 


FIGURE  32  SKEWEO  SAMPLING  AND  INHERENT  ERRORS 


frit 


x(tj  = 


xp(tk> 

xhpl(tk) 

xc)<Vsc> 

xhcl^tk+<c* 

xhp2(VT> 

xhhp2<VT> 

xc2(VtWc> 

xhc2(VT*5c> 

_ xhhc2(tkMWc> 


,x(tk+l) 


Xp(lk+1) 

xhpl(V  * 

xcl(V1+lc' 

xhcl(tk+1+5c) 

xhp2(tk+1+x) 

Xhhp2^k+1+T) 

xc2(V1+t+6c) 

Xhc2(tk+1+X+V 

xhhc2(tk+1+T+fc) 
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The  state  equations  become 


x(tk+l)  ■ F(T,T)x1(t|<)  + G(T.t)w  (tk) 


(5-165) 


where  F(T,t)  is 

fn  = *(T)  + 1|> 


(T-«J 

c r r 
1 2 cLp 

(x+6  -T)  ,T 

f„  = «(T-<r)[*(T-T)»1-  / - EcCp  * ♦,^-p  EcCp] 


^10  = * 


(T-6  ) 


13  " V1 2 "c 


(t+S  -T)  tT~) 

f14  » ♦(T-«c)[*(T-t)4»1 § Hc  + Hc] 

(T-t)  (T‘fc) 

f15  = *(T-6C)*1  2 EcCp  + EcCP 


fi6  = *(2T-fc"T)4,r 


(t+«  -T) 


f17  30 


f18  = *(T-«c)*,iIri  Hc  ♦ +T 


(T-«J 


f19  = *(2T-«  -T)*r 


(t+«  -T) 


4 


f23  = f24  = f25  = f26  = f27  ' f28  " 0 


f35  = f36  = f37 
f42  = 0 


f38  = f39  = 0 


f45  = f46  = f47  = f48  = f49  = 0 


«(T-«  )Mt+«  .T)EC  *(T-UECr 
£ £JL  + vp? 


f53  = 0 


♦(T-«r)vP  (T+6  -T)H  4>(T-5  )H 
f — c I c c , «|i  c c 

T54  2 12 


fee  = * 


(T-«J 


’55  T1  2 c"p 

^ »(T-5c)^(t+«c-T)EcCc 

'56  = 2 


f57  = 0 


f58  = V™c>  Hc 

, *(t-«c)Vt+Vt)  hc 

f59  = 2 


f61  * f62  = f63  a f64  = 0 


11-123 


1 
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The  controller  output  equations  are 


»cl(V5c>  = H1X<V 


(5-166) 


= H2x(  V0) 


(5-167) 


H]  = [EcCp  0 Hc  0 0 0 0 0 0] 

H2-tEcCpf51  EcCpf52  0 EcCpf54  EcCpf55  EcCpf56  Hc 

EcCpf58  EcCpf59^ 


(5-168) 


(5-169) 


0 = EcCp  95 


(5-170) 


Using  the  same  variations  D and  E, 


;A(t)  = (H1-H2)x(t|<)  - owk(tk) 


(5-171) 


for  tk+^+6  <t<tk+l+6c,  k = 0,1,...,  t-6c<r<6c,  and 


?B(tk)  = (H1F-H2)x(tk)  + (H^-nJw  (tk) 


(5-172) 


for  tk+l+6c  < t < tk+l+T+5c>  k = 0,1,...,  t-6c<x<6c. 


5.3.4  Covariance  Analysis 


The  analysis  is  the  same  as  in  variation  D. 

5.3.5  Example 

With  the  data  from  the  first  example  of  the  basic  model,  the  matrix  F 
can  be  calculated  as 


From  (5-171),  (5-172)  and  (5-173)  we  have 


H]  = [-k  0 0 0 0 0 0 0 0 0] 

2 2 2 

H2  * C-k  *_I  o 0 |-(T-6C)  |-(t+6c-T)  0 0 0] 


and 

2 2 2 

(HrH2)  =[0^0  0 ^-(T-6C)  ^-(T+fc-T)  0 0 0] 


k^6  2 2 

I HjF  = [-k[l-  |(T-«C)]  -y^  0 0 |-(2T-t-«c)  y-(T+6c-T)  0 0 0] 


2 2 2 

hif_h2  = [r(T-6c}  r(VT)  0 0 r(T"T)  0 0 ° °3 


L 


HjG  = [-kT] 
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H^G-o  = -k(T-T) 


k4  2 a 2j  W 

PeAss  = 4^T_<V  kT^-kx)  + T 


k4  <T-«C>VT  . k4  (T-)(T-«c)aw2T  _ k4  (T-^)V-  _ k2(T-x)2aw2 


3 _ K C'  W 

eBss  " 4 kT(2-klT 


+ 4 kT(2-kT)  + 


PeAss  and  PeBss  are  P1otted  in  figure  33  is  a function  of  t.  The  diagrams  to 
the  right  of  each  plot  show  tie  times  corresponding  to  the  values  of  the 
controller  outputs  used  to  calculate  eft  and  eg.  The  results  agree  with 
intuition  just  as  in  Case  I.  For  small  t,  the  variance  of  e^  is  small  and 
that  for  eg  is  large.  For  large  •r,  the  complimentary  situation  holds. 


h4(T-8e  )*o£8( 


T|2MT-8C  )' 


FIGURE  33  P, 
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6.0  SOFTWARE  FOR  MODELING  WITH  COVARIANCE  AND  TRANSIENT-RESPONSE  ANALYSIS 


As  in  Reference  1,  the  software  consists  of  a single  FORTRAN  main  program 
and  16  subroutines  for  all  models  (Multirate,  Delay,  Output-Average). 


6.1  Flow  Chart  of  the  Delay  Model 


A flow  chart  of  the  program  of  the  Delay  Model  is  shown  in  Fig.  34.  The 
first  block  shows  the  data  input.  The  variables  are  self  explanatory  except 
for  the  quantities  NT,  NTAU,  NTIME,  and  NWRITE.  Since  a numerical  integra- 
tion is  required  to  compute  V-j(~),  V^T),  V^(T-^),  V-|(T-fc),  Vj(£c)  and 
Vi(T-f  ) where 

t 

V,(t)  = / 4>(t,s)ds 


It  is  necessary  to  quantize  the  time  interval  [0,T]  and  the  user  specifies 
the  degree  of  quantization  by  supplying  NT,  the  number  of  subintervals  in 
[0,T]  which  are  to  be  used  in  the  computation.  The  user  specifies  the  number 
of  evenly-spaced  values  of  t and  f by  providing  the  number  NTAU,  and  NDELAY; 
the  values  of  t and  6c  are  then  computed  within  the  program  as 

NTAU*T 


DELAY 


NDELAY*T 


NTIME  determines  the  number  of  T-incremented  points  to  be  computed;  that 
is,  compute  transient  data  for  the  following  values  of  time. 

TIME  = 0,T NTIME*T 

NWRITE  determines  the  values  of  TIME  for  which  the  data  is  to  be  written, 
that  is,  the  time  response  data  is  to  be  written  every  NWRITE*T  seconds 


Read  in  heading  and  description  of  the  plant  and  the  controller; 

N » N , N . N , An>  B,  , 8,  , C , F , G . H . E , T,  NT, 
p up’  wp  op  c p lp  2p  p c c’  c’  c 

NTAU,  NDELAY,  W,  NTIME,  NWRITE 


Calculate  4>(x),  <t>(6c),  4>( T) , <t>(T-6c) , 4>(t-6c),  ^(t),  ^(t-6c) 
^(T),  \l>i  ( 6 ) and  ^(T-6  ) 


Calculate  ^(t),  i|i^(t-6  ),  ^(T),  ^(s  ) and  ^ ( T- 6C ) 


Form  F(T,t) 


Form  G(T,t) 


Cakulafe  Pxs$,  Pefls$  and  PeBs5 


Calculate  transient  time  response  of  the  states,  controller  outputs, 
and  the  inherent  errors  to  a unit  step  input 


Figure  34.  Flow  Chart  Describing  the  Major  Computations 
of  Program  Skew 
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O.NWRITE  T Sljfr 

The  second  block  specifies  the  calculation  of  f ^ ) % <j>(T),  $(T-<*c), 

■*>(t-'5c  ) . 4»i(T).  't'l  ^ T-^c ) % ^CO.  ) . and 

The  third  block  specified  the  calculation  of  ^ ( i ) , 4>p  ( ^ c ) , ^(T), 

♦2(t-«c).  am1  ^(T-fc). 

In  block  four  F(T,’),  the  discrete-time  state  transition  matrix  for  the 
overall  system,  is  computed  and  written.  As  in  block  four,  block  five  GIT,-1) 
is  computed  and  written. 

In  block  six,  the  steady-state  covariance  of  the  states  (P  ),  P{?  , 

and  Pegs$  are  computed  and  written.  The  technique  used  is  the  same  as  in 
Reference  1. 

The  final  set  of  computations  is  given  in  block  seven.  The  transient  time 
response  of  the  states,  the  controller  outputs,  and  the  inherent  errors  to  a 
unit  step  input  are  computed. 

6.2  Instructions  for  Using  the  Program 

As  in  Reference  1,  let  us  define 

NPM  = maximum  number  of  plant  states 

NUPM  = maximum  number  of  plant  control  inputs 

NWPM  * maximum  number  of  plant  disturbance  inputs 

NOPM  * maximum  number  of  plant  outputs 

NCM  3 maximum  number  of  controller  states. 

Currently  they  are  assigned  the  values  A,  2,  2,  3,  and  2 respectively.  The 

two  dimension  statements  in  the  main  program  must  use  these  same  numerical 

values  in  the  arrays.  Table  2 shows  what  values  to  assign  to  the  given  arrays. 

In  this  table,  the  numerical  values  of  NHM,  NFM,  and  NRRM  are  calculated  from 
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TABLE  2 

Required  Dimensions  of  All  Arrays 


f 

i 

i: 

i 

1. 

i 

i: 

i: 

i 

i: 

i: 

t 


AP(NPM,NPM) 

B1P(NPM,NUPM) 

B2P( NPM.NWPM) 

CP(NOPM.NPM) 

FC(NCM.NCM) 

GC(NCM.NOPM) 

HC(NUPM.NCM) 

EC(NUPM.NOPM) 

ECCP(NUPM.NPM) 

GCCP(NCM.NPM) 

PHITAU(NHM.NHM) 

PHIT(NHM,NHM) 

PHTATD(NHM.NHM) 

PHTDEL(NHM.NHM) 

PSIT(NHM.NHM) 

PSITAU(NHM.NHM) 

PS?TAT(NHM,NHM) 

PS1TTD(NHM,NHM) 

PSITDE(NHM.NHM) 

PS2TAT(NHM,NHM) 

PS2TTD(NHM,NHM) 

PS2TDE(NHM,NHM) 

PS2T(NHM,NHM) 

PS2TAU(NHM,NHM) 

F(NHM,NFM) 

AD(NHM,NHM) 

DD(NHM) 

DELPHI (NHM.NHM) 

HA(NUPM.NFM) 

HB(NUPM,NFM) 

El(NUPM) 

D2(NUPM) 

X(NFM) 

XW(NFM) 


EIF(NHM,4) 

El V(NHM,4) 
EH(NHM.NHM) 
FHST(NHM.NHM) 
KWA(NHM) 

RF(NRRM) 

RR(NRRM) 

I D ( 20 ) 

PT(NHM.NHM) 

INDEX(NHM) 

W(NWPM.NWPM) 

PS(NHM,NHM) 

V(NFM,NFM) 

PXSS(NFM.NFM) 

PEAXSS(NUPM.NUPM) 

PEBXSS(NUPM,NUPM) 

Vl(NPM.NPM) 

VITAU(NPM.NPM) 

VITDEL(NPM.NPM) 

V1TATD(NPM,NPM) 

VITTD(NPM.NPM) 

HI (NUPM.NFM) 

H2(NUPM,NFM) 

AM(NFM.NFM) 

PM(NFM.NFM) 

D(NFM,NFM) 

DW(NFM,NFM) 

PL(NUPM.NUPM) 

YCl(NUPM) 

YC2(NUPM) 

HPH(NUPM,NUPM) 

GWG(NFM.NWPM) 

PLB(NUPM.NUPM) 

G(NFM,NWPM) 


1 


I 


l 

L 

[ 

l 
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NHM  = 2+NPM+NUPM 
NFM  = 2+NPM+3*NCM 
NRRM  = 4+NPM+4 

There  are  two  listings  and  three  examples  of  output  (A-7D  pitch  axis 
flight  control  systems)  for  the  variations  B and  C in  Appendix  A.  The 
flow  chart  in  Figure  34  can  be  used  for  both  of  the  variations  except  that 
for  the  equations  of  F(T,t),  G(T,t),  Pxss,  P.,Ass,  and  PeBss-  For  variation 
B use  the  equations  F(T,t),  G(T,t),  Pxss,  PeAss,  and  PgBss  in  section  4.2 


and  for  variation  C use  the  equations  of  the  matrices  found 


in  section  4.3. 
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7.0  SUMMARY 

As  mentioned  in  Section  1,  there  are  three  models  which  are  developed 
from  the  Basic  Model. 

1 . Multirate  Model 

2.  Delay  Model 

3.  Output-Average  Model. 

The  Multirate  Model  allows  for  separate  sampling  rates  for  the  external 

input  and  the  digital  controllers.  If  n in  this  model  is  equal  to  1,  then  the 

Multirate  Model  is  the  same  as  the  Basic  Model  except  that  the  external  input 

is  sampled.  The  data  from  the  first  example  of  the  Basic  Model  is  applied  to 

study  the  characteristics  of  the  covariance  errors  of  this  model  and  the  limit 

of  the  skew  time  is  0o<  ^ . As  discussed  in  the  section  of  the  Multirate  Model, 

the  covariance  errors  are  the  approximate  value  and  they  are  less  than  the  true 

value.  The  characterist ics  of  the  covariance  errors  follow  the  fact  that, 

they  are  largest  when  the  times  at  which  ycl  and  yc?  are  farthest  apart. 

^ The  Delay  Model  allows  for  computational  delays  due  to  the  time  required 

for  data  conversions  and  controls  output  computations.  There  are  three 

variations  A,  B,  and  C.  In  variation  A,  there  are  three  time  delays 

(6  , <5  and  .)  which  the  sum  of  them  is  less  than  the  skew  time  and  the 
a c d 

computations  of  yc^  and  ycp  are  completed  in  the  same  period  (t^,t^+l).  The 
limit  of  t is  AocT-A  where  A is  the  sum  of  the  delay.  In  variation  B,  there 
is  one  time  delay  ($c)  which  is  less  than  the  skew  time  and  the  computations 
of  yc^  and  yc,  aren't  completed  in  the  same  period.  The  limit  of  ▼ is  ^c<t<T. 

In  variation  C,  there  is  also  one  time  delay  (6  ) which  is  greater  than  the 
skew  time  T and  there  are  two  cases  for  the  computations  of  yc^  and  yc->. 

In  Case  I,  the  computations  of  y^  and  ycp  are  completed  in  the  same  period 
and  in  Case  II,  the  computations  of  yc^  and  aren't  completed  in  the  same 
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period  and  the  computation  of  y ^ are  completed  in  the  following  period. 

There  are  no  difference  for  the  derivation  of  these  two  cases  because  the  in- 
put of  the  plant  is  the  first  output  of  the  controller.  The  difference  is 

the  value  of  the  sum  of  " and  6 . The  limit  of  the  skew  time  are  0<T<f  . 

c c 

The  data  from  the  first  example  of  the  Basic  Model  is  applied  to  these  varia- 
tions and  the  covariance  errors  follow  the  fact  that,  they  are  largest  when 
the  times  at  which  y^  and  y £ are  farthest  apart.  The  values  of  the  covariance 
errors  are  approximate  value  and  less  than  the  true  value  which  are  the  same 
as  the  Multirate  Model.  Only  one  of  the  expressions  of  the  covariance  errors 
depends  on  the  time  delay.  This  is  because  the  data  of  the  first  example  is 
the  special  (Ap=0).  In  the  general  case,  both  of  the  expressions  of  the 
covariance  errors  depend  on  the  time  delay. 

The  Output-Averaging  Model  provides  for  averaging  the  control  outputs 
produced  by  each  of  the  redundant  controllers.  There  are  three  variations, 

D,  E and  F.  In  variation  D,  there  is  no  time  delay  and  the  limit  of  the  skew 
time  is  0<t<T.  In  variation  E,  there  is  one  time  delay  6 which  is  less  than  t. 
The  same  as  variation  B,  the  computations  of  ycl  and  y ^ are  not  completed  in 
the  same  period  and  the  limit  of  the  skew  time  is  fc<x<T.  In  variation  F,  which 
is  the  same  as  variation  C,  the  time  delay  6 is  greater  than  t and  there  are 
two  cases,  Case  I,  the  computations  of  y^  and  y £ are  completed  in  the  same 
period  and  in  Case  II  they  aren't.  Because  the  input  of  the  plant  is  the 
average  of  the  outputs  of  the  controllers  so  the  deviation  of  these  two  cases 

are  different.  The  limit  of  t of  the  first  case  is  O-o'cT-f  and  the  limit  of 

c 

T of  the  second  case  is  T-<$c<t<£c.  The  data  from  the  first  example  of  the 
Basic  Model  is  applied  to  these  variations.  The  covariance  errors  follow  the 
fact  that,  they  are  largest  when  the  times  at  which  y^  and  y £ are  farthest 
apart. 
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Section  6 describes  the  software  which  has  been  developed  to  apply  the 
methods  to  realistic  examples.  A general  description;  a flow  chart  and 
user  instructions  are  given  for  the  Fortran  program.  The  computer  program 
listing  and  the  example  of  the  A-7D  pitch-axis  of  the  Multirate  Model  are 
given  in  Appendix  8.  Appendix  A discusses  and  calculates  the  covariance  of 
the  sample-hold  of  the  white  gaussian  noise. 
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Appendix  A.  The  Correlation  Function  of  the  Sample  and  Zero-Order  Hold  of 
White  Gaussian  Noise 

The  white  noise  is  defined  as  a stationary,  zero-mean  gaussian  process 
with  power  spectrum 

a N 

Vf)  = T f<°° 

where  the  dimensions  of  NQ  are  watts  per  cycle  per  second,  or  joules. 

Actually,  white  noise  (whether  Gaussian  or  not)  must  be  fictitious  be- 
cause its  total  mean  power  would  be 

nj(t)  = / «w(f)  df  = - (C-2) 

which  is  not  meaningful.  Therefore,  we  can't  sample  zero-order  hold  white 
noise.  There  is  the  process  which  is  equivalent  with  sample  zero-order  hold 
process  from  the  fact  that  if  we  pass  the  white  noise  through  a linear  filter 
for  which 

/|H(f ) | 2 df  < °°  (C-3) 

-00 

produces  at  the  filter  output  a stationary,  zero-mean  noise  N(t).  Then  we 
have 

fn(f)  = ^|H(f)|2  (C-4) 

and 

N*(t")  = — jfl  H(f ) 1 2df . (C-5) 

-oo 

which  from  equation  C-3,  is  finite. 

Figure  1 is  the  block  diagram  of  the  equivalent  process  for  sample- 
zero  order  hold  process.  The  filter  output  is  the  average  value  of  the  noise 
in  the  interval  T (sample  period). 


r 
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Figure  1 


Let  z(t)  be  the  filter  output,  then 


Nw(t)  * h(t) 


N(t-  X)h(  X)d  X 


= ~ J N(t-  >)d  X 


The  correlation  function  of  the  filter  output  is 


E[Z2(t)]  = E J J Nw(t-x)hw(t-p)dXdu] 


= ■7  j I R (X-p)dXdu 
0 0 " 


where  Rn(t)  - E[N  (t)N  (t+-r)]  and  from  equation  C-l, 
w w 


Rw(t) 


"z(0)  = -7  / / 5(X-p)dXdP=;  J dXdp 
T o 0 ^ o Z 


7 iT 


Therefore  the  correlation  function  of  the  sample  zero-order  hold  white  Gaussian 
N 

noise  at  t=0  is  • 


noise  at  t=0  is 
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IS  BEST  QUALITY  PRACTICABLE 
raOM  COPY  PURWISHED  TO  ODC 

COMMON  EL  EMX,  MAX  I » MAX J 
DIMENSION  API  4.4)  . B I P 1 4 • 2 ) • 82PI 4 • 2 ) • CP  I 3 *4 ) • FC I 2 • 2 I • 

1 HC( 2 ,2) . EC(2 ,3 ) ,PHIT ( 10, 1 O) ,PS1T( 1 0. 1 0 ) ,F ( 14. 14 ) . 

? DELPHI (10.10).EIP(10*4).EIV(10.4).FH(10.10), 

3 KWA( 10) .RF(20) ,RR( 20) .10(20) ,FCCP (2 .4 ) , GCCP! 2,4 ) . 

4 INDEX! 1 0 ).N(2,2).PS(tO, 10 ) . PS  1 TATI  10. 10) ,GC(  2.3) . 

5 PS1TTDI  10.10)  .PMTATOC  10.  I 0 ) .PhTOEL  UO.IO).AH(IO.IO). 

6 PMTTOI 10 . 10) . VI  I 4.4) .VI TDELI 4, 4) .VITA TO (4. 4). 

7 PS2TDEI AO. 10 ) .PM I 1 4 • 1 4 ) . PS2T T D( 1 0 . 10).V1TT0(4.4). 

8 AM (14,14), G!  14. 2). 0114, 14). OW (14. 14). PL !2, 2). 

9 HI (2,14), M2 (2,14).HA(2.  14). HD I 2 • 14 ) . PL  8(2.2) , 

A PL WPL! 2.2) . PEAXSSI 2 . 2 ) • PEBXSS (2.2) .PXSSI 14.14). 

B FH3T! 10,10 ).DD! 10) .PT( 10 . 1 0) . PS1 TOE( 1 0 . 1 0 ) . 

C PH(TAU<|0,10),PS2TAT(10.10).GWG( 14.14). 

D XI 14) ,XW! 14) .El (2) .E2(2) .EAI2) .Ell (2) .E22I2) . 

E VC1 (2 ).YC2!2 >.PS1TAU( 10. 1 0 ) . V 1 TAU( 4 . 4 ) .HPH( 2 . 2) 

CCCCC  PROVIDE  MAXIMA  POM  CALLED  ARRAYS 
NPM  * 4 
NUPM  = 2 
NWPM  = 2 
NOPM  - 3 
NCM  * 2 

NHM  « 24NPM  4 NOPM 
NPM  * 24 NPM  4 3*NCM 
NRRM  * 44NPM  4 4 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CCCCC  RCAD  INPOT  DATA  CCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CCCCCCCeCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
•RITE (6.2222) 

2222  FORMAT!*!*) 

READ  IS. 900 ) ID 
9 00  FORMAT (20A4 ) 

WRITE !6. 902 ) ID 
902  FORMAT! • 1 • ,20A4> 

READ ! 5 , 906 ) NP . NOP. NMP. NOP . NC 
906  FORMATISl  3) 

HRITEI6.908)  NP. NOP, NWP . NOP . NC 
908  F QRMA  T ( * ONO.  OF  PLANT  STATES  =*.13/ 

t • NO.  OF  PLANT  INPUTS  - *.13/ 

2 • NO.  OF  DISTORBANCE  INPUTS  = ••  13/ 

4 • NO.  OF  PLANT  UOTPOTS  = *.  13/ 

5 • NO.  OF  CONTROLLER  STATES  (EACH  CONTROLLER)  = *.I3) 

WRI TE ( 6 » 9 1 0 ) 

910  FORMAT! • OPLANT  STATE  MATRIX  AP • ) 

110  OO  112  I * 1 • NP 

REA01S.914)  (AP(l.J). J3 1 • NP ) 

112  WRITEI6.913)  (AP1I.J). J* 1 ,NP ) 

913  FORMAT!*  *.8G13.6) 

914  FORMAT ! 4F l 3. 7 ) 

915  FORMAT! 7G13.6 ) 

WRITE  16.916) 

916  FORMAT ( ‘OPLANT  CONTROL  INPUT  MATRIX  — B1P*) 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
inou  COPY  FURNISHED  TO  DDC  


120  OO  122  1 * l.NP 

READ! 5*91 4) 18  l PI 1 . J ) . J= 1 • NUP ) 

122  WRITE  16.913) 1B1P!  I • J ) . J=  1 . NUP  > 

WRITE (6*918) 

918  FORMAT ( • OPL ANT  DISTURBANCE  INPUT  MATRIX B2P • ) 

130  DO  132  1*1. NP 

READ! 5,91 4) <82P( I. J> , J=1 ,NWP» 

132  WRITE!6.913)!82P! I. J). J*1 *NWP) 

WRITE  16*920) 

920  FORMAT  1 'OPLANT  OUTPUT  MATRIX  CP* » 

140  DO  142  1=1, NOP 

READ! 5.91 4>( CP(  I. J>  ,J=1  ,NP) 

142  WRITE  16.913) ! CP ! I • J ) • J= 1 * NP ) 

WRITE  16.922) 

922  FORMAT! '©CONTROLLER  STATE  MATRIX  — FC*) 

150  DO  152  I =1 . NC 

READ!  5. 91  4 MFC!  I.J)  . J=1  .NC  » 

152  WRITE  16.913) ! FC 11.3). J= 1 • NC  ) 

WRI TE 1 6.924) 

924  FORMAT! 'OCONTROLLER  CONTROL  INPUT  MATRIX  GC  *) 

16 a OO  162  1=1 ,NC 

READ! 5,91 4) 1GC! I.J) .J=l.NOP) 

162  WRITE16.913)! GC (I.J). J= 1 • NOP ) 

WRI TE 16.925) 

925  FORMAT! 'OCONTROLLER  OUTPUT  MATRIX  ! STATES)  HC • ) 

170  DO  172  1=1. NUP 

READ! 5.914) 1HC! I,J).J=1.NC> 

172  WRITE  16.913)1 HC ! 1 • J ) • J= 1 • NC ) 

WRITE  16.926) 

926  FORMAT! 'OCONTROLLER  OUTPUT  MATRIX  ! INPUTS)  EC') 

180  DO  182  1=1. NUP 

READ! 5.914) I EC! I.J).J=1 .NOP) 

182  WRI T£!6, 913)1  EC !J.J).J=1, NOP ) 

READ! 5.928) T, NT .NOEL AY .NT AU 
928  FORMAT1F10. 4.515) 

XNT=NT 

XNTAU=NTAU 

XNDEL A=NDELAY 

DELTA=T/XNT 

NIFTAU=NTAU-NDELAY 

NIFT=NT— NDELAY 

WRITE  16.930 )T.NT. NOEL AY .NTAU 
930  FORMA  T! • l T = '.F10.4/ 

1 • NT  * ',15/ 

2 • NDELAY=  » • 1 3// 

3 • NTAU  = '.13// 

4 • T = SAMPLE  RATE.'/ 

5 • NT  = NB.  OF  EVENLY-SPACED  SUBINTERVALS  INTO*/ 

6 • WHICH  T IS  DIVIDED.  •/✓ 

7 • DELTA  = T/NI/NT  = INCREMENT  USED  IN  THE'/ 

0 • NUMERICAL  INTEGRATIONS.'// 

WRI TE 16.931 ) 

931  FORMAT! '©DISTURBANCE  COVARIANCE  MATRIX  W' ) 

184  OO  186  1=1 «NWP 


■ 
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.1 


READC 5,91 4) (M< I , J) , J=l ,NMP) 

186  MR1T£(6.913)(M( I • J ) • J* l •NMP ) 

OO  850  I * l.NMP 

00  850  J * It NMP 
350  M ( I • J ) » M(ltJ)/T 

READ! S.888) I XT1 ME. NXTIME. NXMR IT 
888  FORMAT  < 31 5 ) 

MRITE (6.988)  1 XT  1 ME . NXT I ME • NXMR I T 

988  FORMAT (315) 

N I l*N  I — I 
NHRI TE=NXMR1 T 
NT 1ME=NXT 1ME 

1 T I ME  a I X T I ME 

ccccc  calculate  eccp  and  gccp  ccccccccccccccccccccccccccccc 

OO  T01  I * l.NUP 

DO  701  J * l.NP 

ECCP ( I.J)  a 0.0 
OO  701  K * l.NOP 

701  ECCP(I.J)  * ECCP ( I • J ) ♦ EC( 1 ,K)*CP(K. J) 

OO  702  1 > l.NC 

OO  702  J * l.NP 

GCCP ( I. J)  a 0.0 
OO  702  K a l.NOP 

702  GCCP(I.J)  a GCCP(l.J)  ♦ GC< I .K)*CP(K. J) 
OELMLFaOELTA/2.0 

T 1 *0 

OO  1 1=1. NP 
00  1 J=1.NP 
V l ( I . J)  = 0.0 
PT( 1. J)«0. 0 

1 PS(I.J)=0.0 
I OEL  = 0 

00  2 1=1.NP 

2 Vl(I.l)  = OELHLF 
IMPULS  = 0 
IPRINT=9 

O APPR  X=0 • O 
EPS  a 1 .OE— 7 
NTERMS=6 

1 POLE  =0 
MMAXa— 1 .0 

OO  300  1 * 1,NT 

CALL  MAOO( VI »PT .VI .NP • NPM.NHM. NPM ) 

T 1«TI ♦DELT A 

CALL  EXPK2 ( AP . AH .8 1 P . OAPPRX . OO • DELPHI .E1F.EIV.EPS. 

1 PH1T.FH.FHST.PS1T.  I MPULS • I POLE • I PR  I NT • KMA, NHM • 

2 NRRM.NTERMS.NUP.NUPM.NP.NPM.RF.RR. rr.MMAX. INDEX) 

DO  3 11=1. NP 

OO  3 JJ*1 .NP 

3 PT(I I . J J ) =PMl T(II.JJ) 

CALL  MBVCON ( OELHLF *PT .NP.NHM) 

CALL  MAOOIVl .PT.V1 . NP • NPM . NHM. NPM > 

I OEL« IOEL ♦ 1 

I F ( 1 OEL .EQ.NIFT AO ) GO  TO  4 


p 
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IF  1 I DEL. EQ. NOEL AY  ) CO  TO  7 
IF!  IOEL.EQ.NTAU)  Cl)  TO  10 
IF! IOCL.EQ.NIFT )GO  TO  12 
GO  TO  300 

♦ DO  5 1 1 - 1 . NP 

DO  5 JJ*1.NP 

PHTATDt I I .JJ)*PHITI 1 I.  JJ  ) 

5 VITATOI I I .JJ)  « VIIII.JJ) 

DO  6 1I*1.NP 

DO  6 J J*1 • NUP 

6 P SI  TATI  I I • J J)  *PS1  Ttll.JJ) 

CO  TO  300 

7 DO  a 11*1  . NP 
DO  8 J J* 1 . NP 

PHT  DEL  III .JJ1-PH1T tll.JJ) 

8 VlTOELt I 1 . JJ)=V1 I I I . JJ) 

DO  9 I 1*1 . NP 

DO  9 J J* 1 • NUP 

9 PS1TDE1 I I *JJ) *PS1 Ttll.JJ) 

GO  TO  300 

10  DO  11  I 1*1 .NP 
DO  11  JJ*1»NP 

11  PM1TAUI I I • JJ)*PH1T 11 1 . JJ) 

GO  TO  300 

12  DO  13  11*1. NP 
DO  13  J J* 1 . NP 

P HTTD  tll.JJ)  * PH1TII1.JJ) 

13  V1TTD1I1.JJ)  * VIIII.JJ) 

DO  1A  11*1. NP 

OO  14  J J* 1 • NUP 

14  PS1TTD11I.JJ) *PS1 Ttll.JJ) 

300  CONTINUE 

TAU  * XNTAU4T/XNT 
DELAY  * X NOEL A*  T/XNT 
MR  ITE 16.15)  TAU 

15  FORMAT  1*0* .20X. 'TAU* I NTT AU*T )/NT* • .FIO.S) 

MRITE (6.16)  DELAY 

16  F QRMAT 1 • 0 • .20X, 'DELAY* I XNTTDE4T )/NT* • ,F 10.5) 

C MRITE  PHI  TIT) .PHI  TIT  AU ) .PHI T I T—DEL  AY ) ANDPHI T I ? AU— DEL  AY ) 

MR1TEI6. 171 

17  FORMAT! * OPHl T • ) 

OO  IS  1*1 .NP 

18  MRITE16.91S ) I PHI T 1 1 . J ) • J* 1 . NP ) 

MRITEI6. 19) 

19  FORMAT ('OPHITAU') 

OO  20  1*1 .NP 

20  MR! TE (6.915)1  PH ITAUt  I .J)  • J*l ,NP) 

MR1TEI6.21 ) 

21  FORMAT I • OPHl T I T— DEL AY ) • ) 

DO  22  1*1 .NP 

22  MRITEI6.91S)1 PHT TO I 1 « J ) • J* 1 .NP) 

MRITE I 6.23) 

23  FORMAT  1 '0 PH IT  I TAU— DELAY) • ) 

OO  24  1*1 .NP 
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J34 

CCCCC 


WRlTE(6.915)!PHrAT0( I.J) *Ja 1 . NP ) 

•RITE  PS1 T (T-OELAY ) .PSIT (TAU-DELAY ) AND  PSIT(DELAY) 
WRITEI6.252) 

FORMAT!  »OPSlT  (T  ) • ) 

DO  253  I > I • NP 

WRITE (6,915 >< PSIT ( I. J) . J=I.NUP) 

•RITE (6.25) 

FORMAT! * OP SI T ! T— DELAY ) • ) 

DO  26  ls| ,NP 

■RITE!6 ,9151! PSIT TO! I, Jl .J=l,NUP| 

• R 1TE (6,2  7) 

F ORMAT ( • 0PS1  T (T  AU— DELAY ) • ) 

DO  26  1=1 ,NP 

WRITE  (6.91  5 XPS1TAT!  I,  J)  .J=1.NUP) 

•RITE !6,29) 

FORMAT! • OPS1T! DELAY ) • ) 

DO  30  1=1 .NP 

WRI TE (6.915)1 PSIT DE( I . J) • J=l .NOP) 

CALCULATE  AND  WRITE  PS2T ( TAU-OELAY ) AND  PS2T (DELAY ) 
DO  717  I = 1 • NP 
DO  717  J » 1 • NWP 
PS2TT01I.J)  = 0.0 
DO  717  K * 1 .NP 

PS2TT0! I • J ) = PS2TTD! I , J)  ♦ VI T TD ( i .K ) *B2P . K. J ) 

DO  716  I = l.NP 

DO  718  J = l.NWP 

PS2TOE! I • J ) = 0.0 
DO  718  K = l.NP 

PS2TDE! 1 . J ) = PS2TDE ! I • J)  ♦ Vl TDEL! I .K) *B2P(K. J) 

DO  719  I = l.NP 

DO  719  J = l.NWP 

PS2T  AT! I . J ) = 0.0 
DO  719  K = l.NP 

PS2TAT ( I . J ) = P S2  TAT!  I.J)  ♦ VI T ATO ( I . K ) *B2P( K . J ) 
•RITE16.31 > 

FORMAT!  *0  PS2  T ! TAU-DELAV  ) • ) 

DO  32  1=1 .NP 

WRI  TE (6.915)! PS2T  AT(I.J). Jal . NWP ) 

WR  ITE ( 6.33 ) 

FORMAT! *0  PS2T ( DELAY ) • ) 

DO  34  1=1 .NP 

WRITE16.915)! PS2TDE ( I.J)  »J= l.NWP) 

WRITE (6.333) 

FORMAT! • 0PS2! T-OELAY) • ) 

DO  334  I = l.NP 

WRITE  (6.915  MPS2TTD!  I.J)  • J=1  .NWP) 

CALCULATE  AND  WRITE  G1T.TAU) 

DO  704  I = l.NP 
DO  704  J = l.NWP 
PM(I.J)  * 0.0 
DO  704  K * l.NP 

PM! I.J)  * PM! I.J)  ♦ PHTATD! I .K) 4PS2TDEIK.J) 

DO  3b  1=1 .NP 
DO  36  J=1 .NWP 
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36  AM(I. J1*PM< I.J14PS2TATI I.J) 

C_  FIRST  ROM 

DO  37  1*1 • NP 
DO  37  J*1 .NWP 
G (I • J 1=PS2TTD( I • J ) 

DO  37  K* 1 • NP 

37  Gl  I * J 1*GI  I.  J14PHTTDII  .K 1 4PS2TDEI K . J I 

C SEC UNO  ROW 

DO  38  I* l »NP 
N*NP4 I 

DO  38  J*l # NWP 
36  GIN  .Jl*0.0 

C THIRD  ROW 

DO  39  1*1 .NC 

N*NP4NP4| 

00  39  J=1 . NWP 
39  G I N»  J 1 = 0. 0 

C FOURTH  ROW 

DO  80  1 I * l.NC 
N * NP4NP4NC4I 
DO  801  J **  1 • NWP 

801  GIN.J1  * 0.0 

C FIFTH  ROW 

DO  41  1*1 .NC 

N=NP4NP4NC4NC4l 
OO  41  J=1 .NWP 
GIN. J1*0.0 
00  41  K=1 .NP 

41  GIN. J 1 = 61 N. J» 4GCCPI  I » K ) 4 AM ( K . J ) 

NF  *NP  4NP  4 NC  4NC  4 NC 

WRI TE 16.42) 

42  FORMAT  1*0  GIT.TAU)*) 

DO  43  1*1  .NF 

43  WRITEIb. 91511 Gt I.J) .J*l ,NWP» 

CCCCC  WRITE  PLIT.TAUl 

00  705  I * l.NOP 
DO  705  J * l.NWP 
PLII.Jl  * 0.0 
00  705  K * 1 * NP 

705  PLII.Jl  - PLII.Jl  ♦ fcCCPI 1 .K14AMIK. J! 

WRI TE 16.451 

45  FORMAT! • OPL IT.TAUl  • 1 
DO  46  I * l.NOP 

46  WR1TEI6.91S1I PL  I I • J I . J«I • NWP I 
CCCCC  CALCULATE  ANO  WRITE  FIT.TAUl 

DO  708  I * 1 . NP 
OO  708  J * l.NOP 
0(1  • J 1 * 0.0 
DO  708  K * l.NP 

708  DII.JI  * DII.Jl  4 PHTTDII.KI4PS1TDEIK.J1 
CCCCC  1ST  ROW 

DO  47  I * l.NP 
DO  48  J - l.NP 
FII.J1  * PHI  Til  .Jl 


DO  48  K > 1 .NOP 

48  F 1 1 « J ) - FI  I * J ) ♦PSl T TD 1 1 • K I 4ECCP1  K » J ) 
DO  49  J * l.NP 

M s NP+ J 
F(I.M>  « 0.0 
DO  49  K * l.NUP 

49  F(I.M)  = Ft I.M)+D( I.K)*ECCP(K, J ) 

DO  50  J = l.NC 

M » NP+NP*J 
F(l .M)  a DM(I.J) 

DO  50  K a l.NUP 

50  F(t.M)  a F<  I. Ml  ♦ PS l TTDt I « K ) *HC I K.  J ) 
DO  803  J a l.NC 

M a NP+NP+NC+ J 
FU.M)  a 0.0 

DO  803  K a l.NUP 

803  FII.M)  a F(  I.M)*0(I.K)*HC<K.J> 

DO  52  Ja l ,NC 

M = NP*NP4NC»NC4J 

52  FII.M)  a 0.0 

47  CONTINUE 

CCCCC  2ND  ROW 

DO  53  1 * l.NP 

N a NP+ I 

DO  54  J a l.NP 

54  FIN.J)  a 1. 

DO  55  J a l.NP 
M = NP+J 

55  FIN.M)  a 0.0 

DO  56  J a l.NC 
M a NP+NP+ J 

56  FIN. Ml  a 0.0 

DO  58  J a l.NC 
M a NP4-NP  + NC4-J 

58  FIN.M)  a 0.0 

DO  804  J a l.NC 

M a NP4NP4NC4NC+J 

804  FIN.M!  a o.O 

53  C ONT I NUE 

CCCCC  3RD  ROW 

DO  59  la  l.NC 
N a NP+NP+1 
DO  60  J « l.NP 

60  FIN.JJ  a GCCPII.J) 

DO  61  J a 1 . NP 

M a NP+J 

61  FIN.M)  a o.O 

00  62  J a 1 .NC 
M = NP4NP4J 

62  FIN.M)  a FCII . J) 

DO  64  J a |.NC 

M a NP+NP+NC+ J 
64  FIN.M)  a 0.0 

DO  805  J a l.NC 
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M = NP+NP+NC+NC* J 

805  F IN.M ) * 0.0 
59  CONTINUE 
CCCCC  4 TH  ROW 

OO  806  I « 1 • NC 
N = NP+NP4NC* I 
00  807  J = l.NP 

807  F IN. J ) = 0.0 

OO  808  J = l.NP 

M = NP4J 

808  FIN.M)  = 0.0 

OO  809  J = l.NC 

M * NP+NP+J 

809  FIN. MI  = 1. 

OO  810  J «=  l.NC 

M = NP  + NP  +NC ♦ J 

810  FIN.M)  * 0.0 

OO  811  J = l.NC 

M = NP*NP*NC+NC ♦ J 

811  FIN.M)  = 0.0 

806  CONTINUE 
CCCCC  FIFTH  ROW 

DO  710  1 = l.NP 
OO  710  J = l.NC 
OWII.J)  * 0.0 
OO  710  K = l.NUP 

710  OWII.J)  * OWII.J)  * PS1TATI I . K ) *HC I K . J ) 

00  711  I * l.NP 
DO  71 1 J * l.NUP 

01  I • J ) = 0.0 

OO  711  K « l.NP 

711  OII.J)  = OII.J)  ♦ PHTATDI I . K ) *PS 1 TOE I K. J ) 
00  712  I * l.NP 

OO  712  J = l.NP 
AMII.J)  * 0.0 
OO  712  K = l.NUP 

712  AMII.J)  = AMII.J)  «•  PS 1 T AT  I 1 • K ) 4ECCP I K • J ) 
OO  65  I = l.NP 

00  65  J « l.NP 

65  PTtl.J)  = PHI  TAUI  I • JH-AMt  I • J) 

OO  71 J I = l.NP 
DO  71 J J « l.NC 
AMII.J)  * 0.0 
OO  713  K = l.NUP 

713  AMII.J)  = AMII.J)  ♦ O I 1 • K ) *HC I K • J ) 

OO  714  I * l.NP 

OO  714  J « l.NP 
PMII.J)  = 0.0 
DO  714  K * l.NUP 

714  PMII.J)  s PMII.J)  ♦ Oil  .K)*ECCPIK.J) 

OO  72  Is  l.NC 
N s NP+NP+NC+NC ♦ I 
DO  73  J * 1 .NP 
FIN.  J 1 = 0.0 
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DO  73  K « 1 . NP 

73  FtN.J)  = F(Ni J)*CCCP( I.K)*PTtK< J) 

DO  74  J = t.NP 

M = NPFJ 
FtN.M)  = 0.0 
DO  74  K a t.NP 

74  F(N.M)  a FtN.M)*GCCPt I.K)*PMtK. Jl 

OO  812  J a l . NC 

M * NPFNPFJ 
FtN.M)  a o.O 
DO  812  K = 1 .NP 

812  F(N.M)  a F ( N. H ) ♦ GCCP ( 1 . K » *OW t K . J ) 
OO  75  J a | . NC 
M a NP+NP+NC+  J 
F(N.M)  a 0.0 
DO  75  K a t.NP 

75  F(N.M)  a FtN.M) FGCCPt I.K)tAM(K.J) 

DO  77  J = i.NC 

M a NP+  NPF  NC  ♦ NC  ♦ J 

77  F(N.M)  a FC(I.J) 

72  CONTINUE 

CCCCC  WRITE  F(T.TAU) 

WRITE (6. 78) 

78  F ORMAT ( • OF t T » T AU ) • ) 

DO  79  I a t.NF 

79  WRITE(b.9IS)f F( I. J ) . Ja| .NF ) 

CCCCC  CALCULATE  HI  AND  H2 

DO  80  la  l.NUP 

DO  81  J = 1 ,NP 

81  Hit  I. J)  a ECCPt  I. J) 

DO  82  J a t.NP 

M a NPFJ 

82  Hl( I.M)  a o.O 

DO  83  J a t ,NC 

M a NPFNP+J 

83  Htll.MI  a HC(I.J) 

DO  85  J a | . NC 

M = NP+NP+NCFJ 
85  H 1 ( I . M)  a 0.0 

DO  814  J = 1 • NC 
M a NP+NP+NC+NC+J 
814  Hl(I.M)  a 0.0 

80  CONTINUE 

OO  86  I = 1 .NUP 
DO  87  J = 1 . NP 

H2( I. J)  a 0.0 
DO  87  K a t.NP 

87  H2( I * J)  a H2 ( 1 . J) ♦tCCP t I »K  ) *PT  t K * J) 
DU  88  J a t.NP 

M a NP4J 

H2(  I.M)  a 0.0 

OO  88  K = I .NP 

88  H2U.M)  a H2(  I.M)  + ECCPt 

DO  89  J a t .NC 
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M * NP«-NP*3 
H2 ( I • M ) x 0.0 
DO  89  K x 1 ,NP 

89  H2(|.M)  « H2( I.M)+ECCP( 1 .K)*DW(K. 3) 

00  813  J - l.NC 
M * NP4NP4NC»3 
M2 ( I . M I * 0.0 
OO  813  K « 1 . NP 

813  M2(  1 • M ) x H2( I .Ml ♦ECCP ( 1 .K)4AM(K,3) 

OO  91  J * l.NC 
M = NP*NP*NC*NC*J 

91  M2(  I.M)  * MCI  1. J) 

86  CONTINUE 

CCCCC  CALCULATE  HA.HB  AND  PLB 

00  92  1 * t.NUP 
OO  92  3 * l.NF 

92  HA(1,3)  x HU  J.  J1-M21  I i J» 

00  71S  I X l.NUP 

OQ  715  J * l.NF 
0(1.31  - 0.0 
OO  715  K x 1 , NF 

715  0(1. J)  * 0(1. J1  ♦ HI ( 1 ,K)*F (K, J) 

OO  93  1 * 1 .NUP 

OO  93  J * l.NF 

93  H8(l.3)  » O ( I * J 1 — H2 ( 1 . J 1 

00  716  1 « l.NUP 

OO  716  J x l.NWP 
0(1. J)  x o.O 
00  716  K x l.NF 

716  0(1.31  x 0(1,3)  ♦ HI ( I.K)*G(K,J1 
OO  94  I X l.NUP 

00  94  3 x l.NWP 

94  PLB( 1.3)  x D( 1.31-PL( 1.31 
CCCCC  CALCULATE  PXSS  AND  WRITE 

IT  x 2 

1 MAX  x 30 

CALL  MMMT( G.W.GWG.NF .NWP . NFM.NWPM, NWPM.NFM .0.NFM1 
CALL  MOOCAL ( F. GWG . PXSS. AM. PM. NF • NFM . I MAX • 1 T 1 
WRITE (6 .95) 

95  FORMAT ( • 0 STEADY— STATE  COVARIANCE  OF  STATES*) 

OO  96  1 x l.NF 

96  WRITE(6.915)(PXSS(1 .3) ,3*1 .NF) 

CCCCC  CALCULATE  PEAXSS 

CALL  MMMTIHA.PXSS.HPH.NUP.NF.NUPM. NFM. NFM. NUPM. O, NFM ) 
CALL  MMMT(PLB.W.PLWPL .NUP. NWP. NUPM, NWPM. NWPM. NUPM • 

1 O . NF  Ml 
00  97  1 x l.NUP 

OO  97  3 x l.NUP 

97  PEAXSS( 1,3)  * HPH( I.3)+PLWPL( 1.3) 

WRITE (6.981 

98  FORMAT ( • 0 STEAOV— STATE  COVARIANCE  OF  El*) 

OO  99  1 x l.NUP 

99  WR1TE(6.915)(PEAXSS( 1 • 3 1 • 3» 1 .NUP) 

CCCCC  CALCULATE  AND  WRITE  PEBXSS 


I 
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CALL  MMMT(MM.PXSS.HPM.NUP.NF  .NUPM.  NF  M . NF  M.  NUPM  . D . NFM  ) 
CALL  MMMT  (PLB.  W.PLWPL  .NUP.  NWP.  NUPM.NWPM,  NWPM.NUPM. 

t U.NFM) 

DO  100  I * 1 • NUP 
OU  100  J * I* NUP 

100  PEBXSS(I.J)  ■ HPH(  I . J ) FPLWPL!  I • J > 

WRI TE  <6*101  ) 

101  FORMAT!  *0STEADV— STATE  COVARIANCE  OF  E2M 
00  102  I = I .NUP 

102  WRI  TE  (6.915)!  PfcBXSS  C l.JI.J^l  .NUP) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CCCCC  CALCULATE  X(TIME)  TO  A UNIT-STEP  INPUT  FROM  2ERO 
CCCCC  INITIAL  CONDITIONS  C CCCCCCCCCC CCCCCCC CCCCCCCCC CC CCCC 
CCCCCC CCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
WRITE (6.274)  NT  I ME . NV  RITE 

274  FORMA T< • OUNI T-STEP  TIME  RESPONSE  FOR  T AU  POINT*/ 

1 • NO.  OF  T-INCREMENTEO  POINTS  IN  THE  UNIT-STEP  TIME 

2 RESPONSE  * *.IS/*  TIME  RtSPONSfc  TO  BE  WRITTEN 
J EVERY  • . 15. • *T  SECONDS*) 

TIME  « 0.0 
DO  200  1 « l.NUP 

200  YC2(  I ) * 0.0 

WRITE (6.201  ) T I ME. ( YC2 ( 1 ) . I « l . NUP ) 

201  FORMATt *0T1ME=* .013.6.* YC2( T I ME-T ♦ T AUFOELAY ) = *, 

2 G1J.6I 

DO  202  I * l.NUP 

202  EH  I ) * 0.0 

WRITE (6. 203)1  El ( I ) . I*  1 . NUP  ) 

203  FORMAT! I0X. *El  * *.013.6) 

WRITE  (6.204  ) 

204  F ORMAT ( • 0 * ) 

DO  205  I » l.NF 

205  X(I)  » 0.0 

WRITE(6.206)  T 1 ME  . ( X ! 1 ) • I » 1 . NF  ) 

206  FORMAT ( *0T1ME**.G13.6.*X**.7G13.6./.2SX.7G13.6) 

CCCCC  CALCULATE  YC1 ( T I ME+DELAY ) 

DO  207  I = l.NUP 

Q = 0.0 

DO  208  J * l.NF 

208  O = Q * H1(I.J)*X(J) 

207  YC1 ( I ) * Q 

WRITE (6.209) (VCl(I)*I*t .NUP) 

209  FORMAT! I0X. *VC1 (T1ME4DELAY ) * *.G13.6) 

CCCCC  CALCULATE  E2 

DO  210  1 » l.NUP 

210  C2(I)  « VCI ( I )— VC2 ( I ) 

WRITE (6. 2 I I ) (E2( I ) .1*1 .NUP) 

211  FORMAT! I0X. *E2  « *.013.6) 

CCCCC  CALCULATE  VC2 ( T I MEfT AU+OEL AY ) 

00  21 7 I * l.NUP 
217  YC2 ( l ) * 0.0 

DO  213  1 * l.NUP 
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0*0*0 

00  214  J « I.NF 

214  Q = Q ♦ M2(I.J)*X(J) 

213  VC2 ( 1 ) * YC2 ( 1 ) f 0 

DO  21 5 1 * l.NUP 
Q * 0,0 

DO  216  J * l.NWP 
216  Q = OFPLC  1*  J) 

215  YC2( I 1 * Y C2 ( I)*Q 

WRITE <6.2 18) ( YC2( I ) .1=1 .NUP) 

218  F QRMAT ( 1 OX. * YC2 ( T I ME FT  AUFDELAY ) **.G13.6) 
CCCCC  CALCULATE  El 

OO  219  I * I . NUP 

219  Eld)  = YC1  ( I )-VC2(  I ) 

WRITE <6.203) (El ( 1 > • I - 1 .NUP) 

I WRI T E=0 

THETA  * 0.0 
00  220  JT  = 1 • NT  I ME 
OO  221  I = 1 . NF 
O * 0.0 

OO  222  J * l.NF 

222  0 * QfF  ( I » J ) *X ( J ) 

221  XW< I ) = 0 

OO  223  1*1. NF 
Q * 0.0 

00  224  J * l.NWP 

224  Q = Q FG  ( 1 • J) 

223  X(I)  * X W ( I > + Q 
TINE  * T 1 ME  FT 

CCCCC  CALCULATE  PITCH  ANGLE ( T 1 ME FT/2 .0 ) 

T F£T  A = T HE  TA  FX ( 2 ) * T 
CCCCC  CALCULATE  VCl(TIME) 

OO  225  I * l.NUP 

O * 0.0 

DO  226  J * l.NF 

226  O * 0+Hl(  I.J)*X(J) 

225  YC1 ( I ) * Q 
CCCCC  CALCULATE  E2 

00  227  I * l.NUP 

227  E2(()  * YC1 ( I )— YC2< I ) 

CCCCC  CALCULATE  YC2 < T I MEF TAUFOEL AY ) 

OO  228  I * l.NUP 

228  YC2( I ) * 0.0 

OO  229  I * l.NUP 

Q = 0.0 

OO  23  0 J * l.NF 

230  Q * QFH2 ( I.J)WX(J) 

229  YC2I I ) = YC2 ( I ) FO 

OO  231  1 = l.NUP 

O * 0.0 

DO  232  J * l.NWP 

232  O * a FPL  < 1. J ) 

231  YC2 ( 1 ) * VC2 ( 1 ) F Q 
CCCCC  CALCULATE  El 


CCCCC 


226 

225 

CCCCC 

227 

CCCCC 
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DO  245  I * i.NOP 
245  Eld)  * VC  1(1  1 — VC2  ( 1 ) 

(WRITE  * 1WR1TE+1 

1F( I WR1TE.NE.NWR1 TE  ) GO  TO  220 

(WRITE  * 0 

WRITE <6.206  ) T (ME • (X(I).1=(.  NF ) 

WR I TE ( 6 • 2 09 ) < VC  1 < 1 ) . I * l . NUP > 

WRITE (6.21 1 )( E2< 1 ) . 1=1 .NUP) 

WR1 TE (6.218) (VC2( 1 ) .1=1 .NUP) 

WRITE (6.203) (El ( ( ) • I-t .NUP) 

WR1 TE (6 . 1 24 ) THETA 

124  FORMAT!*  PITCH  ANGLE ( T ( ME*T /2 • 0 ) = *.G13.6) 
22  0 CONTINUE 
STOP 
END 
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' JRG1I  OOPY  FlIWUSHED  TO  11DC 

COMMON  CLEMX.MAXI.MAXJ 

DIMENSION  AP(«.«).QIP<4.2| . B2P(4.2) . CP( 3 • 4 ) ,FC ( 2 . 2 ) • 

1 HC(2.2).EC(2.3).PHIT(10.IO).PS1T(10.10>.F(14.14). 

2 DELPHI  ( 10. 10) .EIP( 10.4 ) • E I V ( tO.4 ) » FH(  10. 10). 

3 KM A ( I O ) . Kf(20 I »RK( < U I • I Di 20 )• ECCP ( 2 • 4 ). GCCP( 2 • 4 ) . 

4 INDEX! 10  ).W(2.2).PS(  lO.I0)«PSlTAT(  10.10). GC ( 2 • J ) • 

s psi r ro( to . to > . phtatoi io.iui .phtdel iio.ioi.ahiio.io). 

6 PNT  TD(1U>10).V|(4»*).V1T  Okl,  I 4.  4 ) ■ VI  T At  D(4 .4)  . 

7 PS2TOE ( 1 0. 10) .PM1  14 .1 4 >.PS2TTD( 10.10).VITT0(4.4). 
d AMI  14.14) . 01 14.2), 01  14. 14) .D«(  14.1 4) .PL 12 .2) . 

9 HI (2.14) . H2( 2. I 4 ). HA (2. 14  >.HR(2. 14 1.PLB12.2). 

A PLWPL(2.2)  .PEAXSS1 2.2)  .PEBXSSI2 .2 ) ,PXSS ( 14 . 1 4 ) . 

B FHST (10. 10). DO IIO).PT(IO.IO) .PSI IDE I 10.10). 

C PHITAU (10.10) • PS2 1AT (10.I0),GWG( 14,14). 

D X! 14 ).XW( 14) .El (2) ,E2(2) ,EA(2) .E 1 I ( 2 ) . E22 ( 2 ) . 

£ VC1  ( 2 ) • YC2  (2  ) • PS  1TAU(10.10),VIT  AU  (4.4)  • HPH(  2.2) 

F PS2TAU1 10. 10) 

CCCCC  PROVIDE  MAXIMA  FOR  CALLED  ARRAYS 
NPM  = 4 


NUPM  = 2 
NWPM  = 2 
NOPM  = 3 
NCM  = 2 

NHM  = 2 4 NPM  ♦ NUPM 
NFM  = 24  NPM  ♦ 3 4 NCM 
NRRM  * 4 4 NPM  *■  4 

cccccccccccccccccccccccccccccccccccccccccccccccccc 

cccccccccccccccccccccccccccccccccccccccccccccccccc 

CCCCC  READ  INPUT  DATA  CCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CCCCCCCCCCCCCCCCCtCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
WRITE (6.2222) 

2222  FORMAT! »1*) 

RE AD ( S. 90  0 ) ID 
900  F ORMAT (20 A4 ) 

WRITE (6.902)  ID 
902  FORMAT  (•  1 • .20A4  ) 

READ!  5 , 906 ) NP. NUP • NWP  » NOP  »NC 
906  FORMAT ( S I 3 ) 

WRITE (6.908)  NP • NUP. NWP. NOP • NC 
908  FORMAT! »ONU.  OF  PLANT  STATES  =*.13/ 

1 • NO.  OF  PLANT  INPUTS  = • .13/ 

2 • NO.  OF  DISTURBANCE  INPUTS  •=  »,  13/ 

4 • NU.  OF  PLANT  OUTPUTS  * ••  13/ 

6 • NO,  OF  CONTROLLER  STATES  (EACH  CONTROLLER)  * *,I3) 

WRI TE( 6,910) 

910  FORMAT (• OPL ANT  STATE  MATRIX  — AP • ) 


1 1 0 

DO  112  I > 

1 

.NP 

READ! S.91 4 

) 

( AP ( I. J)  . 

J»l .NP) 

1 

12 

WRITE (6.91 

3) 

( AP ( 1 • J) 

,J»l,NP 

9 

13 

FORMAT!  • • 

.8613.6) 

9 

1 4 

F ORMAT ( 4F 1 

3. 

7 ) 

9 

IS 

FORMAT!  7G1 

3. 

6 > 

WR1TE(6.91 

6) 

I 
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916  FORMAT  1 • OPLANT  CONTROL  INPUT  MATRIX  — B1P*> 

120  DO  122  I » 1 . NP 

RE AO (5.91 4) (B 1 P! I . J ) . J* 1 .NUP) 

122  *RITE(6.9l3MOlP(I.J).J«l  .NUP) 

WRITE16.918) 

Old  FORMAT (• OPL ANT  DISTURBANCE  INPUT  MATRIX  — B2P * I 
130  OU  132  I * I .NP 

READ (5*914 > (B2PI I • J I .J*1.NNP) 

132  WR|TE16.913)(d2P( I . Jl .3*1  * NWP) 

■ RITE (6.920  ) 

920  FORMAT! • OPLANT  OUTPUT  MATRIX  --  CP*) 

140  DO  142  1*1 .NOP 

READ (5.91  4) ICPI I » J ) . 3*1 ,NP» 

142  WRITE (6«913)(CP( I .J) .3*1 . NP ) 

WRITE (6.922  ) 

922  FORMAT! • OCONTROLLER  STATE  MATRIX  — FC*) 

ISO  DO  1S2  I *1 . NC 

READ ( S . 9 1 4 ) (FCl I.3).3*t .NC  > 

1  52  WR1TE(6.913)(FC(I.3)  »J»1  .NC) 

WRITE  16.924  I 

924  F ORMAT ! • OCONTROLLER  CONTROL  INPUT  MATRIX  — GC  •) 

160  DO  162  1*1. NC 

READ 15,91 4) |GC! I . J ) . J- 1 . NOP > 

162  WRI TE 16.913) I GC 1 1 , J) • J»1 .NOP) 

WRITE  16.925 ) 

925  FORMA T I • OCONTROLLER  OUTPUT  MATRIX  (STATES)  ML • » 

170  DO  172  1*1 .NUP 

RE AO (5, 91 4 ) <HC! 1 . J I .3*1 .NC  1 
172  WHI TE 16.91311 HC 1 I • J ) • J* 1 .NC 1 
WRITE  1 6.926 1 

926  FORMAT  I • OCONTROLLER  OUTPUT  MATRIX  1 INPUTS)  — EC*) 

180  DO  182  1*1 .NUP 

RE  AO  {S.914MECI  1*3  ). 3*1  .NOP  ) 

182  WRI  TE16.913XEC1  I .3)  • J*|  .NOP) 

READ! 5.928) T.NT.NDEL AY. NT AU 
928  FORMATIFIO. 4.515) 

XNT=NT 

XNTAU  * NT  AU 
XNDELA  * NDELAY 
N IFT  * NT-NDELAY 
CELT A*T/XNT 

WRI TE 16.930) T.NT.NDEL AY. NT AU 
930  F ORMAT 1 • 1 T * *.F10.4/ 

1 • NT  = • « 15/ 

2 • NDELAY*  *.13// 

3 • NT  AU  * *.13// 

4 • T * SAMPLE  RATE.*/ 

5 • NT  * N6.  OF  EVENLY-SPACED  SUB1NTERVALS  INTO*/ 

6 * WHICH  T IS  DIVIDED.*// 

7 • DELTA  * T/NI/NT  = INCREMENT  USED  IN  THE*/ 

8 • NUMERICAL  INTEGRAT I QNS. • // 

WRI TE 16.931 ) 

931  FORMAT  I • OU ISTURBANCE  COVARIANCE  MATRIX 

184  DO  186  1*1. NWP 


t 


W*  ) 
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READ  IS. 914) CMC I.J). J=1 ,NMP) 

186  MRI TEC6.913)CMC I. J).J*I.NMP) 

00  850  1 a l.NMP 
OQ  850  J * UN  UP 

850  M C I . J ) * Md.J)/T 

REAO( 5.8881 IX T IME . NXT I ME .NXMR1 T 
888  FORMAT (315) 

MRI TE  <6*988)  I XT 1 ME • NXT I ME • NXMR 1 T 
988  FORMAT (315) 

N 1 1 *N I — 1 
NWRI TE*NXMR1T 
NTIME*NXT IME 

1 TIME* I XT IME 

CCCCC  CALCULATE  ECCP  AND  GCCP  CCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
OO  70  1 I = I.NUP 
OO  701  J * l.NP 

eccpc  i.Ji  * o.o 

DO  70 1 K = 1 . NOP 

701  ECCP(I.J)  * ECCPC l.J)  ♦ ECC I.K)*CP(K, J) 

OO  70 2 I * l.NC 

OO  702  J * l.NP 
GCCPC  I.J)  * 0.0 
OO  702  K = l.NOP 

702  GCCPC  I.J)  = 6C  CPC  I.J)  ♦ GC ( I . K ) *CP( K, J) 

DELHLF  = DELT  A/2  .0 

T I *0 

OO  1 1=1 .NP 

DO  1 J=1.NP 
V 1C  1 • J)=0.0 
PTCU  J>  = 0.0 

1 PSCI.J)  = 0.0 
IOEL=0 

00  2 1=1 .NP 

2 V 1 C 1 . I ) = DEL ML  F 
IMPULS  « O 

1 PR I NT  =9 
OAPPR  X=0. 0 
EPS  = 1.0E-7 
N TERMS* 6 

I POLE =0 
• MAX*— 1 .0 
OO  300  I * 1 . NT 

CALL  MAOOCV1.PT.V1.NP.NPM.NHM.NPM) 

TI*TI ♦DELTA 

CALL  EXPK2C AP.AH.B1P.0APPRX.DD.DELPHI.EIF.EIV.EPS. 

1 PH IT .FH.FHST .PS1 T.  IMPULS. 1P0LE. I PR  I NT .KVA. NHM • 

2 NRRM .NTERMS.NUP.NUPM.NP .NPM.RF . HR . T I . MMAX . INDEX ) 

OO  3 I 1=1 .NP 

00  3 JJ*1 .NP 

3 PTC1I.JJ )=PH1T  C I 1 • J J ) 

CALL  MBYCONCOELHLF.PT.NP.NHM) 

CALL  MAOOC VI .PT .VI . NP • NPM . NHM . NPM ) 

1 DEL* I DEL* I 

I F C 1 DEL . E O*  NOEL  AY ) GO  TO  7 


T 

1 


I 

J 

I 

J 

I 

i 

I 

I 


13 


1 4 
300 


15 

16 

C 

17 

16 

19 

20 

21 

22 
C 

252 

253 
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1F(  IDEL.EQ.NTAU)  GO  TO  10 
I F ( IOEL.EO. Nl FT 1GO  TO  12 
GO  TO  300 
OO  8 1 1-1 .N P 
OO  8 JJ-l.NP 

PHTDEU  II.JJ  1-PH1T  (II.  JJ) 

V 1 TDEL ( I I . JJ)-V1 ( I 1, JJ> 

DO  9 ll-l.NP 

OO  9 JJ-1 . N OP 
P SI  TOE ( II.JJ) — P SIT(II.JJ) 

GO  TO  300 
OO  I 1 I 1*1 .NP 
00  11  JJ-1 .NP 

PHITAUd  l.JJ)  * PHI  T(II.JJ) 

V I T AO (II.JJ)  « VI (II.JJ) 

OO  505  II  * l.NP 

OO  505  JJ  - I .HUP 

PS1TAOC I 1 . JJ)  - PSt  T(II.JJ) 

GO  TO  300 
OO  13  I I — 1 .NP 
OO  13  JJ-l.NP 

PHTTO (II. J J ) - PH  I T ( I l.JJ) 

VlTTO(Il.JJ)  - VKII.JJ) 

OO  14  I I- l.NP 

00  14  J J— 1 .NUP 

PSlTTO( 1 l.J J)-PS1T( II.JJ) 

C ONT  I NUE 
TAU-XNTAU4T/XNT 
OEL  ay  — x no el  a*  t/xn  T 
WRITEI6.15)  TAU 

FORMAT ( *0* .20X . • T AU- ( NTAU*T ) /NT- • .F 1 O. 5 ) 

WRITE (6.16)  DELAY 

FORMAT! *0* . 20X .* DELAY- ( XNDELA* T) /NT- • .FI 0.5) 

WRITE  PH IT(T).PHIT(TAU).PHIT(T  —DELAY ) 

WR1 TE (6.17) 

FORMAT! • OPHI T • ) 

00  16  1-1 .NP 

WRI  TE(6.915)(  PH 1T(I.J). J— 1 • NP ) 

WRITE (6.19) 

FORMAT ( • OPHI TAU* ) 

OO  20  I* l.NP 

WRITE(6.915)( PH IT  AU( l.J) • J— 1 . NP ) 

WR I TE (6.21) 

FORMA T(*0PHIT(T-0EL  AY) • ) 

OO  22  1*1 .NP 

WRITE  (6.91SX  PHTTO  ( l.J)  . J-l  .NP) 

WRITE  PSIT(T-OELAV) .PSIT(TAU)  AND  PS1 T (DELAY ) 

WRI TE (6.252 ) 

FORMAT! *0PS1T( T) • ) 

00  253  I - l.NP 

WRI TE (6.91S)(PS1T(I .J) . J-l. NUP) 

WRITE  (6, 25) 

FORMAT!  *0PS1T  (T-OELAY)*  ) 

00  26  I* l.NP 


* 

r 
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26  WR1  TE  <6.9  ISM  PSITTOt  I . J)  . J*1 ,NUP) 

•RITE (6.27) 

27  F ORM  AT  (*0PS1T(T AU ) • ) 

DO  28  1*1 .NP 

28  WRITE (6.91 5 ) ( PS  1 T AU(  l.J) . J*I.NUP) 

WRITE (6.2V) 

2 V F ORM A T I • OP 6 1 T ( DELAY  ) • ) 

OO  30  1*1 .NP 

JO  WR1 TE (6.915) ( PS1 TOE (l.J).  J* I ,NUP) 

C CALCULATE  AND  WHITE  PS2  T ( T AU)  . PS2  T ( DEL  AY  ) AND  PS?  T ( T-OEL  AY  ) 

OO  717  1 * 1 . NP 
OO  717  J = l.NWP 
PS2TTDO.J)  * 0.0 
DO  717  K * 1 . NP 

717  P S2T  T0(  l.J)  a PS? 1 TD (l.J)  ♦ VITTD(I.K) *B?P ( K » J ) 

00  718  I * 1 * NP 

OO  718  J * l.NWP 
PS2TOEO.J)  * 0.0 
OO  718  K * 1 . NP 

718  PS2TOE (l.J)  * PS?  TOE (l.J)  ♦ VlTO£L(  I . K ) * 82P( K. . J ) 

00  719  l * 1 • NP 

DO  719  J * l.NWP 
PS2  T AU ( l.J)  * 0.0 
OO  719  K * 1 . NP 

719  PS2T  AU(  l.J)  a PS2  T AU(  l.J)  ♦ V 1 T AU  ( 1 .K > *82P ( K. . J ) 

WRITE (6.31 ) 

31  FORMAT! • 0PS2T( TAU) • ) 

OO  32  1=1. NP 

32  WRITE (6.91S)( PS2  T AU( l.J)  • J * 1 . NW  P ) 

WR!TE(6. 33) 

33  F ORMA  T ( • 0 PS2 T ( DEL A V ) • ) 

00  34  1*1 .NP 

34  WR  1TE(6.91S)( PS2TOE ( l.J)  . J*1 .NWP) 

WRITE (6.333) 

333  F ORMA  T( • OPS?( T — OELAV ) • ) 

OO  334  1 = 1 • NP 

334  WR1 TE (6.91S)(  PS2TTD( l.J) • J«1 .NWP) 

CCCCC  CALCULATE  AND  WRITE  G(T.TAU) 

C FIRST  ROW 

OO  37  1*1 .NP 

OO  37  J* l.NWP 

Gd • Jl=PS2TTO( 1 .J) 

00  37  K=1.NP 

37  G(t.J)=G<l.J) +PHT TO ( I . K ) 4PS2TOE ( K , J ) 

C SECOND  ROW 

00  36  1=1 .NP 
N*NP» I 

OO  38  J* l.NWP 

38  G (N. J 1=0.0 

C TH1R0  ROW 

OO  39  1=1. NC 

NsNP+NP*! 

00  39  J*1 .NWP 
G(N. J )*0. O 


il 

II 


39 


I 


r 
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C FOURTH  ROM 

OO  80 » l « l.NC 
N * NP+NP+NOl 
OO  801  J * I . NMP 

801  GIN  . J ) ■ 0.0 

C FIFTH  ROM 

OO  * 1 I«1 *NC 

N = NP4NP4NC4NCF1 
DO  4 1 J*l • NMP 
SIN<JI<0<0 
OO  4 1 K* 1 < NP 

41  G I N < J ) « GIN. J14GCCPI 1 «K)4PS2T AU(K, J) 

NF  * NP4NP+NC4-NONC 

MRITE  (6.42) 

42  FORMAT  1*0  GIT. TAU) • ) 

OO  43  1=1. NF 

43  MRITE |6.915)IG| 1 .J) • J= 1 .NMP ) 

CCCCC  MRITE  PL1T.TAU) 

DO  705  1 = l.NUP 
OO  705  J ~ 1 .NMP 
PLII.J)  « 0.0 
OO  705  K =>  l.NP 

705  PLI1.J1  « PLII.J)  ♦ ECCPI l.K)*PS2TAUlK. J) 
MRITE 16.45) 

45  F ORMAT I • OPL I T . T AU ) • ) 

DO  46  1 = l.NUP 

46  MRITEI6. 915)1  PL  I I . J ) • J* 1 • NMP ) 

CCCCC  CALCULATE  AND  MRITE  FIT.TAU) 

OO  708  I = l.NP 
OO  708  J * 1 • NUP 
OII.J)  = 0.0 
OO  708  K * l.NP 

708  OII.J)  * Oll.J)  ♦ PHTTOI I ,K)*PS1TD€ IK. J ) 
CCCCC  1 ST  ROM 

0047  1 * l.NP 
00  48  J « l.NP 
FI  1. J)  = PHI T 1 1 . J ) 

OO  48  K ■ 1 .NUP 

46  FII.J)  * FI  I. J ) 4PS1 TTO ( I • K ) 4ECCPI K. J) 

OO  49  J * l.NP 
M = NP4J 
FI  l.M ) * 0.0 
DO  49  K « l.NUP 

49  FII.N)  * Ft  l.M) 4D I I .K ) *6CCP I K. J ) 

OO  802  J = l.NC 

M * NP4NPFJ 
FI  l.M  ) - 0.0 
OO  802  K * l.NUP 

802  FII.M)  * Ft  1 . M ) 4PS1 TTO  | 1 . K ) 4HC I K • J) 

DO  50  J > l.NC 

M * NP4NP4NCFJ 
FI  l.M  ) » 0.0 
DO  50  K > 1 *NUP 

50  FII.M)  = Ft  l.M)+DI I .KI4HCIK. J) 


50 


rj  ..... 

- i 
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00  52  J*1 .NC 

M a NP+NP+NC+NC+ J 

mi 

[ 1 52 

F|  1,M>  * 0.0 

C ONT 1 NOE 

ccccc 

2 NO  ROW 

DO  53  1=1 .NP 

N » NPt I 

DO  54  J a 1 • NP 

54 

F (N. J 1=1. 

OO  55  J t | • NP 

M = NP+J 

55 

FIN. Ml  a 0.0 

00  56  J a | ,NC 

M * NPfNP+J 

56 

F(N.M)  a 0.0 

DO  58  J a t . NC 

M = NP+NP+NC+J 

58 

F(N.M)  a o.O 

OO  80  3 J = l.NC 

M a NP+NP+NC+  NC  + J 

803 

FIN.  Ml  = 0.0 

53 

CONTINUE 

CCCCC 

3RD  ROW 

OO  59  I a i.nc 

N a NP+NP+ 1 

OO  60  J =1 .NP 

60 

F I N . J ) a GCCP ( I . J 1 

00  61  J a 1 ,NP 

M = NPt J 

61 

[ ; 

F {N.M  1 a 0.0 

00  62  J a i .NC 

M a NP+NP+J 

62 

FIN. M)  a FC(t.J) 

. 

DO  64  J a l.NC 

M a NPtNPtNCtJ 

64 

F I N.M  ) a o.O 

OO  804  J a l.NC 

M a NP  + NP  +NC  + NC ♦ J 

804 

FIN. M)  a 0.0 

59 

CONT INUE 

CCCCC 

4 TM  ROW 

00  805  I a l.NC 

N a NP+NP+NC+ 1 

00  806  J a l.NP 

306 

FIN.Jl  a o.O 

OO  807  J a l.NP 

M = NP+J 

807 

FIN. Ml  a 0.0 

OO  808  J a l.NC 

M = NP+NPtJ 

808 

FIN.N)  a l. 

00  809  J a |,NC 

M * NP+NP +NC+  J 

809 

FIN. Ml  a o.O 

.1 
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DO  810  J ■ l.NC 
M = NP*NP*NC4NC ♦ J 

810  FCN.M)  « 0.0 
805  CONTINUE 

DO  710  I * l.NP 

DO  710  J * l.NC 

AMCI.J)  * 0.0 
DO  710  K - l.NUP 

710  AMCI.J)  * AH< I.J)  ♦ PSITAUC l.K)*HCCK.J) 

DO  714  1 = l.NP 

DO  714  J * 1.  NP 

PM(I.J)  = 0.0 
DO  714  K ■ l.NUP 

714  PMCI.J)  * PMC l.J)  ♦ PS1 T AU ( 1 • K ) 4ECCP (K. J ) 
OO  72  I * l.NC 
N * NP+NP4NC4NC+1 
DO  73  J * l.NP 
F (N. J )-0.0 
DO  73  K * 1 .NP 

73  FCN.J)  b F < N. J 1 ♦6CCPC 1 » K ) 4PH1 T AUC  K. J 1 
DO  74  J a l.NP 

M = NP4J 
F(N.M)  » C.O 
DO  74  K * l .NP 

74  F(N.M)  * FlN.M)*GCCP( l.K)*PM(K. Jl 

OO  81 1 J « l.NC 

M * NP+NP+J 

811  FCN.M)  * 0.0 
DO  75  J * l.NC 
M b NP4NP4NC4J 
FCN.M)  » 0.0 
DO  75  K « l.NP 

75  FCN.M)  b FCN.M)+GCCP( I.K)*AMCK. J) 

DO  77  J ■ l.NC 

M b NP  + NP+NC  + NC  + J 

77  FCN.M)  = FCC1.J) 

72  CONTINUE 

CCCCC  MRITE  FCT.TAU) 

NRITEC6.78) 

78  FORMATC*OF(T.TAU)« ) 

OO  79  1 * l.NF 

79  MRITEC6.91SJCFC I. , NF ) 

CCCCC  CALCULATE  HI  AND  H2 

OO  80  f b l.NUP 

OO  81  J * l.NP 

81  MIC l.J)  « ECCPCI.J) 

OO  82  J B 1 .NP 

M * NP4J 

82  H1CI.M)  « 0.0 

OO  83  J = l.NC 

M « NP*-NPAJ 

83  HIII.M)  3 HCU.J) 

OO  85  J b l.NC 

M b NP+NP+NC+ J 


I 

I 


85  Hl(I.M)  > 0.0 

OO  012  J = l.NC 

M * NP+NP+NCFNCFJ 
012  HI(l.M)  a 0.0 

90  CONTINUE 

OO  86  I * I *N UP 

OO  87  J * I »NP 

H2(I.J)  * 0.0 

OO  87  K * I #NP 

87  H2(  I.J)  * H2  ( I . J ) 4- EC  CP  ( I • K )*PH  I T AU  ( K . J ) 

OO  88  J » I . NP 

M a NP»J 
H21I.M)  < 0.0 

OO  88  K a l.NP 

88  M2 ( I . M ) a H2( l.M)*ECCP( I.K)*PM(K.J> 

OO  813  J * 1 • NC 

M * NP+NP ♦ J 
813  H2II.M)  > 0.0 

OO  89  J a l.NC 

M a NP*NPFNCFJ 
H2 ( I * M ) a 0.0 
OO  89  K a l.NP 

89  H2II.M)  a H2 ( I » N ) ♦EC CPI 1 • K ) * AM (K . J ) 

OO  91  J a l.NC 

M a NP  + NP  +NC  + NC  + J 

91  H2(I.M)  a MC ( I.J) 

86  CONTINUE 

CCCCC  CALCULATE  HA. HB  AND  PL8 

OO  92  I a i , NUP 

OO  92  J a l.N F 

92  HA(I.J)  a HI  I I . J)-H2( 1. J ) 

00  71 5 I a l . NUP 

00  715  J = 1 . NF 
0(1  .J)  a 0.0 
OO  715  K a l.NF 

715  0(1. J)  a 0(1. J)  ♦ H1(I.K)*F(K.J) 

OO  93  I a 1 .NUP 

OO  93  J a I ,NF 

93  HB(I.J)  a 0(1 • J ) -H2 (I.J) 

00  716  I a l.NUP 

OO  716  J a l.NWP 

0(1. J)  a o.o 

OO  716  K a | . NF 

716  0(1. J)  a O ( I.J)  ♦ HI  ( l.K)*G(K. J) 

DO  94  1 a l.NUP 

00  94  J « 1 .N WP 

94  PL8( I.J)  « Dtl.J)-PL(l.J) 

CCCCC  CALCULATE  PXSS  AND  WRITE 

1 T a 2 
IMAX  a 30 

C ALL  MMM  T (G.V.GWG. NF . NWP. NF  M. NWPM. NWPM.NFM • O. NFM 1 
CALL  MODCAL (F .GWG.PXSS. AM.PM.NF .NFM. IMAX. IT) 

WR  ITE  (6.95  ) 

FORMAT  ( *OSTEADY~STATE  COVARIANCE  OF  STATES*) 


95 
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DO  96  I « 1 »NF 

96  MR I TE16.915)! PXSS ! I . J ) . J* 1 » NF ) 

CCCCC  CALCULATE  PEAXSS 

CALL  MMMT (HA.  PXSS.HPH.NUP.Nr .NUPM. NFM.NFM.NUPM. D.NFM) 
CALL  MMMT  ( PLB  • W • PL  WPL  . NUP  . NMP*  NUPM.  NWPM.  NWPM.  NUPM. 
t D.NFM) 

OO  97  1*1 .NUP 
DO  97  J = 1 .NUP 

97  PEAXSS!  I . J ) * HPHC  1 . J)  FPLWPL!  I . J) 

WRITE  16*96) 

96  FORMAT! • OSTEADV-STATE  COVARIANCE  OF  El* I 
DO  99  1*1 .NUP 

99  WRITE <6,9  IS) (PEAXSS! I , J) ,J=l .NUP) 

CCCCC  CALCULATE  AND  WRITE  PEBXSS 

CALL  MMMT  <HB. PXSS. HPH. NUP. NF. NUPM. NFM. NFM, NUPM . D . NFM I 
CALL  MMMT! PLB • W. PL WPL .NUP. NWP, NUPM. NWPM. NWPM. NUPM. 

I  D.NFM) 

OO  100  I = l.NUP 
00  100  J = l.NUP 

100  PEBXSS!  I • J)  * HPH!  I.  J) FPL WPL!  I , J ) 

WRI TE (6. 101) 

101  FORMAT! • OSTEADV-STATE  COVARIANCE  OF  E2*) 

DO  102  I = l.NUP 

102  WRITE16.91S)! PEBXSS! I. J) ,J=1.NUP) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CCCCC  CALCULATE  X(T1ME)  TO  A UNIT-STEP  INPUT  FROM  ZERO 
CCCCC  INITIAL  CONDITIONS  CCCCCCCCCCCCCCCCCCCCCCCCGCCCCCCCC 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
WRITEI6.27A)  NT  I ME .N WRITE 

274  FORMAT! • OUNIT-STEP  TIME  RESPONSE  FOR  TAU  POINT*/ 

1 • NO.  OF  T-INCREMENTEO  POINTS  IN  THE  UNIT-STEP  TIME 

2 RESPONSE  * *.I 5/*  TIME  RESPONSE  TO  BE  WRITTEN 

3 EVERY  *.!S.**T  SECONDS') 

TIME  = 0.0 

DO  200  I «=  l.NUP 

200  YC2! I ) * 0.0 

WRITE (6 .201)  T I ME • ( VC2 ( I ) .1=1 .NUP) 

201  FORMAT! *OTJME  = * .G13.6. • YC2! TIME— TFT  AUFDELAY)  = *. 

2  613.6) 

DO  202  I = l.NUP 

202  El(I)  » 0.0 

WRITE !6, 203) ! El ( I ) .1=1. NUP ) 

203  FORMAT! 10X. *E 1 « *.613. 6) 

WRITE (6.204) 

204  FORMAT! *0*) 

DO  20S  I = I.NF 

205  X(I)  = 0.0 

WR ITE ! 6. 206 ) T 1ME . ( X( I ) . 1= 1 .NF ) 

206  FORMAT! *0T1ME  = * • G1 3 .6. * X* • • 7G 1 3.6. /.25X. 7G13.6) 

CCCCC  CALCULATE  YC 1 ( T I ME* DELAY) 

DO  207  I = l.NUP 

Q = 0.0 
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OO  208  J = l.NF 

208  Q » Q *■  H1(1»J)*X(J) 

20  7 YCIC I ) * Q 

WHITE (6.209) ( YC I ( I ) • Is  I *NUP ) 

209  F ORMAT  (iOXt'YCl  (TI  ME  ♦DEL  AY  ) = *.G13.6) 
CCCCC  CALCULATE  E2 

00  210  I a X.NUP 

210  E2(I)  = YC  1(1  )~VC2(  I ) 

WRITE (6.21 1 )1E2( 1 ). 1=1. NOP) 

211  F ORMAT (10X«'£2  a *.G13.6> 

CCCCC  CALCULATE  YC2 ( T IME 4TAU40E LAV > 

DO  217  I a 1 ,NUP 

217  YC211)  a o.O 

00213  I a i , NUP 

Q * 0,0 

OO  214  J a l.NF 

214  Q = Q ♦ H2 ( I • J ) 4X1 J ) 

213  YC2( I > a YC21 1)40 

OO  215  I a l.NOP 

Q a 0.0 

OO  216  J a t.NWP 
216  Q * Q4PL1 I . J) 

215  YC21 I ) a VC21 II  ♦ Q 

WRITE (6,21 8)1 YC21  I).I*1  .NUP) 

218  FORMAT! 10X. • YC2 ( T I ME4T AU4DELAY > a*,G13.6> 
CCCCC  CALCULATE  El 

00  219  I a I.  NUP 

219  EMI)  a YCt  ( I ) -VC2 11) 

WRITE (6.203)1 £1! ll.Iat.NUP) 

1 WRI T E=0 
THETA  = 0.0 

OO  220  JT  a 1 .NT IME 
OO  221  1 a l.NF 

Q a 0.0 

OO  22  2 J a l.NF 

222  O a QFF( I . J)4X< J) 

221  XW( I ) a Q 

OO  223  I * 1 . NF 
Q a 0.0 

OO  224  J a l.NWP 

224  Q a QFG! I . J) 

223  X ( 1 1 a x»(  1 )4Q 
TIME  a TIME FT 

CCCCC  CALCULATE  PITCH  ANGLE ( T I ME4T/2 . 0 1 
THETA  = THETA  + X(2>*T 
CCCCC  CALCULATE  YC11TIME) 

OO  225  I = I .NUP 
Q a 0.0 

OO  226  J a 1 , NF 
226  Q a QFH1 { I. J)*X( J) 

225  YCt ( I ) a Q 
CCCCC  CALCULATE  E2 

OO  227  I a 1 . NUP 
E2C1I  a YC1 ( I 1-YC2! I ) 


227 
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ccccc  calculate  yc2i t ime*tao*oelay> 

OO  228  I * 1 * NOP 

228  YC2( I > - 0.0 

OO  229  I ■ I .NUP 

Q * 0.0 

OO  230  J * l.NF 

230  0 • Q+H2 ( I.J)4X<J) 

229  Y C2{ I ) » YC2 ( I ) ♦Q 
OO  231  I = l.NUP 
Q * o.O 

OO  232  J * l.NUP 
232  Q * Q*PL< l,Jt 

231  YC2< 1 1 * YC2( 1 > ♦ Q 
CCCCC  CALCULATE  El 

00  2A » I a l.NUP 

245  El(l)  « YCld)-VC2d> 

1 WRITE  « I WR I TE4 1 
IFI1WR1TE.NE.NWRITE)  GO  TO  220 
1 WRITE  ■ 0 

WRITE(6.206)  TIME.IXd 1.1-1 ,NF> 

WRITE (6.209)  ( YC1 ( 1 > • 1*1 .NUP) 

WRI TE 16.21 1 1IE2II ) . I * 1 . NOP ) 

WRITE  16.21  8)  I YC2d  > .1*1  .NUP) 

WRITE  (6.203  MEM  I ).  !»  l.NUP) 

WRITE16.I24)  THETA 

124  FORMAT!*  PITCH  ANGLE ( T IME4T/2. 0)  * *.013.6) 

220  CONTINUE 
STOP 
END 
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OF  OI  STUR  JANCF  INPUTS  = 1 
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COUPLING  ACTIVITIES 

This  research  originated  under  an  AFOSR  program,  namely,  the  1975 
USAF-ASEE  Summer  Faculty  Research  Program  (June  9 - Aug.  15).  Professor 
Slivinsky  was  one  of  22  participants  and  was  assigned  to  AFFDL/FGL  (DAIS). 
He  was  assisted  by  Major  Vincent  J.  Darcy.  Together  they  devised  a set 
of  differential  equations  to  model  a dual -redundant,  closed-loop  flight 
control  system  with  a simple  voting  algorithm.  A software  package  was 
written  to  allow  parametric  analyses  of  the  effects  of  design  parameters 
on  both  transient  and  steady-state  inherent  errors. 

AFOSR  sponsored  the  continuation  of  this  research  with  a grant  to 
University  of  Missouri -Columbia  (AF0SR-76-2968,  1 Feb.  76  - 31  Jan.  77). 
Professor  Slivinsky  (at  UMC)  and  Major  Darcy  (at  AFFDL)  verified  the  model 
and  software  and  applied  both  to  a version  of  the  A-7D  flight  control 
system.  They  found  that  inherent  errors  were  on  the  order  of  3%  of 
command  inputs.  While  small,  these  errors  are  not  negligible,  as  their 
characteristics  must  be  known  to  specify  the  algorithms  that  isolate 
failed  signals  and  determine  the  best  output  from  among  the  unfailed, 
redundant  signals.  This  work  resulted  in  an  AFFDL  Technical  Report 
(AFFDL-TR-76-16  dated  April  1976)  and  a paper  in  the  IEEE  Transactions 
on  Automatic  Control  (February  1977). 

AFOSR  renewed  the  grant  for  a second  year  (AF0SR-76-2968A,  1 Feb.  77  - 
31  Jan.  78),  to  develop  more  powerful  closed-loop  models  and  software, 
and  to  study  signal-selection  algorithms.  The  last  five  months  of  this 
grant  overlap  with  the  1977-78  Sabbatical  Leave  of  Professor  Slivinsky 
at  AFFDL. 
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For  the  last  two  years,  the  Flight  Dynamics  Laboratory  supported 
work  related  to  the  AFOSR-sponsored  research.  During  the  Summer  of 
1976,  Professor  SI i vinsky  spent  six  weeks  as  a Consultant  (University 
of  Dayton  Contract  F33615-76-C-3076)  at  AFFDL/FGL  working  under  Dr.  A. 
DeThomas.  The  major  result  was  the  report  "Some  Guidelines  for  Testing 
AFFDL  DAIS  Concepts,"  which  was  used  in  writing  the  DAIS  Concepts  Test  Plan. 

For  the  period  Nov.  76  - July  77,  Professor  Slivinsky  continued  as 
a Consultant,  working  approximately  one  day  per  week.  This  work  was 
mainly  to  follow  the  progress  of  the  DAIS  program,  offer  consultation, 
and  special  studies. 

Since  Sept.  77,  Professor  Slivinsky  has  been  on  Sabbatical  Leave  at 
AFFDL/FGL  (DAIS).  For  the  period  Sept.  77  - Uan.  78,  he  spends  one  day 
per  week  on  the  AFOSR-sponsored  research  and  the  remainder  on  AFFDL 
related  research  and  development.  His  specific  tasks  include  writing 
the  flight  control  software  for  the  AFFDL  DAIS  Flight  Engineering 
Facility  and  studying  problems  associated  with  using  higher  order 
languages  in  digital  flight  control.  AFFDL  provides  half-time  support 
during  this  period. 


